Nanopore functionality control method and system

Abstract

The invention relates to the technical field of nanopore detection, in particular to a nanopore functionality control method and system.The nanopore functionality control method comprises the following steps of extracting path information based on liquid bridge anchor point distribution, constructing a dynamic configuration state, generating an interface response map, recognizing a stability signal, correcting confinement field distribution and monitoring a temperature-sensitive response behavior. And forming an on-off state transition trigger mark. According to the method, by introducing sequence control and staggered behavior recognition of a liquid migration path, optimized organization of a liquid bridge distribution dynamic configuration is achieved, in combination with interface response mapping of a displacement peak index, the precision and efficiency of liquid position regulation and control are effectively improved, and by means of path and boundary response convergence extraction and periodic direction change analysis, the accuracy and efficiency of liquid position regulation and control are improved. The dynamic pre-judgment of the stable state of the boundary of the nano orifice is realized, and the stability and controllability under the scene of molecular on-off regulation and structural dynamic recognition are improved through liquid coverage area expansion and tension path compression.

Description

Technical Field

[0001] The present invention relates to the field of nanopore detection technology, and in particular to a nanopore functionality control method and system. Background Art

[0002] The field of nanopore detection technology involves the use of nanopores to precisely identify and analyze substances. Core issues include the fabrication and functionalization of nanopores, as well as molecular-level detection and detection through nanopores. Nanopore detection technology is widely used in fields such as molecular biology, chemical analysis, environmental monitoring, and medical diagnostics. In molecular biology, nanopores enable single-molecule-level measurement and analysis of molecules such as DNA, RNA, and proteins. In nanocarrier applications, nanopore technology is used for molecular detection and carrier control, enabling targeted therapy and delivery of nanoparticle-based drug carriers. In environmental monitoring, nanopore technology can be used to detect trace pollutants in air and water. In medical diagnostics, the precise detection capabilities of nanopores can be used for early disease screening and genetic analysis. Nanopore functional control methods involve optimizing and controlling the performance of nanopores by manipulating their geometry, material composition, or electrical environment. These methods employ methods such as adjusting the size, shape, surface modification, and electric field of the nanopore to meet specific application requirements. In the field of molecular biology, the permeability of specific molecules or ions can be regulated by controlling the size of nanopores, thereby improving the sensitivity and accuracy of molecular detection. The function of nanopores can be controlled by chemical modification, addition of surfactants, electric field regulation, etc., to ensure accuracy and efficiency in applications.

[0003] Although existing technologies can achieve the screening and detection of molecules by controlling the size of nanopores and adjusting the electric field environment, they lack the identification mechanism for interlaced behavior during dynamic liquid migration and interface morphology response, and cannot achieve real-time regulation of the interaction between complex liquid bridge structures. They do not have the ability to judge periodic trends in the continuity analysis of interface responses, which can easily lead to unstable changes in boundary states. When it comes to the response of temperature-sensitive materials, they fail to establish a coordinated judgment standard between structural thickness and response ratio, resulting in an increase in the deviation in on-off state recognition. In actual applications, in temperature-controlled drive and high-sensitivity detection scenarios, problems such as judgment lag or response interruption are prone to occur, affecting the repeatability and accuracy of the overall performance of the nanopore. Summary of the Invention

[0004] In order to solve the problems existing in the prior art, such as the lack of an identification mechanism for interlaced behavior in the process of dynamic liquid migration and interface morphology response, the inability to achieve real-time regulation of the interaction between complex liquid bridge structures, the lack of periodic trend judgment capability in the continuity analysis of interface response, which easily leads to unstable changes in the boundary state, and the failure to establish a coordinated judgment standard between the structure thickness and the response ratio when involving the response of temperature-sensitive materials, resulting in an increase in the deviation in the recognition of the on-off state. In actual applications, judgment lags or response interruptions are easily generated in temperature-controlled drive and high-sensitivity detection scenarios, affecting the repeatability and accuracy of the overall performance of the nanopore, the present invention provides a nanopore functional control method and system. The technical solution is as follows: In one aspect, a method for controlling nanopore functionality is provided, the method comprising: S1: Based on the initial distribution state of nanopore liquid bridge anchor points in the porous array, the liquid volume, contact interface length and adhesion migration limit are collected to screen the staggered behavior segments in the migration path and generate the initial dynamic configuration state of the liquid bridge network; S2: extracting the projection trajectory from the edge of the liquid bridge to the nanopore according to the initial dynamic configuration state of the liquid bridge network, identifying the displacement change of the end point of the trajectory, performing positional translation on the liquid distribution, and generating a nanopore interface morphology response shift map; S3: calling the nanopore interface morphology response deviation map, extracting the pore edge curvature sequence, the liquid bridge path length change amplitude and the molecular contact frequency fluctuation range, performing continuity detection on the curvature change trend and extracting the periodic turning points, performing projection sequence sorting on the path length change amplitude, identifying the segment that converges with the curvature change, and generating the nanopore boundary stability response signal; S4: According to the nanopore boundary stability response signal, the liquid coverage area and tension path value of the liquid bridge anchor point in the instability zone are extracted, the coverage area is amplified according to the set expansion ratio, and the conduction path is compressed to the constraint interval to generate a limited correction field distribution map of the boundary area.

[0005] As a further solution of the present invention, the initial dynamic configuration state of the liquid bridge network includes the anchor network topology, the path priority connection sequence, and the staggered path distribution pattern; the nanopore interface morphology response offset map includes the interface response offset contour, the displacement dynamic amplitude range, and the trajectory terminal offset mapping result; the nanopore boundary stability response signal includes the stability period identification, the path trend coincidence section indication, and the contact frequency fluctuation amplitude parameter; the boundary area confined correction field distribution map includes the tension compression mapping map, the coverage area expansion area, and the confined field strength change block.

[0006] As a further solution of the present invention, the steps for obtaining the initial dynamic configuration state of the liquid bridge network are specifically as follows: S101: Based on the initial distribution state of the nanopore liquid bridge anchor points in the porous array, the liquid volume and the contact interface length in the nanopore liquid bridge are collected, and a corresponding mapping is established between the liquid volume formed between the anchor points and the contact interface length. The liquid attachment migration limit between each pair of anchor points is recorded to generate a liquid bridge characteristic value table; S102: Recalling the liquid adhesion migration limit and contact interface length data between anchor point pairs in the liquid bridge characteristic value table, calculating the corresponding projection trajectory length based on the two-dimensional plane projection coordinates between adjacent nanopore liquid bridge anchor points, sorting the anchor point pairs in ascending order, and removing paths that do not meet the migration conditions in the sorted sequence based on the liquid adhesion migration limit between the anchor point pairs to obtain a sequence of migratable path lengths; S103: Based on the order and position distribution of the path starting points in the migratable path length sequence, the path starting points are controlled in sequence and the interlaced behavior state of the path migration is tracked, the position segments where the cross-overlapping sections appear in the path are identified and marked, and the initial dynamic configuration state of the liquid bridge network is generated.

[0007] As a further solution of the present invention, the steps for obtaining the nanopore interface morphology response shift map are specifically as follows: S201: extracting a two-dimensional projection trajectory from the edge of the liquid bridge to the center of the corresponding nano-orifice according to the initial dynamic configuration state of the liquid bridge network, recording a path coordinate sequence of the projection trajectory in the initial coordinate system, and uniformly mapping the trajectory to a unified coordinate system to generate a projection path sequence from the liquid bridge to the orifice; S202: Based on the liquid bridge to orifice projection path sequence, identifying the coordinate value of the trajectory endpoint on the nano-orifice surface, comparing the endpoint coordinates of multiple time sections between the same anchor points, calculating the trajectory endpoint displacement equilibrium characteristic value, extracting the peak value change and recording the coordinate index position thereof, and generating the trajectory endpoint displacement index distribution value; The displacement balance characteristic value of the trajectory endpoint is calculated using the formula: ; in, Represents anchor pair 、 The displacement equilibrium eigenvalue of the trajectory endpoint between Represents an anchor point In the time section The horizontal end point coordinate value of Represents an anchor point In the time section The horizontal end point coordinate value of Represents an anchor point In the time section The vertical end point coordinate value, Represents an anchor point In the time section The vertical end point coordinate value, Represents an anchor point In the time section The local displacement increment, Represents an anchor point 、 The average value of the displacement increment in the local direction in all time sections, represents the number of time sections, The sequence index value representing the time section, 、 Represents the selected anchor point number; S203: Mapping the spatial distribution area of ​​the original liquid around the nanopore according to the corresponding index position in the displacement index distribution value of the trajectory endpoint, uniformly translating the liquid distribution area in a two-dimensional coordinate system, determining the translation direction according to the direction indicated by the displacement index, and determining the translation amplitude using the displacement increment to generate a liquid distribution translation coordinate set; S204: calling the liquid distribution translation coordinate set, reorganizing the liquid distribution area according to the corresponding nanopore number, and reconstructing the interface contour graphic under the standard pore contour boundary, plotting the boundary contour area corresponding to the liquid distribution difference before and after the displacement, and obtaining the nanopore interface morphology response offset map.

[0008] As a further solution of the present invention, the step of obtaining the nanopore boundary stability response signal is specifically as follows: S301: calling the nanopore interface morphology response offset map, extracting the continuous numerical value of the pore edge curvature along the boundary coordinate direction, collecting the length change value of the corresponding liquid bridge path and the molecular contact frequency sequence of the liquid bridge end contact area, recording the numerical fluctuation interval on the time axis, and generating a curvature path frequency joint series; S302: Based on the hole edge curvature sequence in the curvature path frequency joint series, the curvature value change relationship between consecutive coordinate points is detected in the order of boundary coordinates. For the positions where the adjacent positive and negative incremental directions switch, the periodic direction change points are extracted, and the index sections where the curvature trend reverses are marked to obtain the edge curvature direction change index set; S303: Call the liquid bridge path length change amplitude in the edge curvature direction change index set, rearrange the coordinates and normalize the direction of the path change amplitude according to the index segment, select the path segment with the same change trend as the direction change segment as the convergence segment, count the difference range of the molecular contact frequency fluctuation in the corresponding segment, and obtain the nanopore boundary stability response signal.

[0009] As a further solution of the present invention, the step of obtaining the boundary region limited correction field distribution map is specifically as follows: S401: extracting the liquid bridge anchor point number marked as the unstable section according to the nano-orifice boundary stability response signal, recording the liquid coverage area under the corresponding anchor point and the length of the tension conduction path to the orifice boundary, and generating a liquid bridge anchor point coverage tension value group; S402: Call the liquid bridge anchor point coverage tension value group, perform numerical amplification processing on the coverage area according to the set expansion ratio, calculate the amplified area characteristic value, and at the same time converge the tension path length to the set tension conduction constraint range, uniformly map it to the two-dimensional nanopore boundary coordinate system, re-plot the tension distribution area contour of the anchor point position in the boundary area, and obtain the boundary area confined correction field distribution map.

[0010] As a further solution of the present invention, the characteristic value of the amplified area is calculated using the formula: ; in, represents the characteristic value of the amplified area, represents the initial area covered by the foundation liquid bridge anchor point, Representative Anchor point tension transmission gain factor, Representative Item anchor point tension value, Representative Item Anchor point path length compression damping coefficient, Representative Item anchor path length, Represents the equivalent radius length formed by mapping the tension path in the two-dimensional plane, represents the average conducted disturbance displacement correction value, Represents the number of anchor points.

[0011] As a further solution of the present invention, the method further includes step S5: S5: calling the boundary area limited domain correction field distribution map, detecting the temperature change rate of the covering thermosensitive composite material and the material response window interval value, extracting the pore edge membrane thickness and the pore mouth variable area ratio data, comparing the execution rate response hysteresis and the periodic continuity between the thermal boundary deformation time, judging whether the continuous response behavior of the thermal trigger behavior has met the closing condition, and generating the nanopore on-off state transition trigger mark; The nanopore on-off state transition trigger mark includes a trigger threshold index, a response closure cycle, and an on-off behavior identification site.

[0012] As a further solution of the present invention, the step of obtaining the trigger mark of the nanopore on-off state transition is specifically: S501: calling the boundary region limited domain correction field distribution map, monitoring the time series temperature change curve of the area covering the thermosensitive composite material, calculating the temperature change rate of the area per unit time, synchronously extracting the window interval value, arranging the temperature rate and response window value pairs according to the coordinate number, and generating a thermosensitive response rate window pairing table; S502: Based on the thermal response rate window pairing table, extract the film thickness data at the nanopore edge position and the variable area ratio data in the pore mouth area, identify the corresponding structure between the pore edge structure size and the response ability, compare the rate response lag time with the duration of thermal boundary deformation, identify the periodic continuous change section in time series order, and obtain the response period continuous section index sequence; S503: Call the response cycle continuous segment index sequence, identify whether the rate response lag in the cycle segment completes the return judgment within the material response window interval, screen the cycle segments that meet the synchronous closure of the deformation duration and the response lag time, mark the corresponding cycle closure information by number, and generate a nanopore on-off state transition trigger mark.

[0013] In another aspect, the nanopore functionality control system is used to perform the above-mentioned nanopore functionality control method, and the system comprises: The path recognition module obtains the coordinate values ​​of the nanopore liquid bridge anchor points in the porous array, the liquid volume value, the contact interface length value, and the attachment migration limit value. It calculates the sorting sequence based on the projected trajectory length between the liquid bridge anchor points. Based on the sorting sequence, it controls the arrangement order of the starting points of the liquid migration path to generate the initial dynamic configuration state of the liquid bridge network. The interface response module extracts the two-dimensional projection trajectory data from the edge of the liquid bridge to the nanopore according to the initial dynamic configuration state of the liquid bridge network, determines the change direction between the trajectory end point and the displacement vector of the orifice area, and generates a nanopore interface morphology response offset map; The structural stability module calculates the position of the continuity inflection point within the curvature change range according to the nanopore interface morphology response deviation map, selects the boundary dividing point between the structural stability range and the instability range, and obtains the nanopore boundary stability response signal; The confinement correction module calls the nanopore boundary stability response signal, extracts the liquid coverage area value and the tension path length value of the liquid bridge anchor point in the instability zone, adjusts the coverage area range according to a fixed expansion ratio, and generates a confinement correction field distribution map of the boundary area; The thermal trigger module detects the temperature change rate value and the material response time window value of the thermosensitive composite material based on the limited correction field distribution map of the boundary area, calculates the period ratio between the temperature change lag time and the pore deformation time, and obtains the trigger mark of the nanopore on-off state transition.

[0014] The beneficial effects brought about by the technical solution provided by the embodiment of the present invention include at least: By introducing sequential control of liquid migration paths and identification of staggered behavior, the optimized organization of the dynamic configuration of liquid bridge distribution is achieved. Combined with the interface response mapping of the displacement peak index, the accuracy and efficiency of liquid position control are effectively improved. Through the convergence extraction and periodic direction change analysis of the path and boundary response, the dynamic prediction of the stable state of the nanopore boundary is achieved. By expanding the liquid coverage area and compressing the tension path, an adjustable confined correction field is constructed and the local response enhancement of the tension concentration area is achieved. The thermally triggered closing behavior judgment under temperature-sensitive conditions is realized in the changes of the nanopore edge membrane thickness and area ratio. The precise marking and response cycle control of the nanopore on-off state are completed, which significantly improves the stability and controllability in the scenarios of molecular on-off regulation and structural dynamic recognition. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] Figure 1 It is a schematic diagram of the workflow of the present invention; Figure 2 It is a system flow chart of the present invention. DETAILED DESCRIPTION

[0016] The technical solution of the present invention is described below in conjunction with the accompanying drawings.

[0017] In the embodiments of the present invention, words such as "exemplarily" and "for example" are used to indicate examples, illustrations, or explanations. Any embodiment or design described as an "exemplary" in the present invention should not be interpreted as being preferred or advantageous over other embodiments or designs. Rather, the use of the word "exemplary" is intended to present concepts in a concrete manner. Furthermore, in the embodiments of the present invention, "and / or" can mean both or either of the two.

[0018] In order to make the technical problems, technical solutions and advantages to be solved by the present invention clearer, a detailed description will be given below with reference to the accompanying drawings and specific embodiments.

[0019] See also Figure 1 The embodiment of the present invention provides a method for controlling the functionality of a nanopore. The processing flow of the method may include the following steps: S1: Based on the initial distribution state of nanopore liquid bridge anchor points in the porous array, the liquid volume, contact interface length and adhesion migration limit are collected, the projected trajectory lengths between adjacent nanopore liquid bridge anchor points are called and sorted, the starting order of the liquid migration path is controlled, the interleaving behavior segments in the migration path are screened, and the initial dynamic configuration state of the liquid bridge network is generated; S2: Based on the initial dynamic configuration state of the liquid bridge network, the projection trajectory from the liquid bridge edge to the nanopore is extracted, the displacement change at the end point of the trajectory is identified, the peak displacement position index of the corresponding pore interface is calculated, the liquid distribution is positionally shifted, and a nanopore interface morphological response shift map is generated; S3: Call the nanopore interface morphology response deviation map to extract the pore edge curvature sequence, the liquid bridge path length change amplitude and the molecular contact frequency fluctuation range, perform continuity detection on the curvature change trend and extract the periodic turning points, perform projection sequence sorting on the path length change amplitude, identify the segment that converges with the curvature change, and generate the nanopore boundary stability response signal; S4: According to the nanopore boundary stability response signal, the liquid coverage area and tension path value of the liquid bridge anchor point in the instability zone are extracted, the coverage area is amplified according to the set expansion ratio, and the conduction path is compressed to the constraint interval to generate a confined correction field distribution map of the boundary area; S5: Call the boundary area limited correction field distribution map, detect the temperature change rate of the thermosensitive composite material and the material response window interval value, extract the pore edge membrane thickness and the pore mouth variable area ratio data, compare the execution rate response hysteresis and the periodic continuity between the thermal boundary deformation time, determine whether the continuous response behavior of the thermal trigger behavior has completed the closing condition, and generate the nanopore on-off state transition trigger mark; The initial dynamic configuration state of the liquid bridge network includes the anchor network topology, path priority connection sequence, and staggered path distribution pattern. The nanopore interface morphology response offset map includes the interface response offset contour, displacement dynamic amplitude range, and trajectory terminal offset mapping results. The nanopore boundary stability response signal includes the stability period identification, path trend matching section indication, and contact frequency fluctuation amplitude parameter. The boundary area confined correction field distribution map includes the tension compression mapping map, coverage area expansion area, and confined field strength change block. The nanopore on-off state transition trigger mark includes the trigger threshold index, response closure cycle, and on-off behavior identification site.

[0020] The specific steps for obtaining the initial dynamic configuration state of the liquid bridge network are: S101: Based on the initial distribution state of the nanopore liquid bridge anchor points in the porous array, the liquid volume and the contact interface length in the nanopore liquid bridge are collected, and a corresponding mapping is established between the liquid volume formed between the anchor points and the contact interface length. The liquid attachment migration limit between each pair of anchor points is recorded to generate a liquid bridge characteristic value table; A coordinate grid system for the entire array structure is established to identify the position of each nanopore and mark and number them. On this basis, for each pair of anchor points that can form a liquid bridge connection, the liquid volume in the liquid bridge and the length of the contact interface are collected. To this end, the liquid bridge formation process in the array needs to be recorded frame by frame through high-resolution image capture, and then the liquid bridge contour is extracted by combining image processing tools (such as ImageJ or OpenCV). The liquid volume can be calculated by voxel integration, that is, by calculating the volume projection of the liquid bridge in the image and converting the volume based on the thickness of the array material. The contact interface length can be obtained by edge detection algorithms (such as Canny edge detection) to obtain the contour line and then estimate the length using the Euclidean distance formula. In a 10×10 nanopore array, for nanopores numbered (2, 3) and (2, 4), the liquid bridge interface contour is an approximate elliptical structure, and the equivalent diameter of the liquid bridge after image recognition is 0.3 μm, and the interface length is approximately π×0.15×0.1≈0.047 μm. After performing the same processing on the anchor point pairs, a set of two-variable mapping data sets containing liquid volume (in μm³) and the corresponding contact interface length (in μm) is formed. Then, a liquid volume-interface length table is constructed, and the liquid adhesion migration limit of each pair of anchor points under different humidity conditions is recorded. This limit can be set using the boundary value experiment method, that is, gradually increasing the interface length or liquid volume until the liquid bridge breaks or migration stops. For example, a migration critical test is performed on the (2, 3)-(2, 4) anchor point pair. Under the conditions of 25°C and 50% RH, stable adhesion is still achieved when the liquid bridge volume is 0.08 μm³ and the interface length is 0.045 μm. However, migration failure occurs when the liquid bridge volume reaches 0.06 μm³ and 0.041 μm. Therefore, the migration limits are set to 0.07 μm³ and 0.043 μm and included in the numerical table to generate the liquid bridge characteristic numerical table.

[0021] S102: Recalling the liquid adhesion migration limit and contact interface length data between anchor point pairs in the liquid bridge characteristic value table, calculating the corresponding projection trajectory length based on the two-dimensional plane projection coordinates between adjacent nanopore liquid bridge anchor points, sorting the anchor point pairs in ascending order, and removing the paths that do not meet the migration conditions in the sorted sequence based on the liquid adhesion migration limit between the anchor point pairs to obtain a sequence of migratable path lengths; The geometric position is determined according to the two-dimensional projection coordinates of adjacent anchor points in the nanopore array. The distance between anchor points (1, 1) and (1, 2) in the array is set to d = 0.5 μm, and the projection track length is μm. This process requires traversing adjacent anchor point pairs, calculating the trajectory lengths between anchor point pairs that can form liquid bridges in turn, and then sorting the anchor point pairs in ascending order of length, such as (1, 1)-(1, 2): 0.5 μm, (2, 2)-(2, 3): 0.6 μm, ... After sorting, liquid attachment migration condition screening is performed on each anchor point pair, and conditional judgment is performed on each pair of paths (i, j). If the migration limit value of the pair of anchor points in the characteristic value table (such as liquid bridge volume ≥ 0.07 μm³, interface length ≥ 0.045 μm) is not met, then this path is eliminated, and the remaining paths such as (1, 1)-(1, 2), (2, 2)-(2, 3), (2, 3)-(2, 4), etc. constitute the migration path length sequence { , , ,…}, Two conditional judgments need to be performed during the screening process: one is to calculate whether the volume and interface of the actual liquid bridge formation meet the limit standards, and the other is to calculate whether the path sorting process conforms to the ascending structure. For example, although the path (1, 3)-(1, 4) is short, the interface length is 0.035μm, which is lower than the limit of 0.043μm, so it is eliminated to obtain the transferable path length sequence.

[0022] S103: Based on the order and position distribution of the path starting points in the migratable path length sequence, the path starting points are sequentially controlled and the interleaving behavior of the path migration is tracked. The position segments where the overlapping sections appear in the path are identified and marked, thereby generating the initial dynamic configuration state of the liquid bridge network. The starting point control strategy is set in sequence. The edge nanopore is used as the starting point. The upper left corner of the array (1, 1) is set as the initial starting point. The shortest reachable path is selected for tracking. During the path migration process, the path needs to be analyzed section by section to see if there is any intersecting behavior, that is, whether the two paths have overlapping intersection areas in the two-dimensional projection. The overlap judgment is based on whether the coordinates of the projection trajectory are overlapped in any section. If the paths (1, 1)-(1, 2)-(1, 3) and (1, 1)-(2, 1)-(3, 1) intersect at point (1, 1), the intersection area is marked, and the coordinate matching is performed on the path. If there is an intersection in the path segment, the intersection area is marked. If it appears in two paths at the same time, it is judged as an intersection segment and the path segment number is recorded. For example, if (2, 2)-(2, 3) and (3, 2)-(2, 3) intersect at (2, 3), they are marked as segments. At the intersection point, after the complete path tracing is completed, the intersection and overlapping segment numbers and coordinates are recorded to generate the initial dynamic configuration state of the liquid bridge network.

[0023] The specific steps for obtaining the nanopore interface morphology response shift map are as follows: S201: extracting a two-dimensional projection trajectory from the edge of the liquid bridge to the center of the corresponding nano-orifice according to the initial dynamic configuration state of the liquid bridge network, recording the path coordinate sequence of the projection trajectory in the initial coordinate system, and uniformly mapping the trajectory to a unified coordinate system to generate a projection path sequence from the liquid bridge to the orifice; The formation and movement of liquid bridges in a porous array are visually tracked. A standard coordinate framework for the array is constructed in the initial state, with a two-dimensional coordinate system set with the array's central orifice as the origin. Each orifice position is assigned a unique coordinate number. High-speed microscopy is used to record the dynamic behavior of the liquid bridge throughout its entire process. Image recognition tools are used to extract the liquid bridge paths from the video frame sequence. The liquid bridge motion trajectory formed between the anchor point and the orifice is extracted as a trajectory data stream consisting of several path coordinate points. Each trajectory contains a two-dimensional path from the central anchor point to a specific orifice. The path data is converted into trajectory information within a unified coordinate system. By comparing the relative offsets of the path points with the original anchor point positions in each frame, the paths are mapped and unified in the standard coordinate system. This operation facilitates the comparability of liquid bridge trajectories formed at different time points and in different experimental batches, laying the foundation for subsequent processing. In observing nanoporous materials with regular array structures, the two-dimensional trajectories formed by each liquid bridge path during its movement can be classified, standardized, and archived, generating a sequence of projected paths from the liquid bridge to the orifice.

[0024] S202: Based on the liquid bridge to orifice projection path sequence, identify the coordinate value of the trajectory endpoint on the nano-orifice surface, compare the endpoint coordinates of multiple time sections between the same anchor points, calculate the trajectory endpoint displacement equilibrium characteristic value, extract the peak change and record the coordinate index position, and generate the trajectory endpoint displacement index distribution value; The displacement equilibrium characteristic value of the trajectory end point is calculated using the formula: ; in, Represents anchor pair 、 The displacement equilibrium eigenvalue of the trajectory endpoint between Represents an anchor point In the time section The horizontal end point coordinate value of Represents an anchor point In the time section The horizontal end point coordinate value of Represents an anchor point In the time section The vertical end point coordinate value, Represents an anchor point In the time section The vertical end point coordinate value, Represents an anchor point In the time section The local displacement increment, Represents an anchor point 、 The average value of the displacement increment in the local direction in all time sections, represents the number of time sections, The sequence index value representing the time section, 、 Represents the selected anchor point number; Formula is used to calculate anchor pair 、 The displacement equilibrium eigenvalues ​​of the end point coordinate sequence at different time sections are given by: It represents the average square difference of the anchor point pair in the horizontal coordinate direction, reflecting the consistency of the horizontal displacement of the two anchor points; the second term It represents the average value of the product of the displacement difference of the anchor point pair in the vertical coordinate direction and the incremental deviation of the local displacement, reflecting the degree of coordinated change of the two anchor points in the vertical displacement. By summing the absolute values ​​of these two items and then taking the square root, we can obtain an index that comprehensively reflects the consistency of the spatial displacement of the two anchor points. The trajectory endpoint displacement equilibrium eigenvalue is a numerical indicator that measures the consistency of the endpoint coordinate changes in space for the same anchor point pair at multiple time sections. It reflects the displacement differences in the horizontal and vertical directions of the anchor point pair and the degree of their fluctuation coordination. To further illustrate how to obtain each parameter, the following lists the specific calculation process and sample data: 、 、 、 :Through high-precision image recognition equipment, at each time section Extract anchor points and The horizontal and vertical end point coordinates, in the time section When, anchor point The horizontal coordinate of , the vertical coordinate is ; Anchor The horizontal coordinate of , the vertical coordinate is ; : Indicates the anchor point In the time section The local displacement increment is obtained by calculating the vertical coordinate difference between adjacent time sections. and Up, anchor point The vertical coordinates are and ,but ; : Indicates the anchor point and The average value of the local displacement increment on the time section is set at On a time section, the anchor point The local displacement increments are 、 、 ,but ; : represents the total number of time sections, which is determined by experimental design. If 3 time section data are collected in the experiment, then ; Substitute the above sample data into the formula for calculation: ; ; Comprehensive calculation: ; This result shows that the anchor and The displacement equilibrium eigenvalue in the selected time section is , indicating that the consistency of the two anchor points in spatial displacement is high; The benefit of this formula is that, by introducing the product of the vertical displacement difference of the anchor point pair and the local displacement increment deviation, it can more comprehensively reflect the degree of coordinated change of the two anchor points in spatial displacement and improve the accuracy of the displacement equilibrium eigenvalue. To further illustrate the parameter acquisition process, the experimental data are listed below: Table 1: Experimental data table

[0025] As shown in Table 1, the anchor points are listed in three time sections. and The horizontal and vertical end coordinates, as well as the anchor point The local direction displacement increment of the trajectory can be calculated according to the above formula through the data, and the displacement equilibrium characteristic value of the trajectory end point can be calculated according to the above formula.

[0026] S203: Mapping the spatial distribution area of ​​the original liquid around the nanopore according to the corresponding index position in the displacement index distribution value at the end point of the trajectory, uniformly translating the liquid distribution area in a two-dimensional coordinate system, determining the translation direction according to the direction indicated by the displacement index, and determining the translation amplitude using the displacement increment to generate a liquid distribution translation coordinate set; The position offset of the end point of the liquid bridge at different time sections is converted into the transformation information of the original liquid spatial distribution. The translation direction and amplitude contained in the displacement index are used as the spatial transformation reference. The liquid distribution pattern around the original nanopore is migrated as a whole to construct a new coordinate set of the translated liquid area. During the operation, each end point offset in the displacement index needs to be read one by one, and the represented spatial direction and movement amplitude are applied to the corresponding liquid distribution pattern. The liquid area around each orifice is spatially displaced according to the transformation to obtain a new liquid area layout, which is set at the orifice in the 5th row and 4th column of the array. The liquid area is originally concentrated in the center, but after the trajectory end point displacement conversion, it is offset to the edge of the orifice. In this case, the translation transformation needs to be accurately recorded and summarized to form a coordinate set, which will serve as the basis for the next step of spatial structure reconstruction. After integrating the translation coordinate information of the anchor point, the liquid distribution translation coordinate set is generated.

[0027] S204: calling the liquid distribution translation coordinate set, reorganizing the liquid distribution area according to the corresponding nanopore number, and reconstructing the interface contour graph under the standard pore contour boundary, plotting the boundary contour area corresponding to the liquid distribution difference before and after the displacement, and obtaining the nanopore interface morphology response offset map; Each nanopore is numbered and calibrated accordingly, the liquid distribution areas before and after migration are reorganized, and the interface contour graphics are reconstructed against the standard pore contour. The liquid distribution in the original state is spatially matched with its corresponding translation state graphics, and the boundary contours in the two states are drawn. The difference between the envelope of the boundary line and the outer contour is calculated to generate the interface morphology boundary change area caused by the change in liquid distribution. In the graphic reconstruction process, the liquid boundary contours before and after migration need to be marked separately, and the difference area is highlighted in a filled manner. The liquid movement process experienced by each nanopore is fully characterized by parameters such as the degree of boundary layer movement, direction and degree of contour overlap. In actual observations, it can be found that some pores form asymmetric migration graphics due to uneven surface structure. Such differences can be identified and annotated through the offset map, and a complete pore response graphic structure set is constructed to obtain the nanopore interface morphology response offset map.

[0028] The specific steps for obtaining the nanopore boundary stability response signal are as follows: S301: Invoking the nanopore interface morphology response offset map, extracting the continuous numerical value of the pore edge curvature along the boundary coordinate direction, and collecting the length change value of the corresponding liquid bridge path and the molecular contact frequency sequence of the liquid bridge end contact area, recording the numerical fluctuation interval on the time axis, and generating a curvature path frequency joint series; The curvature value of the pore edge is extracted in a unified coordinate system. Samples are taken point by point along the boundary coordinate direction. The boundary curvature value of each coordinate point is continuously recorded using contour tracing technology, and the values ​​are combined to form a complete edge curvature sequence. At the same time, the length changes of each liquid bridge path at different time periods are extracted from the normalized liquid bridge trajectory. The change in coordinates between the starting and end points of the path is summarized to obtain the change in liquid bridge path length, and the fluctuation of the value over time is recorded. The contact behavior of molecules in the contact area between the end of the liquid bridge and the nanopore surface is statistically analyzed. The contact frequency change sequence per unit time is recorded using an orbit tracing algorithm or molecular dynamics simulation method, and a contact frequency data set that varies with time is gradually accumulated. The data is synchronized with the curvature sequence and path change value, and classified according to the time axis. The numerical fluctuation range of each parameter at different time nodes is identified. The edge curvature, path length change and molecular contact frequency data are combined to generate a joint series of curvature path frequencies.

[0029] S302: Based on the hole edge curvature sequence in the curvature path frequency joint sequence, the curvature value change relationship between consecutive coordinate points is detected in the order of boundary coordinates. For the positions where the adjacent positive and negative incremental directions switch, the periodic direction change points are extracted, and the index segments where the curvature trend reverses are marked to obtain the edge curvature direction change index set; The curvature value change trend between consecutive coordinate points along the boundary contour coordinate sequence is checked one by one, and the positive and negative incremental changes in the curvature values ​​of two adjacent points are recorded. If the curvature value of a certain segment shows a continuous increasing trend and then begins to turn into a decreasing trend at a certain point, it indicates that this point is a trend reversal node, that is, a periodic turning point. In this operation, the boundary coordinate points are scanned in order of number, and the difference judgment is performed on each pair of consecutive curvature values ​​and their change direction is recorded. When a point has a direction opposite to the previous and next changes, it is marked as a trend reversal position and its index number is marked in the boundary sequence. At the same time, the continuous reversal points are connected into index segments. The 35th point in a certain orifice boundary is set as the critical point from increase to decrease, and the 46th point is set from decrease to increase, indicating that 35 to 46 is a complete curvature turning segment. The index segment represents the local contour area where the curvature change trend is reversed. The area often corresponds to the concentrated area of ​​microscopic deformation of the orifice and plays a key role in the subsequent path and molecular behavior analysis. The edge curvature turning index set is obtained.

[0030] S303: Calling the liquid bridge path length change amplitude in the edge curvature direction change index set, rearranging the coordinates of the path change amplitude and normalizing the direction according to the index segment, selecting the path segment with the same change trend as the direction change segment as the convergence segment, and calculating the difference range of the molecular contact frequency fluctuation in the corresponding segment to obtain the nanopore boundary stability response signal; The liquid bridge path length change amplitude data corresponding to each direction-changing segment is called, and the spatial coordinates of the path change values ​​are rearranged according to the segment range of the curvature change. That is, the path start and end coordinates are adjusted to the coordinate position consistent with the boundary index, unifying the analysis starting point of the path and curvature change. At the same time, the direction of the path length change is normalized, and the data in the increasing or decreasing direction are uniformly labeled to facilitate subsequent comparison and screening. The consistency of the path change direction and the trend of the direction-changing segment is determined. If the segment where the curvature changes from increasing to decreasing is accompanied by the path changing from shortening to growing, it can be considered that the path segment is consistent with the curvature change trend and will be screened as a convergent segment. The difference range of molecular contact frequency is extracted in the convergent segment, that is, the extreme value difference of the frequency sequence is extracted within several frames before and after the change, and the difference is classified and counted to obtain the molecular contact strength response range of each convergent segment, which is used for interface state modeling analysis and abnormal warning identification, and to obtain the nanopore boundary stability response signal.

[0031] The specific steps for obtaining the boundary area limited correction field distribution map are as follows: S401: extracting the liquid bridge anchor point number marked as the unstable section based on the nano-pore boundary stability response signal, recording the liquid coverage area under the corresponding anchor point and the length of the tension conduction path to the pore boundary, and generating a liquid bridge anchor point coverage tension value group; Clearly identify the segments marked as unstable. These segments correspond to areas where the liquid bridge exhibits structural anomalies or dramatic fluctuations in contact frequency during the formation of specific anchor points. The signal sequence needs to be parsed segment by segment, and the anchor point numbers marked as unstable are extracted. An index list of corresponding anchor points is then established. Based on this, the area of ​​the liquid region covered by each unstable anchor point is calculated. This operation relies on the closed figure formed by the contour boundary of the liquid region. The liquid coverage area value is obtained through contour tracing and pixel area conversion. The length of the tension conduction path formed from the anchor point to the orifice boundary is recorded. Specifically, the actual path length is calculated along the shortest continuous channel between the liquid bridge and the orifice boundary. Path length extraction requires reference to the projection change trajectory of the liquid bridge in the three-dimensional structure and is estimated based on the path projection in two-dimensional coordinates. At each unstable anchor point, the corresponding liquid coverage area and tension path length are simultaneously obtained. These two data are classified and stored with the anchor point number to generate a liquid bridge anchor point coverage tension value set.

[0032] S402: Calling the liquid bridge anchor point coverage tension value group, numerically amplifying the coverage area according to the set expansion ratio, calculating the amplified area characteristic value, and converging the tension path length to the set tension conduction constraint range. The tension path length is uniformly mapped to the two-dimensional nanopore boundary coordinate system, and the tension distribution area contour of the anchor point position in the boundary area is re-plotted to obtain the boundary area confined correction field distribution map; The characteristic value of the amplified area is calculated using the formula: ; in, represents the characteristic value of the amplified area, represents the initial area covered by the foundation liquid bridge anchor point, Representative Anchor point tension transmission gain factor, Representative Item anchor point tension value, Representative Item Anchor point path length compression damping coefficient, Representative Item anchor path length, Represents the equivalent radius length formed by mapping the tension path in the two-dimensional plane, represents the average conducted disturbance displacement correction value, Represents the number of anchor points; In this formula, the base area It represents the projected area of ​​the liquid bridge anchor point before expansion, which is obtained from the measured data of the liquid bridge morphology and calculated by optical imaging combined with image segmentation to extract the area contour; the tension value The micro tension sensor is embedded in the anchor structure to measure the tension fluctuation average value under different working conditions and process it through amplitude equalization; the path length The linear distance between points is obtained along the tension direction by using a laser displacement meter; the coefficient 、 The tension gain and damping compression weight coefficients represent the positive push ratio of tension to area expansion and the negative compression ratio of path length to area increase, respectively. The two settings are based on the process test calibration values. The settings are as follows: the experimental samples are divided into three groups, and the area expansion rate under independent tension changes is tested (to obtain the tension effect coefficient ), and the area shrinkage under the change of path length (obtaining the compression effect coefficient ), using the linear fitting coefficient as the source of the coefficient value; The correction value for the displacement of the tension transmission disturbance is obtained by quantifying the position change of the anchor point after being disturbed by the surrounding environment in multiple experiments. It is obtained by taking the square of the difference between the maximum displacement under disturbance and its static mean value. is the equivalent radius of the tension path, which is calculated by the radius of the circumscribed circle formed by the spatial mapping of the end of the anchor tension path; The amplified area characteristic value refers to the effective coverage area under the action of anchor point tension expansion. It is the area correction result obtained after considering factors such as tension strength, path complexity and disturbance correction. It reflects the maximum effective influence range that the liquid bridge anchor point can actually achieve in space due to physical driving. Parameter acquisition and actual calculation examples: Taking the experimental sample A as an example, the basic area is extracted by image recognition in the initial state. ; There are 4 anchor points in total, and the number of effective conduction anchor points is , one anchor point was removed due to insufficient tension; The tension values ​​measured by the tension sensor are ; The path lengths were measured to be ; The weight coefficients determined by the experiment are: (Fitted by tension effect rate); (normalized by path compression sensitivity experiments); Disturbance correction value The maximum offset of each anchor point is obtained by monitoring with a laser interferometer and then calculated. ; Path projection circumcircle radius ; Substitute each item into the formula and calculate as follows: ; ; ; The sum is: ; Substituting into the formula: ; The results show that the area eigenvalue increased to 73.7 μm², and after the tension-path joint parameter correction, the area expansion has reached 5.9 times the base area. This area eigenvalue will be used in the subsequent steps to redraw the boundary of the tension distribution area. The benefit of the formula is that it constructs the area contribution term by introducing the absolute difference between the tension gain coefficient and the path compression coefficient, and performs normalized square root correction in combination with the tension path space projection and the disturbance correction term. This allows the area gain processing to not only reflect the actual tension distribution difference, but also effectively integrate the path complexity and the offset caused by the disturbance, providing more reasonable data support for subsequent boundary redrawing.

[0033] Table 2: Area increase parameter monitoring table

[0034] Table 2 lists the parameter data of each anchor point before area expansion and the corresponding corrected calculated values, which serve as the components of the summation term in the formula.

[0035] The specific steps for obtaining the trigger mark of the nanopore on-off state transition are as follows: S501: Calling the boundary area limited domain correction field distribution map, monitoring the time series temperature change curve of the area covering the thermosensitive composite material, calculating the temperature change rate of the area per unit time, synchronously extracting the window interval value, and arranging the temperature rate and response window value pairs according to the coordinate number to generate a thermosensitive response rate window pairing table; It is necessary to determine the area of ​​thermosensitive composite material covered in the figure and perform time series monitoring of its continuous temperature changes. Non-contact high-resolution infrared thermal imaging is used to record the surface temperature field of the material in frames, and the temperature values ​​at continuous time points are collected at a set frame rate to form a complete temperature change curve sequence. For each calibrated two-dimensional coordinate number position, the temperature change rate per unit time is calculated. The method is to take the temperature difference between continuous time points and divide it by the time interval to obtain the temperature rate value of each area. The temperature rate sequence is then segmented and extracted according to the preset time window parameters. The extracted values ​​are used as representative responses within the window interval. A combined data set is established with the coordinate number as the index, the temperature rate as the first value, and the temperature change within the window interval as the second value to generate a thermosensitive response rate window pairing table.

[0036] S502: Based on the thermal response rate window pairing table, the film thickness data at the nanopore edge position and the variable area ratio data in the pore mouth area are extracted to identify the corresponding structure between the pore edge structure size and the response ability. The rate response lag time is compared with the duration of thermal boundary deformation. The periodic continuous change segments are identified in the time series order to obtain the response period continuous segment index sequence; The membrane thickness information at the nanopore edge position corresponding to the pore area is extracted, and the membrane structural dimension data at each pore edge position is obtained using nano-tomography or atomic force microscopy measurement equipment. In combination with the ratio of the variable area to the overall area in the pore area under different thermal response states, a corresponding structural table of pore edge size and response ability is constructed. The structure and rate response data are correlated and compared, focusing on the parts with thicker membranes or higher variable area ratios to see if there is a significant shift in their thermal response rate and hysteresis characteristics. The rate response lag time, that is, the time interval between the thermal stimulus input and the material's obvious temperature response, is retrieved from the thermal response curve. The duration of each thermal interface deformation is also recorded. Both are used as key response parameters and arranged in time series order. The lag time is compared to see if it is shorter or longer than the response duration. The segments with repeated rhythmic changes are retrieved to form continuous segments with periodic thermal response characteristics. Each period segment is indexed and a response period continuous segment index sequence is output.

[0037] S503: Calling the response cycle continuous segment index sequence, identifying whether the rate response hysteresis in the cycle segment has completed the return judgment within the material response window interval, selecting the cycle segment that satisfies the synchronous closure of the deformation duration and the response hysteresis time, marking the corresponding cycle closure information by number, and generating the nanopore on-off state transition trigger mark; It is necessary to judge whether the cycle meets the conditions for the rate response lag to complete reset within the cycle segment, that is, whether the process of the material from being heated to generating a response and then returning to the original state is closed within the cycle within the complete cycle. During the operation, it is necessary to compare the thermal response window interval to judge whether the lag time is always contained in the material response window. If the lag behavior within the cycle exceeds the response interval, it is judged as an unclosed cycle segment, otherwise it is a synchronously closed cycle segment. In the cycle segment, the cycle segment that meets the synchronous closure of the lag reset and the deformation duration is screened out, and its cycle number is specially marked. A closed mark is added after the cycle segment to form a complete nanopore response cycle labeling system, which provides a control signal for each cycle that meets the material thermal response stable closure conditions, provides a logical basis and time node basis for the pore mouth to realize the open and close switching state under micro-heat drive, and generates a nanopore on-off state transition trigger mark.

[0038] See also Figure 2 , a nanopore functional control system, the system comprising: The path recognition module obtains the coordinate values ​​of the nanopore liquid bridge anchor points in the porous array, the liquid volume value, the contact interface length value, and the attachment migration limit value. It calculates the sorting sequence based on the projected trajectory length between the liquid bridge anchor points. Based on the sorting sequence, it controls the arrangement order of the starting points of the liquid migration path to generate the initial dynamic configuration state of the liquid bridge network. The interface response module extracts the two-dimensional projection trajectory data from the liquid bridge edge to the nanopore according to the initial dynamic configuration state of the liquid bridge network, determines the change direction between the trajectory endpoint and the displacement vector of the pore area, and generates a nanopore interface morphology response offset map; The structural stability module calculates the position of the continuity inflection point within the curvature change range based on the nanopore interface morphology response deviation map, selects the boundary dividing point between the structural stability range and the instability range, and obtains the nanopore boundary stability response signal; The confinement correction module calls the nanopore boundary stability response signal, extracts the liquid coverage area value and tension path length value of the liquid bridge anchor point in the instability zone, adjusts the coverage area range according to a fixed expansion ratio, and generates a confinement correction field distribution map of the boundary area; The thermal trigger module detects the temperature change rate value and the material response time window value of the thermosensitive composite material based on the limited correction field distribution map of the boundary area, calculates the period ratio between the temperature change lag time and the pore deformation time, and obtains the trigger mark of the nanopore on-off state transition.

[0039] The above description is merely a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any modifications or substitutions that can be easily conceived by a person skilled in the art within the technical scope disclosed in the present invention should be included in the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be based on the scope of protection of the claims.

Claims

1. A method for controlling nanopore functionality, characterized in that: The following steps are involved: S1: Based on the initial distribution state of nanopore liquid bridge anchor points in the porous array, the liquid volume, contact interface length and adhesion migration limit are collected to screen the staggered behavior segments in the migration path and generate the initial dynamic configuration state of the liquid bridge network; S2: extracting the projection trajectory from the edge of the liquid bridge to the nanopore according to the initial dynamic configuration state of the liquid bridge network, identifying the displacement change of the end point of the trajectory, performing positional translation on the liquid distribution, and generating a nanopore interface morphology response shift map; S3: calling the nanopore interface morphology response deviation map, extracting the pore edge curvature sequence, the liquid bridge path length change amplitude and the molecular contact frequency fluctuation range, performing continuity detection on the curvature change trend and extracting the periodic turning points, performing projection sequence sorting on the path length change amplitude, identifying the segment that converges with the curvature change, and generating the nanopore boundary stability response signal; S4: According to the nanopore boundary stability response signal, the liquid coverage area and tension path value of the liquid bridge anchor point in the instability zone are extracted, the coverage area is amplified according to the set expansion ratio, and the conduction path is compressed to the constraint interval to generate a limited correction field distribution map of the boundary area.

2. The method for controlling nanopore functionality according to claim 1, wherein: The initial dynamic configuration state of the liquid bridge network includes the anchor point network topology, the path priority connection sequence, and the staggered path distribution pattern; the nanopore interface morphology response offset map includes the interface response offset contour, the displacement dynamic amplitude interval, and the trajectory terminal offset mapping result; the nanopore boundary stability response signal includes the stability period identifier, the path trend coincidence section indication, and the contact frequency fluctuation amplitude parameter; the boundary area confined correction field distribution map includes the tension compression mapping map, the coverage area expansion area, and the confined field strength change block.

3. The method for controlling nanopore functionality according to claim 1, wherein: The steps for obtaining the initial dynamic configuration state of the liquid bridge network are specifically as follows: S101: Based on the initial distribution state of the nanopore liquid bridge anchor points in the porous array, the liquid volume and the contact interface length in the nanopore liquid bridge are collected, and a corresponding mapping is established between the liquid volume formed between the anchor points and the contact interface length. The liquid attachment migration limit between each pair of anchor points is recorded to generate a liquid bridge characteristic value table; S102: Recalling the liquid adhesion migration limit and contact interface length data between anchor point pairs in the liquid bridge characteristic value table, calculating the corresponding projection trajectory length based on the two-dimensional plane projection coordinates between adjacent nanopore liquid bridge anchor points, sorting the anchor point pairs in ascending order, and removing paths that do not meet the migration conditions in the sorted sequence based on the liquid adhesion migration limit between the anchor point pairs to obtain a sequence of migratable path lengths; S103: Based on the order and position distribution of the path starting points in the migratable path length sequence, the path starting points are controlled in sequence and the interlaced behavior state of the path migration is tracked, the position segments where the cross-overlapping sections appear in the path are identified and marked, and the initial dynamic configuration state of the liquid bridge network is generated.

4. The method for controlling nanopore functionality according to claim 1, wherein: The steps for obtaining the nanopore interface morphology response shift map are specifically as follows: S201: extracting a two-dimensional projection trajectory from the edge of the liquid bridge to the center of the corresponding nano-orifice according to the initial dynamic configuration state of the liquid bridge network, recording a path coordinate sequence of the projection trajectory in the initial coordinate system, and uniformly mapping the trajectory to a unified coordinate system to generate a projection path sequence from the liquid bridge to the orifice; S202: Based on the liquid bridge to orifice projection path sequence, identifying the coordinate value of the trajectory endpoint on the nano-orifice surface, comparing the endpoint coordinates of multiple time sections between the same anchor points, calculating the trajectory endpoint displacement equilibrium characteristic value, extracting the peak value change and recording the coordinate index position thereof, and generating the trajectory endpoint displacement index distribution value; The displacement balance characteristic value of the trajectory endpoint is calculated using the formula: ； in, Represents anchor pair 、 The displacement equilibrium eigenvalue of the trajectory endpoint between Represents an anchor point In the time section The horizontal end point coordinate value of Represents an anchor point In the time section The horizontal end point coordinate value of Represents an anchor point In the time section The vertical end point coordinate value, Represents an anchor point In the time section The vertical end point coordinate value, Represents an anchor point In the time section The local displacement increment, Represents an anchor point 、 The average value of the displacement increment in the local direction in all time sections, represents the number of time sections, The sequence index value representing the time section, 、 Represents the selected anchor point number; S203: Mapping the spatial distribution area of ​​the original liquid around the nanopore according to the corresponding index position in the displacement index distribution value of the trajectory endpoint, uniformly translating the liquid distribution area in a two-dimensional coordinate system, determining the translation direction according to the direction indicated by the displacement index, and determining the translation amplitude using the displacement increment to generate a liquid distribution translation coordinate set; S204: calling the liquid distribution translation coordinate set, reorganizing the liquid distribution area according to the corresponding nanopore number, and reconstructing the interface contour graphic under the standard pore contour boundary, plotting the boundary contour area corresponding to the liquid distribution difference before and after the displacement, and obtaining the nanopore interface morphology response offset map.

5. The method for controlling nanopore functionality according to claim 4, wherein: The steps for obtaining the nanopore boundary stability response signal are specifically as follows: S301: calling the nanopore interface morphology response offset map, extracting the continuous numerical value of the pore edge curvature along the boundary coordinate direction, collecting the length change value of the corresponding liquid bridge path and the molecular contact frequency sequence of the liquid bridge end contact area, recording the numerical fluctuation interval on the time axis, and generating a curvature path frequency joint series; S302: Based on the hole edge curvature sequence in the curvature path frequency joint series, the curvature value change relationship between consecutive coordinate points is detected in the order of boundary coordinates. For the positions where the adjacent positive and negative incremental directions switch, the periodic direction change points are extracted, and the index sections where the curvature trend reverses are marked to obtain the edge curvature direction change index set; S303: Call the liquid bridge path length change amplitude in the edge curvature direction change index set, rearrange the coordinates and normalize the direction of the path change amplitude according to the index segment, select the path segment with the same change trend as the direction change segment as the convergence segment, count the difference range of the molecular contact frequency fluctuation in the corresponding segment, and obtain the nanopore boundary stability response signal.

6. The method for controlling nanopore functionality according to claim 5, wherein: The steps for obtaining the boundary area limited correction field distribution map are specifically as follows: S401: extracting the liquid bridge anchor point number marked as the unstable section according to the nano-orifice boundary stability response signal, recording the liquid coverage area under the corresponding anchor point and the length of the tension conduction path to the orifice boundary, and generating a liquid bridge anchor point coverage tension value group; S402: Call the liquid bridge anchor point coverage tension value group, perform numerical amplification processing on the coverage area according to the set expansion ratio, calculate the amplified area characteristic value, and at the same time converge the tension path length to the set tension conduction constraint range, uniformly map it to the two-dimensional nanopore boundary coordinate system, re-plot the tension distribution area contour of the anchor point position in the boundary area, and obtain the boundary area confined correction field distribution map.

7. The method for controlling nanopore functionality according to claim 6, wherein: The characteristic value of the amplified area is calculated using the formula: ； in, represents the characteristic value of the amplified area, represents the initial area covered by the foundation liquid bridge anchor point, Representative Anchor point tension transmission gain factor, Representative Item anchor point tension value, Representative Item Anchor point path length compression damping coefficient, Representative Item anchor path length, Represents the equivalent radius length formed by mapping the tension path in the two-dimensional plane, represents the average conducted disturbance displacement correction value, Represents the number of anchor points.

8. The method for controlling nanopore functionality according to claim 1, wherein: The method further comprises step S5: S5: calling the boundary area limited domain correction field distribution map, detecting the temperature change rate of the covering thermosensitive composite material and the material response window interval value, extracting the pore edge membrane thickness and the pore mouth variable area ratio data, comparing the execution rate response hysteresis and the periodic continuity between the thermal boundary deformation time, judging whether the continuous response behavior of the thermal trigger behavior has met the closing condition, and generating the nanopore on-off state transition trigger mark; The nanopore on-off state transition trigger mark includes a trigger threshold index, a response closure cycle, and an on-off behavior identification site.

9. The method for controlling nanopore functionality according to claim 8, wherein: The steps for obtaining the trigger mark of the nanopore on-off state transition are specifically as follows: S501: calling the boundary region limited domain correction field distribution map, monitoring the time series temperature change curve of the area covering the thermosensitive composite material, calculating the temperature change rate of the area per unit time, synchronously extracting the window interval value, arranging the temperature rate and response window value pairs according to the coordinate number, and generating a thermosensitive response rate window pairing table; S502: Based on the thermal response rate window pairing table, extract the film thickness data at the nanopore edge position and the variable area ratio data in the pore mouth area, identify the corresponding structure between the pore edge structure size and the response ability, compare the rate response lag time with the duration of thermal boundary deformation, identify the periodic continuous change section in time series order, and obtain the response period continuous section index sequence; S503: Call the response cycle continuous segment index sequence, identify whether the rate response lag in the cycle segment completes the return judgment within the material response window interval, screen the cycle segments that meet the synchronous closure of the deformation duration and the response lag time, mark the corresponding cycle closure information by number, and generate a nanopore on-off state transition trigger mark.

10. A nanopore functional control system, characterized in that: The system is used to implement the nanopore functionality control method according to any one of claims 1 to 9, and the system comprises: The path recognition module obtains the coordinate values ​​of the nanopore liquid bridge anchor points in the porous array, the liquid volume value, the contact interface length value, and the attachment migration limit value. It calculates the sorting sequence based on the projected trajectory length between the liquid bridge anchor points. Based on the sorting sequence, it controls the arrangement order of the starting points of the liquid migration path to generate the initial dynamic configuration state of the liquid bridge network. The interface response module extracts the two-dimensional projection trajectory data from the edge of the liquid bridge to the nanopore according to the initial dynamic configuration state of the liquid bridge network, determines the change direction between the trajectory end point and the displacement vector of the orifice area, and generates a nanopore interface morphology response offset map; The structural stability module calculates the position of the continuity inflection point within the curvature change range according to the nanopore interface morphology response deviation map, selects the boundary dividing point between the structural stability range and the instability range, and obtains the nanopore boundary stability response signal; The confinement correction module calls the nanopore boundary stability response signal, extracts the liquid coverage area value and the tension path length value of the liquid bridge anchor point in the instability zone, adjusts the coverage area range according to a fixed expansion ratio, and generates a confinement correction field distribution map of the boundary area; The thermal trigger module detects the temperature change rate value and the material response time window value of the thermosensitive composite material based on the limited correction field distribution map of the boundary area, calculates the period ratio between the temperature change lag time and the pore deformation time, and obtains the trigger mark of the nanopore on-off state transition.

Images