A method and system for controlling the functionality of nanopores

Abstract

This invention relates to the field of nanopore detection technology, specifically to a method and system for functional control of nanopores, comprising the following steps: extracting path information based on the distribution of liquid bridge anchor points, constructing a dynamic configuration state, generating an interface response spectrum, identifying stability signals, correcting the confined field distribution, monitoring temperature-sensitive response behavior, and forming a trigger marker for on / off state transitions. In this invention, by introducing sequential control and staggered behavior recognition of liquid migration paths, the dynamic configuration of the liquid bridge distribution is optimized. Combined with interface response mapping based on displacement peak indices, the accuracy and efficiency of liquid position control are effectively improved. Through the extraction of convergence between path and boundary responses and periodic direction analysis, dynamic prediction of the stable state of the nanopore opening boundary is achieved. Through the expansion of liquid coverage area and compression of tension paths, the stability and controllability in scenarios of molecular on / off control and dynamic structural recognition are enhanced.

Description

Technical Field

[0001] This invention relates to the field of nanopore detection technology, and in particular to a method and system for controlling the functionality of nanopores. Background Technology

[0002] Nanopore detection technology involves the precise identification and analysis of substances using nanoscale pores. Core aspects include the fabrication and functionalization of nanopores, as well as molecular-level detection and assays. Nanopore detection technology is widely used in 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 nanomedicines. 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 gene analysis. Nanopore functional control methods refer to optimizing and controlling the performance of nanopores by adjusting their geometry, material composition, or electrical environment. This involves adjusting the size, shape, surface modification, and electric field of the nanopores 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 means of chemical modification, addition of surfactants, and electric field modulation, so as to ensure accuracy and efficiency in applications.

[0003] While existing technologies can screen and detect molecules by controlling the size of nanopores and adjusting the electric field environment, they lack a mechanism for identifying interlacing behavior during dynamic liquid migration and interface morphology responses. They cannot achieve real-time control of the interactions between complex liquid bridge structures, lack the ability to judge periodic trends in the continuous analysis of interface responses, and are prone to unstable changes in boundary states. When it comes to the response of temperature-sensitive materials, they have failed to establish a collaborative judgment standard between structural thickness and response ratio, resulting in increased bias in the identification of on / off states. In practical applications, they are prone to problems such as judgment lag or response interruption in temperature-controlled and high-sensitivity detection scenarios, affecting the repeatability and accuracy of the overall performance of nanopores. Summary of the Invention

[0004] To address the shortcomings of existing technologies, such as the lack of identification mechanisms for interlacing behaviors during dynamic liquid migration and interface morphology responses, the inability to achieve real-time control of interactions between complex liquid bridge structures, the lack of periodic trend judgment capabilities in the continuity analysis of interface responses leading to unstable changes in boundary states, and the failure to establish a collaborative judgment standard between structural thickness and response ratio when involving temperature-sensitive material responses, resulting in increased bias in on / off state identification, and the potential for judgment lag or response interruption in practical applications under temperature-controlled and high-sensitivity detection scenarios, thus affecting the repeatability and accuracy of the overall performance of nanopores, this invention provides a method and system for functional control of nanopores. The technical solution is as follows:

[0005] On the one hand, a method for controlling the functionality of nanopores is provided, the method comprising:

[0006] S1: Based on the initial distribution state of nanoporous liquid bridge anchor points in the porous array, the liquid volume, contact interface length and adhesion migration limit are collected, the intersecting behavior segments in the migration path are screened, and the initial dynamic configuration state of the liquid bridge network is generated.

[0007] S2: Based on the initial dynamic configuration state of the liquid bridge network, extract the projection trajectory from the edge of the liquid bridge to the nanopore, identify the displacement change at the endpoint of the trajectory, perform position translation on the liquid distribution, and generate a nanopore interface morphology response offset map.

[0008] S3: Call the nanopore interface morphology response offset map, extract the pore edge curvature sequence, the liquid bridge path length change amplitude and the molecular contact frequency fluctuation range, perform continuous detection on the curvature change trend and extract periodic change points, perform projection sequence organization on the path length change amplitude, identify the segment that is similar to the curvature change, and generate the nanopore boundary stability response signal.

[0009] S4: Based on the stability response signal of the nanopore boundary, extract the liquid coverage area and tension path value of the liquid bridge anchor point in the unstable region, increase the coverage area according to the set expansion ratio, and compress the conduction path to the constraint interval to generate a boundary region confined correction field distribution map.

[0010] As a further aspect 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 profile, the displacement dynamic amplitude range, and the trajectory terminal offset mapping result; the nanopore boundary stability response signal includes the stability period identifier, the path trend matching section indicator, and the contact frequency fluctuation amplitude parameter; and the boundary region confined correction field distribution map includes the tension compression mapping map, the coverage area expansion region, and the confined field strength change block.

[0011] As a further aspect of the present invention, the step of obtaining the initial dynamic configuration state of the liquid bridge network specifically includes:

[0012] S101: Based on the initial distribution state of the anchor points of the nanoporous liquid bridge in the porous array, the liquid volume and contact interface length in the nanoporous 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 adhesion migration limit between each pair of anchor points is recorded, and a liquid bridge characteristic value table is generated.

[0013] S102: Call the liquid adhesion migration limit and contact interface length data between anchor point pairs in the liquid bridge feature value table, calculate the corresponding projection trajectory length based on the two-dimensional plane projection coordinates between adjacent nanopore liquid bridge anchor points, sort the anchor point pairs in ascending order, and eliminate paths in the sorted sequence that do not meet the migration conditions by combining the liquid adhesion migration limit between anchor point pairs, and obtain the migration path length sequence.

[0014] S103: Based on the path starting point order and position distribution in the migrated path length sequence, control the path starting point in sequence and track the staggered behavior state of path migration, identify the position segments of the path where there are intersecting and overlapping sections and filter and mark them to generate the initial dynamic configuration state of the liquid bridge network.

[0015] As a further aspect of the present invention, the step of obtaining the nanopore interface morphology response shift spectrum specifically includes:

[0016] S201: Based on the initial dynamic configuration state of the liquid bridge network, extract the two-dimensional projection trajectory from the edge of the liquid bridge to the center of the corresponding nanopore, record the path coordinate sequence of the projection trajectory in the initial coordinate system, and map the trajectory to a unified coordinate system to generate the projection path sequence from the liquid bridge to the pore.

[0017] S202: Based on the projection path sequence from the liquid bridge to the orifice, identify the coordinate value of the trajectory endpoint on the surface of the nano-orifice, compare the endpoint coordinates of multiple time sections between the same anchor point, calculate the trajectory endpoint displacement equilibrium characteristic value, extract the peak change amount and record the coordinate index position, and generate the trajectory endpoint displacement index distribution value.

[0018] The trajectory endpoint displacement equalization characteristic value is obtained using the formula:

[0019] ;

[0020] in, Represents anchor point pair , The trajectory endpoint displacement equilibrium characteristic value between them Represents anchor point In time section The horizontal endpoint coordinates, Represents anchor point In time section The horizontal endpoint coordinates, Represents anchor point In time section The vertical endpoint coordinates, Represents anchor point In time section The vertical endpoint coordinates, Represents anchor point In time section The local directional displacement increment, Represents anchor point , The average value of the local displacement increment over all time sections. Represents the number of time segments. The sequence index value representing the time section. , Represents the selected anchor point number;

[0021] S203: Based on the corresponding index position in the displacement index distribution value of the trajectory endpoint, map the spatial distribution area of ​​the original liquid around the nanopore, perform a uniform translation of the liquid distribution area in the two-dimensional coordinate system, determine the translation direction according to the direction indicated by the displacement index, determine the translation amplitude using the displacement increment, and generate a liquid distribution translation coordinate set.

[0022] S204: Call the liquid distribution translation coordinate set, reorganize the liquid distribution area according to the corresponding nanopore number, reconstruct the interface contour graphic under the standard pore contour boundary, plot the boundary contour area corresponding to the difference in liquid distribution before and after displacement, and obtain the nanopore interface morphology response offset map.

[0023] As a further aspect of the present invention, the step of obtaining the nanopore boundary stability response signal specifically includes:

[0024] S301: Call the nanopore interface morphology response offset map, extract the continuous values ​​of the pore edge curvature along the boundary coordinate direction, and collect the length change value of the corresponding liquid bridge path and the molecular contact frequency sequence of the contact area at the end of the liquid bridge. Record the numerical fluctuation range on the time axis to generate the curvature path frequency joint sequence.

[0025] S302: Based on the hole edge curvature sequence in the curvature path frequency joint sequence, detect the curvature value change relationship between continuous coordinate points in the boundary coordinate order, extract the periodic change point at the position where the adjacent positive and negative incremental directions change, mark the index segment where the curvature trend reverses, and obtain the edge curvature change index set.

[0026] S303: Call the edge curvature direction change index to concentrate the liquid bridge path length change range, rearrange the coordinates and normalize the direction of the path change range according to the index segment, select the path segment with the same change trend as the convergence segment, and count the difference range of molecular contact frequency fluctuation in the corresponding segment to obtain the nanopore boundary stability response signal.

[0027] As a further aspect of the present invention, the step of obtaining the boundary region confined correction field distribution map specifically includes:

[0028] S401: Based on the stability response signal of the nanopore boundary, extract the liquid bridge anchor point number marked as the unstable section, record the liquid coverage area under the corresponding anchor point and the tension transmission path length to the pore boundary, and generate a liquid bridge anchor point coverage tension value group.

[0029] S402: Call the tension value group of the liquid bridge anchor point coverage, perform numerical amplification processing on the coverage area according to the set expansion ratio, calculate the characteristic value of the amplified area, 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, redraw the tension distribution area outline of the anchor point position in the boundary area, and obtain the boundary area confinement correction field distribution map.

[0030] As a further aspect of the present invention, the amplification area characteristic value is expressed by the formula:

[0031] ;

[0032] in, The characteristic value representing the increase in area, This represents the initial area covered by the basic liquid bridge anchor point. Representing the Anchor point tension conduction gain factor Representing the Anchor tension value, Representing the Anchor point path length compression damping coefficient Representing the Item anchor path length, This represents the equivalent radius length formed by the mapping of the tension path in a two-dimensional plane. This represents the average conducted interference displacement correction value. This represents the number of anchor points.

[0033] As a further aspect of the present invention, the method further includes step S5:

[0034] S5: Call the boundary region confined correction field distribution map, detect the temperature change rate and material response window interval value of the covered thermosensitive composite material, extract the thickness of the pore edge film and the proportion of the variable region of the pore opening, compare the periodic continuity between the execution rate response lag and the thermal boundary deformation time, determine whether the continuous response behavior of the thermal triggering behavior has completed the closing condition, and generate a nanopore on / off state transition trigger mark.

[0035] The nanopore on / off state transition trigger marker includes a trigger threshold index, a response closure period, and an on / off behavior identification site.

[0036] As a further aspect of the present invention, the step of obtaining the nanopore on / off state transition trigger marker specifically includes:

[0037] S501: Call the boundary region confined correction field distribution map, monitor the time series temperature change curve of the area covered by the thermosensitive composite material, calculate the temperature change rate per unit time of the area, extract the window interval value synchronously, organize the temperature rate and response window value pair according to the coordinate number, and generate a thermosensitive response rate window pairing table.

[0038] S502: Based on the thermosensitive response rate window pairing table, extract the film thickness data at the nanopore edge position and the variable region ratio data in the pore area, identify the corresponding structure between the pore edge structure size and the response capability, compare the rate response lag time and the duration of thermal boundary deformation, identify periodic continuous change segments in time sequence, and obtain the response period continuous segment index sequence.

[0039] S503: Call the continuous segment index sequence of the response cycle, identify whether the rate response lag within the period segment has completed the return judgment within the material response window interval, filter the period segments that satisfy the synchronous closure of deformation duration and response lag time, mark the corresponding period closure information according to the number, and generate the nanopore on / off state transition trigger mark.

[0040] On the other hand, the nanopore functional control system is used to execute the above-described nanopore functional control method, the system comprising:

[0041] The path recognition module obtains the coordinates, liquid volume, contact interface length, and adhesion migration limit of the nanoporous liquid bridge anchor points in the porous array. It calculates the sorting sequence based on the length of the projected trajectory between the liquid bridge anchor point pairs, controls the arrangement order of the liquid migration path starting points based on the sorting sequence, and generates the initial dynamic configuration state of the liquid bridge network.

[0042] The interface response module extracts two-dimensional projection trajectory data from the edge of the liquid bridge to the nanopore based on the initial dynamic configuration state of the liquid bridge network, determines the direction of change between the trajectory endpoint and the displacement vector of the pore area, and generates a nanopore interface morphology response offset map.

[0043] The structural stabilization module calculates the position of the continuity inflection point within the curvature change range based on the nanopore interface morphology response offset spectrum, filters the boundary point between the structurally stable range and the unstable range, and obtains the nanopore boundary stability response signal.

[0044] The confinement correction module calls the boundary stability response signal of the nanopore, extracts the liquid coverage area value and tension path length value of the liquid bridge anchor point in the unstable region, adjusts the coverage area range according to a fixed expansion ratio, and generates a confinement correction field distribution map of the boundary region.

[0045] Based on the boundary region confined correction field distribution map, the thermal triggering module detects the temperature change rate value and material response time window value of the thermosensitive composite material, calculates the period ratio between the temperature change lag time and the orifice deformation time, and obtains the nanopore on / off state transition trigger mark.

[0046] The beneficial effects of the technical solutions provided in the embodiments of the present invention include at least the following:

[0047] By introducing sequential control and staggered behavior recognition of liquid migration paths, the dynamic configuration of liquid bridge distribution is optimized. Combined with interface response mapping of displacement peak index, the accuracy and efficiency of liquid position control are effectively improved. Through path and boundary response convergence extraction and periodic change analysis, the dynamic prediction of the stable state of 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 thermal triggering closure behavior judgment under temperature-sensitive conditions is achieved in the change of nanopore edge membrane thickness and regional ratio. The accurate marking and response period control of nanopore on / off state are completed, significantly improving the stability and controllability in molecular on / off control and structural dynamic recognition scenarios. Attached Figure Description

[0048] Figure 1 This is a schematic diagram of the workflow of the present invention;

[0049] Figure 2 This is a system flowchart of the present invention. Detailed Implementation

[0050] The technical solution of the present invention will now be described with reference to the accompanying drawings.

[0051] In embodiments of the present invention, words such as "exemplarily," "for example," etc., are used to indicate that something is an example, illustration, or description. Any embodiment or design described as "exemplary" in the present invention should not be construed as being more preferred or advantageous than other embodiments or designs. Specifically, the use of the word "exemplary" is intended to present the concept in a concrete manner. Furthermore, in embodiments of the present invention, the meaning expressed by "and / or" can be both, or either one.

[0052] To make the technical problems, technical solutions and advantages of the present invention clearer, a detailed description will be given below in conjunction with the accompanying drawings and specific embodiments.

[0053] Please see Figure 1 This invention provides a method for controlling the functionality of nanopores, the processing flow of which may include the following steps:

[0054] S1: Based on the initial distribution state of nanoporous 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 nanoporous liquid bridge anchor points are called and sorted, the starting order of liquid migration path is controlled, the intersecting behavior segments in the migration path are screened, and the initial dynamic configuration state of the liquid bridge network is generated.

[0055] S2: Based on the initial dynamic configuration of the liquid bridge network, extract the projection trajectory from the edge of the liquid bridge to the nanopore, identify the displacement change at the endpoint of the trajectory, calculate the peak displacement position index of the corresponding pore interface, perform position translation on the liquid distribution, and generate a nanopore interface morphology response offset map.

[0056] S3: Call the nanopore interface morphology response offset map, extract the pore edge curvature sequence, the liquid bridge path length change amplitude and the molecular contact frequency fluctuation range, perform continuous detection on the curvature change trend and extract periodic change points, perform projection sequence processing on the path length change amplitude, identify the segment that converges with the curvature change, and generate the nanopore boundary stability response signal.

[0057] S4: Based on the stability response signal of the nanopore boundary, extract the liquid coverage area and tension path value of the liquid bridge anchor point in the unstable region, increase the coverage area according to the set expansion ratio, and compress the conduction path to the constraint interval to generate a boundary region confined correction field distribution map.

[0058] S5: Call the boundary region confined correction field distribution map, detect the temperature change rate and material response window interval value of the covered thermosensitive composite material, extract the thickness of the pore edge film and the proportion of the variable region of the pore opening, compare the periodic continuity between the execution rate response lag and the thermal boundary deformation time, determine whether the continuous response behavior of the thermal triggering behavior has completed the closing condition, and generate the nanopore on / off state transition trigger mark.

[0059] The initial dynamic configuration 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 profile, displacement dynamic amplitude range, and trajectory terminal offset mapping result. The nanopore boundary stability response signal includes stability period identifier, path trend matching section indicator, and contact frequency fluctuation amplitude parameter. The boundary region confined correction field distribution map includes tension compression mapping map, coverage area expansion region, and confined field strength change block. The nanopore on / off state transition trigger marker includes trigger threshold index, response closure period, and on / off behavior identification site.

[0060] The specific steps for obtaining the initial dynamic configuration state of the liquid bridge network are as follows:

[0061] S101: Based on the initial distribution state of the anchor points of the nanoporous liquid bridge in the porous array, the liquid volume and contact interface length in the nanoporous 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 adhesion migration limit between each pair of anchor points is recorded, and a liquid bridge characteristic value table is generated.

[0062] A coordinate grid system for the entire array structure was established to identify and label the location of each nanopore. Based on this, for each pair of anchor points that could form a liquid bridge, the liquid volume and contact interface length within the liquid bridge were collected. This required high-resolution image capture to record the liquid bridge formation process frame-by-frame, followed by image processing tools (such as ImageJ or OpenCV) to extract the liquid bridge contours. The liquid volume was calculated using voxel integration, i.e., by calculating the volume projection of the liquid bridge in the image and performing volume conversion based on the array material thickness. The contact interface length was estimated using an edge detection algorithm (such as Canny edge detection) to obtain the contour lines and then 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 approximately elliptical, and the equivalent diameter of the liquid bridge after image recognition is 0.3. μm, 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 bivariate mapping datasets containing liquid volume (unit μm³) and corresponding contact interface length (unit μ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 by 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℃ and 50%RH, when the liquid bridge volume is 0.08μm³ and the interface length is 0.045μm, it can still adhere stably, while migration fails when it reaches 0.06μm³ and 0.041μm. Therefore, its migration limit is set to 0.07μm³ and 0.043μm and included in the numerical table to generate a liquid bridge characteristic numerical table.

[0063] S102: Call the liquid adhesion migration limit and contact interface length data between anchor point pairs in the liquid bridge feature value table, calculate the corresponding projection trajectory length based on the two-dimensional plane projection coordinates between adjacent nanopore liquid bridge anchor points, sort the anchor point pairs in ascending order, and eliminate paths that do not meet the migration conditions in the sorted sequence based on the liquid adhesion migration limit between anchor point pairs to obtain the migration path length sequence.

[0064] The geometric position is determined based on the two-dimensional projection coordinates of adjacent anchor points in the nanopore array. Anchor points (1, 1) and (1, 2) are set to have a distance of d = 0.5 μm in the array. The length of the projected trajectory is then... μm, this process requires traversing adjacent anchor point pairs, calculating the trajectory length between anchor point pairs that can form a liquid bridge 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 adhesion migration condition screening is performed on each anchor point pair. Condition judgment is performed on each pair of paths (i,j). If the migration limit of the anchor point pair in the feature value table (such as liquid bridge volume ≥ 0.07μm³, interface length ≥ 0.045μm) is not met, then this path is removed. The remaining paths, such as (1,1)-(1,2), (2,2)-(2,3), (2,3)-(2,4), etc., form a sequence of migrateable path lengths { , , During the screening process, two conditions need to be determined: first, whether the volume and interface of the actual liquid bridge formation meet the limit standard; second, whether the path sorting process conforms to the ascending order structure. For example, although the path (1,3)-(1,4) is short, its interface length is 0.035μm, which is lower than the limit of 0.043μm, so it is rejected to obtain the sequence of transferable path lengths.

[0065] S103: Based on the order and location distribution of the starting points of the paths in the sequence of migrated path lengths, control the starting points of the paths in sequence and track the staggered behavior of the path migration, identify the locations of overlapping sections in the paths and filter and mark them, and generate the initial dynamic configuration state of the liquid bridge network.

[0066] The starting point control strategy is set sequentially, starting from the edge nanopores. 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, it is necessary to analyze whether the path intersects segment by segment, that is, whether the two paths overlap in the two-dimensional projection. The overlap is judged based on whether the coordinates of the projected trajectory overlap in any segment. 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. Coordinate matching is performed on the path. If there is an intersection in the path segment... If a segment appears in both paths, it is considered 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), then it is marked as a segment. At the intersection, 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.

[0067] The specific steps for obtaining the nanopore interface morphology response shift spectrum are as follows:

[0068] S201: Based on the initial dynamic configuration of the liquid bridge network, extract the two-dimensional projection trajectory from the edge of the liquid bridge to the center of the corresponding nanopore, record the path coordinate sequence of the projection trajectory in the initial coordinate system, and map the trajectory to a unified coordinate system to generate the projection path sequence from the liquid bridge to the pore.

[0069] Based on the visualization and tracking of the formation and movement of liquid bridges in a porous array, a standard coordinate framework for the array is constructed in the initial state. A two-dimensional coordinate system is set with the central aperture of the array as the origin, and each aperture position corresponds to a unique coordinate number. The dynamic behavior of the liquid bridges is recorded throughout the process using a high-speed microscopic imaging device. The liquid bridge paths in the video frame sequence are extracted using image recognition tools. The movement trajectory of the liquid bridge formed between the anchor point and the aperture is extracted into a trajectory data stream composed of several path coordinate points. Each trajectory contains a two-dimensional movement path formed from the central anchor point to a specific aperture. The path data is transformed into trajectory path information within a unified coordinate system. By comparing the relative offset between the path points and the original anchor point positions in each frame, the mapping and unification of the path in the standard coordinate system is completed. This operation helps to achieve comparability of liquid bridge trajectories formed at different time points and in different experimental batches, laying the foundation for subsequent processing. When observing nanoporous materials with regular array structures, the two-dimensional trajectories formed by each liquid bridge path during movement can be classified, standardized, and archived to generate a sequence of projected paths from the liquid bridge to the aperture.

[0070] S202: Based on the projection path sequence from the liquid bridge to the orifice, identify the coordinate value of the trajectory endpoint on the surface of the nano-orifice, compare the endpoint coordinates of multiple time sections between the same anchor point, calculate the trajectory endpoint displacement equilibrium characteristic value, extract the peak change amount and record the coordinate index position, and generate the trajectory endpoint displacement index distribution value.

[0071] The trajectory endpoint displacement equilibrium characteristic value is calculated using the following formula:

[0072] ;

[0073] in, Represents anchor point pair , The trajectory endpoint displacement equilibrium characteristic value between them Represents anchor point In time section The horizontal endpoint coordinates, Represents anchor point In time section The horizontal endpoint coordinates, Represents anchor point In time section The vertical endpoint coordinates, Represents anchor point In time section The vertical endpoint coordinates, Represents anchor point In time section The local directional displacement increment, Represents anchor point , The average value of the local displacement increment over all time sections. Represents the number of time segments. The sequence index value representing the time section. , Represents the selected anchor point number;

[0074] The formula is used to calculate anchor point pairs. , The displacement equilibrium characteristic value of the endpoint coordinate sequence at different time sections, in this formula, the first term The second term represents the average squared difference between the anchor points in the horizontal coordinate direction, reflecting the consistency of the horizontal displacement between the two anchor points; The average value of the product of the displacement difference of the anchor point pair in the vertical coordinate direction and the local directional displacement increment deviation reflects the degree of coordinated change of the two anchor points in the vertical displacement. By summing the absolute values ​​of these two items and taking the square root, we obtain an index that comprehensively reflects the consistency of the two anchor points in spatial displacement.

[0075] The trajectory endpoint displacement equilibrium characteristic value is a numerical index that measures the consistency of the endpoint coordinate changes in space for the same anchor point pair across multiple time sections. It reflects the displacement differences and fluctuation coordination of the anchor point pair in the horizontal and vertical directions.

[0076] To further illustrate how each parameter is obtained, the following details the calculation process and example data:

[0077] , , , : Using high-precision image recognition equipment, at each time segment Above, extract anchor points and The horizontal and vertical endpoint coordinates in the time section At that time, anchor point The horizontal coordinate is The vertical coordinate is Anchor point The horizontal coordinate is The vertical coordinate is ;

[0078] : Indicates anchor point In time section The local displacement increment is obtained by calculating the vertical coordinate difference between adjacent time sections. and Above, anchor point The vertical coordinates are respectively and ,but ;

[0079] : Indicates anchor point and The average value of the local directional displacement increment on the time section is set as follows: At each time cross-section, the anchor point The local directional displacement increments are respectively , , ,but ;

[0080] : Represents the total number of time segments, determined through experimental design. If data from 3 time segments are collected in the experiment, then ;

[0081] Substitute the example data above into the formula to calculate:

[0082] ;

[0083] ;

[0084] Comprehensive calculation:

[0085] ;

[0086] This result indicates that the anchor point and The displacement equilibrium characteristic value within the selected time section is This indicates a high degree of consistency between the two anchor points in spatial displacement;

[0087] The advantage of the 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 characteristic value.

[0088] To further illustrate the parameter acquisition process, the experimental data are listed below:

[0089] Table 1: Experimental Data Table

[0090]

[0091] Table 1 lists the anchor points at three time sections. and The horizontal and vertical endpoint coordinates, and the anchor point The local directional displacement increment can be used to calculate the trajectory endpoint displacement equilibrium characteristic value using the above formula.

[0092] S203: Based on the corresponding index position in the trajectory endpoint displacement index distribution value, map the spatial distribution area of ​​the original liquid around the nanopore, perform a uniform translation of the liquid distribution area in the two-dimensional coordinate system, determine the translation direction according to the direction indicated by the displacement index, determine the translation amplitude using the displacement increment, and generate a liquid distribution translation coordinate set.

[0093] The positional offset of the liquid bridge endpoint under 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 transferred as a whole to construct a new translational liquid region coordinate set. In the operation, each endpoint offset in the displacement index needs to be read one by one, and the spatial direction and movement amplitude represented are applied to the corresponding liquid distribution pattern. The liquid region around each pore is spatially displaced according to this transformation to obtain a new liquid region layout. The pore in the 5th row and 4th column of the array is set. Its liquid region is concentrated in the center in the original state, but after the trajectory endpoint displacement transformation, it is offset to the edge of the pore. This translation transformation needs to be accurately recorded and summarized to form a coordinate set. This coordinate set 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.

[0094] S204: Call the liquid distribution translation coordinate set, reorganize the liquid distribution area according to the corresponding nanopore number, and reconstruct the interface contour graphic under the standard pore contour boundary. Plot the boundary contour area corresponding to the difference in liquid distribution before and after displacement to obtain the nanopore interface morphology response offset map.

[0095] Each nanopore is numbered and labeled accordingly. The liquid distribution areas before and after migration are reconstructed, and the interface contour is reconstructed by referring to the standard pore contour. The liquid distribution in the original state is spatially matched with its corresponding translational state, and the boundary contours of the two states are drawn. By calculating the difference between the boundary line envelope and the outer contour, the interface morphology boundary change area caused by the change in liquid distribution can be generated. In the process of graphic reconstruction, the liquid boundary contours before and after migration need to be plotted separately, and the difference area is highlighted by filling. By using parameters such as the degree of boundary layer movement, direction and the degree of contour overlap, the liquid movement process experienced by each nanopore is fully characterized. In actual observation, some pores can be found to form asymmetric migration patterns due to uneven surface structure. Such differences can be identified and labeled by the migration spectrum. A complete pore response graphic structure set is constructed to obtain the nanopore interface morphology response migration spectrum.

[0096] The specific steps for obtaining the stability response signal at the nanopore boundary are as follows:

[0097] S301: Call the nanopore interface morphology response offset map, extract the continuous values ​​of the pore edge curvature along the boundary coordinate direction, and collect the length change value of the corresponding liquid bridge path and the molecular contact frequency sequence of the contact area at the end of the liquid bridge. Record the numerical fluctuation range on the time axis to generate the curvature path frequency joint sequence.

[0098] Curvature values ​​at the orifice edge are extracted under a unified coordinate system. Sampling is performed point-by-point along the boundary coordinate direction, and the boundary curvature values ​​at each coordinate point are continuously recorded using contour tracing technology. These values ​​are then combined to form a complete edge curvature sequence. Simultaneously, the length changes of each liquid bridge path at different time periods are extracted from the normalized liquid bridge trajectory. By summarizing the coordinate changes between the path's start and end points, the liquid bridge path length change value is obtained, and the numerical fluctuations over time are recorded. The contact behavior of molecules in the contact area between the liquid bridge end and the nanopore surface is statistically analyzed. The contact frequency change sequence per unit time is recorded using orbital tracking algorithms or molecular dynamics simulation methods, gradually accumulating to form a contact frequency dataset that changes over time. The data is synchronized with the curvature sequence and path change values, categorized by time axis, and the numerical fluctuation range of each parameter at different time points is identified. Finally, the edge curvature, path length changes, and molecular contact frequency data are combined to generate a joint curvature-path-frequency sequence.

[0099] S302: Based on the hole edge curvature sequence in the joint sequence of curvature path frequency, the curvature value change relationship between continuous coordinate points is detected in the order of boundary coordinates. For the position where adjacent positive and negative incremental directions change, the periodic change point is extracted, and the index segment where the curvature trend reverses is marked to obtain the edge curvature change index set.

[0100] The curvature value change trend between consecutive coordinate points is checked sequentially along the boundary contour coordinates. The positive and negative increments of the curvature value between adjacent points are recorded. If a certain segment of curvature value shows a continuous increasing trend and then starts to decrease at a certain point, it indicates that this point is a trend reversal node, i.e., a periodic change point. In this operation, the boundary coordinate points are scanned in numerical order, and the difference of each pair of consecutive curvature values ​​is judged and its change direction is recorded. When a point is opposite to the previous and subsequent change direction, the point is marked as the trend reversal position, and its index number is marked in the boundary sequence. At the same time, consecutive reversal points are connected into index segments. The 35th point in a certain orifice boundary is set as the critical point from increasing to decreasing, and the 46th point is set as the point from decreasing to increasing. This indicates that 35 to 46 is a complete curvature change segment. The index segment represents the local contour area where the curvature change trend is reversed. The area often corresponds to the concentrated part of the orifice micro-deformation and plays a key role in the subsequent path and molecular behavior analysis, thus obtaining the edge curvature change index set.

[0101] S303: Call the edge curvature change index to concentrate the liquid bridge path length change range, rearrange the coordinates and normalize the direction of the path change range according to the index segment, select the path segment with the same change trend as the convergence segment, and count the difference range of molecular contact frequency fluctuation in the corresponding segment to obtain the nanopore boundary stability response signal.

[0102] The algorithm retrieves the data on the change in length of the liquid bridge path within each direction-changing segment and rearranges the spatial coordinates of the path change values ​​according to the range of the curvature change segment. This involves adjusting the start and end coordinates of the path to match the boundary index, unifying the starting point for path and curvature change analysis. Simultaneously, it normalizes the direction of path length change, uniformly labeling data in increasing or decreasing directions for easier comparison and filtering later. It also determines the consistency between the path change direction and the trend of the direction-changing segment. If a segment where curvature changes from increasing to decreasing is accompanied by a change in path length from shortening to increasing, then that segment is considered to have a consistent curvature change trend and is selected as a convergent segment. Within the convergent segment, the difference range of molecular contact frequencies is extracted, i.e., the extreme difference of the frequency sequence is extracted within several frames before and after the change. These differences are then categorized and statistically analyzed to obtain the molecular contact intensity response range for each convergent segment. This is used for interface state modeling analysis and anomaly warning identification, and to obtain the nanopore boundary stability response signal.

[0103] The specific steps for obtaining the boundary region confined correction field distribution map are as follows:

[0104] S401: Based on the stability response signal of the nanopore boundary, extract the liquid bridge anchor point number marked as the unstable section, record the liquid coverage area under the corresponding anchor point and the tension transmission path length to the pore boundary, and generate a liquid bridge anchor point coverage tension value group.

[0105] The process involves clearly identifying the unstable segments, which correspond to areas exhibiting structural anomalies or drastic fluctuations in contact frequency during the formation of specific anchor points in the liquid bridge. The signal sequence must be analyzed segment by segment to extract the anchor point numbers marked as unstable, creating an index list of these anchor points. Based on this, the area covered by each unstable anchor point is calculated. This operation relies on the closed shape formed by the liquid region's contour boundary, obtaining the liquid coverage area value through contour tracking and pixel area conversion. The length of the tension transmission path from the anchor point to the orifice boundary is recorded, specifically by calculating the actual path length along the shortest continuous channel between the liquid bridge and the orifice boundary. The path length extraction requires reference to the liquid bridge's projection trajectory in the three-dimensional structure and estimation based on the path projection in two-dimensional coordinates. For each unstable anchor point, both the corresponding liquid coverage area and tension path length are simultaneously obtained. These two data points, along with the anchor point numbers, are categorized and stored to generate a group of liquid bridge anchor point coverage tension values.

[0106] S402: Call the tension value group of the liquid bridge anchor point coverage, perform numerical amplification on the coverage area according to the set expansion ratio, calculate the characteristic value of the amplified area, 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, redraw the tension distribution area outline of the anchor point position in the boundary area, and obtain the boundary area confinement correction field distribution map.

[0107] The characteristic value of the area of ​​increase is calculated using the following formula:

[0108] ;

[0109] in, The characteristic value representing the increase in area, This represents the initial area covered by the basic liquid bridge anchor point. Representing the Anchor point tension conduction gain factor Representing the Anchor tension value, Representing the Anchor point path length compression damping coefficient Representing the Item anchor path length, This represents the equivalent radius length formed by the mapping of the tension path in a two-dimensional plane. This represents the average conducted interference displacement correction value. Represents the number of anchor points;

[0110] In this formula, the base area The projected area of ​​the liquid bridge anchor point before expansion is represented by the area calculated from the measured data of the liquid bridge morphology using optical imaging combined with image segmentation; tension value. Tension fluctuations are measured using a miniature tension sensor embedded in the anchor point structure, and the average value of tension fluctuations under different working conditions is collected and processed through amplitude equalization; path length. The linear distance between points is obtained along the tension direction using a laser displacement meter; coefficient , The tension gain and damping compression weighting coefficients represent the positive driving force of tension on area expansion and the negative compression weight of path length on area increase, respectively. These values ​​are set based on process test calibration values, with the following reference: 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 contraction amplitude under changes in path length (to obtain the compression coefficient) The linear fitting coefficients are used as the source of coefficient values. The displacement correction for tension transmission disturbance is obtained by quantifying the positional changes of the anchor point after being disturbed by the surrounding environment in multiple experiments. It is obtained by averaging the squared difference between the maximum displacement under disturbance and its static mean value. The radius of the tension path is calculated by spatial mapping at the end of the tension path at the anchor point, which forms the radius of the circumcircle.

[0111] The amplified area characteristic value refers to the effective coverage area under the tension expansion effect of the anchor point. It is the area correction result obtained after considering factors such as tension intensity, 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 effect.

[0112] Parameter acquisition and practical examples:

[0113] Taking experimental sample A as an example, its basic area was extracted by image recognition in its initial state. ;

[0114] There are a total of 4 anchor points, and the number of effective transmission anchor points is [number missing]. One anchor point was rejected due to insufficient tension;

[0115] The tension values ​​measured by the tension sensor are as follows: ;

[0116] The path lengths were measured as follows: ;

[0117] The weighting coefficients determined in the experiment are as follows:

[0118] (Fitted by tension influence rate);

[0119] (Derived from normalized path compression sensitivity experiments);

[0120] Disturbance correction value The maximum offset of each anchor point was obtained by monitoring with a laser interferometer and then calculated. ;

[0121] Path projection circumcircle radius ;

[0122] Substituting each term into the formula, the calculation is as follows:

[0123] ;

[0124] ;

[0125] ;

[0126] Summing yields:

[0127] ;

[0128] Substitute into the formula:

[0129] ;

[0130] The results show that the area feature value is increased to 73.7 μm², and the area expansion, after correction by the tension-path joint parameters, has reached 5.9 times the base area; this area feature value will be used for boundary redrawing of the tension distribution area in subsequent steps.

[0131] The advantage of the formula is that by introducing the absolute difference between the tension gain coefficient and the path compression coefficient to construct an area contribution term, and combining the tension path space projection and the disturbance correction term to perform normalization square root correction, the area increase processing not only reflects the difference in the actual tension distribution, but also effectively integrates the path complexity and the offset caused by the disturbance, providing more reasonable data support for subsequent boundary redrawing.

[0132] Table 2: Monitoring Table of Area Increase Parameters

[0133]

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

[0135] The specific steps for obtaining the on / off state transition trigger marker for nanopores are as follows:

[0136] S501: Call the boundary region confined correction field distribution map, monitor the time series temperature change curve of the area covered by the thermosensitive composite material, calculate the temperature change rate per unit time of the area, extract the window interval value simultaneously, organize the temperature rate and response window value pairs according to the coordinate number, and generate a thermosensitive response rate window pairing table.

[0137] It is necessary to identify the thermosensitive composite material area covered in the map and perform continuous temperature change time series monitoring on it. Non-contact high-resolution infrared thermal imaging is used to record the temperature field of the material surface in frames, and 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 position, the temperature change rate per unit time is calculated by dividing the temperature difference between continuous time points by the time interval to obtain the temperature rate value of each region. Then, the temperature rate sequence is segmented and extracted according to the preset time window parameters. The extracted values ​​are used as representative responses within the window interval. A combined dataset is established with coordinate number as index, temperature rate as the first value, and temperature change within the window interval as the second value, generating a thermosensitive response rate window pairing table.

[0138] S502: Based on the thermosensitive response rate window pairing table, extract the film thickness data at the nanopore edge position and the variable region ratio data in the pore area, identify the corresponding structure between the pore edge structure size and the response capability, compare the rate response lag time and the duration of thermal boundary deformation, identify periodic continuous change segments in time sequence, and obtain the response period continuous segment index sequence.

[0139] The film thickness information at the nanopore edge positions corresponding to the pore region is extracted. Nanoscale tomography or atomic force microscopy is used to acquire the film structure size data at each pore edge position. Combined with the ratio of the variable region to the overall region in the pore region under different thermal response states, a corresponding structure table of pore edge size and response capability is constructed. The correlation between structure and rate response data is compared, focusing on the parts with thicker films or a high proportion of variable regions to check for significant shifts in their thermal response rate and hysteresis characteristics. The rate response hysteresis time (the time interval between the thermal stimulus input and the material's significant temperature response) is retrieved from the thermosensitive response curve. The duration of each thermal interface deformation is recorded, and both are used as key response parameters, arranged in chronological order. The hysteresis time is compared to see if it is shorter or longer than the response duration. Segments with repetitive rhythmic changes are retrieved to form continuous segments with periodic thermal response characteristics. Each periodic segment is indexed and numbered, and the continuous segment index sequence of the response cycle is output.

[0140] S503: Call the continuous segment index sequence of the response cycle, identify whether the rate response lag within the period segment has completed the return judgment within the material response window interval, filter the period segments that meet the synchronous closure of deformation duration and response lag time, mark the corresponding period closure information by number, and generate nanopore on / off state transition trigger markers.

[0141] It is necessary to determine whether the cycle meets the condition of rate response lag completion within the cycle segment. That is, within the complete cycle, 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. In operation, it is necessary to compare with the thermosensitive response window range to determine whether the lag time is always contained within the material response window. If the lag behavior within the cycle exceeds the response range, it is judged as an unclosed cycle segment; otherwise, it is a synchronously closed cycle segment. In the cycle segment, the cycle segments that meet the synchronous closure of lag completion and deformation duration are selected and their cycle numbers are specially marked. A closure mark is added after the cycle segment to form a complete nanopore response cycle labeling system. This provides a control signal for each cycle that meets the material thermal response stable closure condition, provides a logical basis and time node basis for the opening and closing state switching of the pore under micro-thermal drive, and generates a nanopore on / off state transition trigger mark.

[0142] Please see Figure 2 A nanopore functional control system, the system comprising:

[0143] The path recognition module obtains the coordinates, liquid volume, contact interface length, and adhesion migration limit of the nanoporous liquid bridge anchor points in the porous array. It calculates the sorting sequence based on the length of the projected trajectory between the liquid bridge anchor point pairs, controls the arrangement order of the liquid migration path starting points based on the sorting sequence, and generates the initial dynamic configuration state of the liquid bridge network.

[0144] The interface response module extracts the two-dimensional projection trajectory data from the edge of the liquid bridge to the nanopore based on the initial dynamic configuration state of the liquid bridge network, determines the direction of change between the trajectory endpoint and the displacement vector of the pore area, and generates a nanopore interface morphology response offset map.

[0145] The structural stabilization module calculates the position of the continuity inflection point within the curvature change range based on the morphological response offset spectrum of the nanopore interface, selects the boundary point between the structurally stable range and the unstable range, and obtains the nanopore boundary stability response signal.

[0146] The confinement correction module calls the boundary stability response signal of the nanopore, extracts the liquid coverage area value and tension path length value of the liquid bridge anchor point in the unstable region, adjusts the coverage area range according to a fixed expansion ratio, and generates a confinement correction field distribution map of the boundary region.

[0147] The thermal triggering module detects the temperature change rate and material response time window value of the thermosensitive composite material based on the boundary region confined correction field distribution map, calculates the period ratio between the temperature change lag time and the orifice deformation time, and obtains the nanopore on / off state transition trigger mark.

[0148] 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 variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in the present invention should be included within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.

Claims

1. A method for controlling the functionality of nanopores, characterized in that, Includes the following steps: S1: Based on the initial distribution state of nanoporous liquid bridge anchor points in the porous array, the liquid volume, contact interface length and adhesion migration limit are collected, the intersecting 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, extract the projection trajectory from the edge of the liquid bridge to the nanopore, identify the displacement change at the endpoint of the trajectory, perform position translation on the liquid distribution, and generate a nanopore interface morphology response offset map. S3: Call the nanopore interface morphology response offset map, extract the pore edge curvature sequence, the liquid bridge path length change amplitude and the molecular contact frequency fluctuation range, perform continuous detection on the curvature change trend and extract periodic change points, perform projection sequence organization on the path length change amplitude, identify the segment that is similar to the curvature change, and generate the nanopore boundary stability response signal. S4: Based on the stability response signal of the nanopore boundary, extract the liquid coverage area and tension path value of the liquid bridge anchor point in the unstable region, increase the coverage area according to the set expansion ratio, and compress the conduction path to the constraint interval to generate a boundary region confined correction field distribution map.

2. The method for controlling the functionality of nanopores according to claim 1, characterized in that, 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 profile, displacement dynamic amplitude range, and trajectory terminal offset mapping result. The nanopore boundary stability response signal includes stability period identifier, path trend matching section indicator, and contact frequency fluctuation amplitude parameter. The boundary region confined correction field distribution map includes tension compression mapping map, coverage area expansion region, and confined field strength change block.

3. The method for controlling the functionality of nanopores according to claim 1, characterized in that, The specific steps for obtaining the initial dynamic configuration state of the liquid bridge network are as follows: S101: Based on the initial distribution state of the anchor points of the nanoporous liquid bridge in the porous array, the liquid volume and contact interface length in the nanoporous 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 adhesion migration limit between each pair of anchor points is recorded, and a liquid bridge characteristic value table is generated. S102: Call the liquid adhesion migration limit and contact interface length data between anchor point pairs in the liquid bridge feature value table, calculate the corresponding projection trajectory length based on the two-dimensional plane projection coordinates between adjacent nanopore liquid bridge anchor points, sort the anchor point pairs in ascending order, and eliminate paths in the sorted sequence that do not meet the migration conditions by combining the liquid adhesion migration limit between anchor point pairs, and obtain the migration path length sequence. S103: Based on the path starting point order and position distribution in the migrated path length sequence, control the path starting point in sequence and track the staggered behavior state of path migration, identify the position segments of the path where there are intersecting and overlapping sections and filter and mark them to generate the initial dynamic configuration state of the liquid bridge network.

4. The method for controlling the functionality of nanopores according to claim 1, characterized in that, The specific steps for obtaining the nanopore interface morphology response shift map are as follows: S201: Based on the initial dynamic configuration state of the liquid bridge network, extract the two-dimensional projection trajectory from the edge of the liquid bridge to the center of the corresponding nanopore, record the path coordinate sequence of the projection trajectory in the initial coordinate system, and map the trajectory to a unified coordinate system to generate the projection path sequence from the liquid bridge to the pore. S202: Based on the projection path sequence from the liquid bridge to the orifice, identify the coordinate value of the trajectory endpoint on the surface of the nano-orifice, compare the endpoint coordinates of multiple time sections between the same anchor point, calculate the trajectory endpoint displacement equilibrium characteristic value, extract the peak change amount and record the coordinate index position, and generate the trajectory endpoint displacement index distribution value. The trajectory endpoint displacement equalization characteristic value is obtained using the formula: ； in, Represents anchor point pair , The trajectory endpoint displacement equilibrium characteristic value between them Represents anchor point In time section The horizontal endpoint coordinates, Represents anchor point In time section The horizontal endpoint coordinates, Represents anchor point In time section The vertical endpoint coordinates, Represents anchor point In time section The vertical endpoint coordinates, Represents anchor point In time section The local directional displacement increment, Represents anchor point , The average value of the local displacement increment over all time sections. Represents the number of time segments. The sequence index value representing the time section. , Represents the selected anchor point number; S203: Based on the corresponding index position in the displacement index distribution value of the trajectory endpoint, map the spatial distribution area of ​​the original liquid around the nanopore, perform a uniform translation of the liquid distribution area in the two-dimensional coordinate system, determine the translation direction according to the direction indicated by the displacement index, determine the translation amplitude using the displacement increment, and generate a liquid distribution translation coordinate set. S204: Call the liquid distribution translation coordinate set, reorganize the liquid distribution area according to the corresponding nanopore number, reconstruct the interface contour graphic under the standard pore contour boundary, plot the boundary contour area corresponding to the difference in liquid distribution before and after displacement, and obtain the nanopore interface morphology response offset map.

5. The method for controlling the functionality of nanopores according to claim 4, characterized in that, The specific steps for obtaining the nanopore boundary stability response signal are as follows: S301: Call the nanopore interface morphology response offset map, extract the continuous values ​​of the pore edge curvature along the boundary coordinate direction, and collect the length change value of the corresponding liquid bridge path and the molecular contact frequency sequence of the contact area at the end of the liquid bridge. Record the numerical fluctuation range on the time axis to generate the curvature path frequency joint sequence. S302: Based on the hole edge curvature sequence in the curvature path frequency joint sequence, detect the curvature value change relationship between continuous coordinate points in the boundary coordinate order, extract the periodic change point at the position where the adjacent positive and negative incremental directions change, mark the index segment where the curvature trend reverses, and obtain the edge curvature change index set. S303: Call the edge curvature direction change index to concentrate the liquid bridge path length change range, rearrange the coordinates and normalize the direction of the path change range according to the index segment, select the path segment with the same change trend as the convergence segment, and count the difference range of molecular contact frequency fluctuation in the corresponding segment to obtain the nanopore boundary stability response signal.

6. The method for controlling the functionality of nanopores according to claim 5, characterized in that, The specific steps for obtaining the boundary region confined correction field distribution map are as follows: S401: Based on the stability response signal of the nanopore boundary, extract the liquid bridge anchor point number marked as the unstable section, record the liquid coverage area under the corresponding anchor point and the tension transmission path length to the pore boundary, and generate a liquid bridge anchor point coverage tension value group. S402: Call the tension value group of the liquid bridge anchor point coverage, perform numerical amplification processing on the coverage area according to the set expansion ratio, calculate the characteristic value of the amplified area, 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, redraw the tension distribution area outline of the anchor point position in the boundary area, and obtain the boundary area confinement correction field distribution map.

7. The method for controlling the functionality of nanopores according to claim 6, characterized in that, The characteristic value of the increased area is expressed by the formula: ； in, The characteristic value representing the increase in area, This represents the initial area covered by the basic liquid bridge anchor point. Representing the Anchor point tension conduction gain factor Representing the Anchor tension value, Representing the Anchor point path length compression damping coefficient Representing the Item anchor path length, This represents the equivalent radius length formed by the mapping of the tension path in a two-dimensional plane. This represents the average conducted interference displacement correction value. This represents the number of anchor points.

8. The method for controlling the functionality of nanopores according to claim 1, characterized in that, The method further includes step S5: S5: Call the boundary region confined correction field distribution map, detect the temperature change rate and material response window interval value of the covered thermosensitive composite material, extract the thickness of the pore edge film and the proportion of the variable region of the pore opening, compare the periodic continuity between the execution rate response lag and the thermal boundary deformation time, determine whether the continuous response behavior of the thermal triggering behavior has completed the closing condition, and generate a nanopore on / off state transition trigger mark. The nanopore on / off state transition trigger marker includes a trigger threshold index, a response closure period, and an on / off behavior identification site.

9. The method for controlling the functionality of nanopores according to claim 8, characterized in that, The specific steps for obtaining the on / off state transition trigger marker of the nanopore are as follows: S501: Call the boundary region confined correction field distribution map, monitor the time series temperature change curve of the area covered by the thermosensitive composite material, calculate the temperature change rate per unit time of the area, extract the window interval value synchronously, organize the temperature rate and response window value pair according to the coordinate number, and generate a thermosensitive response rate window pairing table. S502: Based on the thermosensitive response rate window pairing table, extract the film thickness data at the nanopore edge position and the variable region ratio data in the pore area, identify the corresponding structure between the pore edge structure size and the response capability, compare the rate response lag time and the duration of thermal boundary deformation, identify periodic continuous change segments in time sequence, and obtain the response period continuous segment index sequence. S503: Call the continuous segment index sequence of the response cycle, identify whether the rate response lag within the period segment has completed the return judgment within the material response window interval, filter the period segments that satisfy the synchronous closure of deformation duration and response lag time, mark the corresponding period closure information according to the number, and generate the nanopore on / off state transition trigger mark.

10. A nanopore functional control system, characterized in that, The system is used to implement the nanopore functional control method according to any one of claims 1-9, the system comprising: The path recognition module obtains the coordinates, liquid volume, contact interface length, and adhesion migration limit of the nanoporous liquid bridge anchor points in the porous array. It calculates the sorting sequence based on the length of the projected trajectory between the liquid bridge anchor point pairs, controls the arrangement order of the liquid migration path starting points based on the sorting sequence, and generates the initial dynamic configuration state of the liquid bridge network. The interface response module extracts two-dimensional projection trajectory data from the edge of the liquid bridge to the nanopore based on the initial dynamic configuration state of the liquid bridge network, determines the direction of change between the trajectory endpoint and the displacement vector of the pore area, and generates a nanopore interface morphology response offset map. The structural stabilization module calculates the position of the continuity inflection point within the curvature change range based on the nanopore interface morphology response offset spectrum, filters the boundary point between the structurally stable range and the unstable range, and obtains the nanopore boundary stability response signal. The confinement correction module calls the boundary stability response signal of the nanopore, extracts the liquid coverage area value and tension path length value of the liquid bridge anchor point in the unstable region, adjusts the coverage area range according to a fixed expansion ratio, and generates a confinement correction field distribution map of the boundary region. Based on the boundary region confined correction field distribution map, the thermal triggering module detects the temperature change rate value and material response time window value of the thermosensitive composite material, calculates the period ratio between the temperature change lag time and the orifice deformation time, and obtains the nanopore on / off state transition trigger mark.

Images