Modularized and Parametric Modeling Technology for Finite Element Simulations of Underground Engineering under Complicated Geological Conditions

The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling efficiency of underground engineering,a modularized and parametric modeling cloud server is developed by using Python codes.The basic framework of the cloud server is as follows:input the modeling parameters into the web platform,implement Rhino software and FLAC3D software to model and run simulations in the cloud server,and return the simulation results to the web platform.The modeling program can automatically generate instructions that can run the modeling process in Rhino based on the input modeling parameters.The main modules of the modeling program include modeling the 3D geological structures,the underground engineering structures,and the supporting structures as well as meshing the geometric models.In particular,various cross-sections of underground caverns are crafted as parametricmodules in themodeling program.Themodularized and parametric modeling program is used for a finite element simulation of the underground powerhouse of the Shuangjiangkou Hydropower Station.This complicatedmodel is rapidly generated for the simulation,and the simulation results are reasonable.Thus,this modularized and parametric modeling program is applicable for three-dimensional finite element simulations and analyses.

Mathematical Modeling of HIV Investigating the Effect of Inconsistent Treatment

HIV is a retrovirus that infects and impairs the cells and functions of the immune system. It has caused a great challenge to global public health systems and leads to Acquired Immunodeficiency Syndrome (AIDS), if not attended to in good time. Antiretroviral therapy is used for managing the virus in a patient’s lifetime. Some of the symptoms of the disease include lean body mass and many opportunistic infections. This study has developed a SIAT mathematical model to investigate the impact of inconsistency in treatment of the disease. The arising non-linear differential equations have been obtained and analyzed. The DFE and its stability have been obtained and the study found that it is locally asymptotically stable when the basic reproduction number is less than unity. The endemic equilibrium has been obtained and found to be globally asymptotically stable when the basic reproduction number is greater than unity. Numerical solutions have been obtained and analyzed to give the trends in the spread dynamics. The inconsistency in treatment uptake has been analyzed through the numerical solutions. The study found that when the treatment rate of those infected increases, it leads to an increase in treatment population, which slows down the spread of HIV and vice versa. An increase in the rate of treatment of those with AIDS leads to a decrease in the AIDS population, the reverse happens when this rate decreases. The study recommends that the community involvement in advocating for consistent treatment of HIV to curb the spread of the disease.

Mathematical Modeling of the Co-Infection Dynamics of HIV and Tuberculosis Incorporating Inconsistency in HIV Treatment

A non-linear HIV-TB co-infection has been formulated and analyzed. The positivity and invariant region has been established. The disease free equilibrium and its stability has been determined. The local stability was determined and found to be stable under given conditions. The basic reproduction number was obtained and according to findings, co-infection diminishes when this number is less than unity, and persists when the number is greater than unity. The global stability of the endemic equilibrium was calculated. The impact of HIV on TB was established as well as the impact of TB on HIV. Numerical solution was also done and the findings indicate that when the rate of HIV treatment increases the latent TB increases while the co-infected population decreases. When the rate of HIV treatment decreases the latent TB population decreases and the co-infected population increases. Encouraging communities to prioritize the consistent treatment of HIV infected individuals must be emphasized in order to reduce the scourge of HIV-TB co-infection.

Predictive modeling for postoperative delirium in elderly patients with abdominal malignancies using synthetic minority oversampling technique

BACKGROUND Postoperative delirium,particularly prevalent in elderly patients after abdominal cancer surgery,presents significant challenges in clinical management.AIM To develop a synthetic minority oversampling technique(SMOTE)-based model for predicting postoperative delirium in elderly abdominal cancer patients.METHODS In this retrospective cohort study,we analyzed data from 611 elderly patients who underwent abdominal malignant tumor surgery at our hospital between September 2020 and October 2022.The incidence of postoperative delirium was recorded for 7 d post-surgery.Patients were divided into delirium and non-delirium groups based on the occurrence of postoperative delirium or not.A multivariate logistic regression model was used to identify risk factors and develop a predictive model for postoperative delirium.The SMOTE technique was applied to enhance the model by oversampling the delirium cases.The model’s predictive accuracy was then validated.RESULTS In our study involving 611 elderly patients with abdominal malignant tumors,multivariate logistic regression analysis identified significant risk factors for postoperative delirium.These included the Charlson comorbidity index,American Society of Anesthesiologists classification,history of cerebrovascular disease,surgical duration,perioperative blood transfusion,and postoperative pain score.The incidence rate of postoperative delirium in our study was 22.91%.The original predictive model(P1)exhibited an area under the receiver operating characteristic curve of 0.862.In comparison,the SMOTE-based logistic early warning model(P2),which utilized the SMOTE oversampling algorithm,showed a slightly lower but comparable area under the curve of 0.856,suggesting no significant difference in performance between the two predictive approaches.CONCLUSION This study confirms that the SMOTE-enhanced predictive model for postoperative delirium in elderly abdominal tumor patients shows performance equivalent to that of traditional methods,effectively addressing data imbalance.

Background removal from global auroral images:Data-driven dayglow modeling

Global images of auroras obtained by cameras on spacecraft are a key tool for studying the near-Earth environment.However,the cameras are sensitive not only to auroral emissions produced by precipitating particles,but also to dayglow emissions produced by photoelectrons induced by sunlight.Nightglow emissions and scattered sunlight can contribute to the background signal.To fully utilize such images in space science,background contamination must be removed to isolate the auroral signal.Here we outline a data-driven approach to modeling the background intensity in multiple images by formulating linear inverse problems based on B-splines and spherical harmonics.The approach is robust,flexible,and iteratively deselects outliers,such as auroral emissions.The final model is smooth across the terminator and accounts for slow temporal variations and large-scale asymmetries in the dayglow.We demonstrate the model by using the three far ultraviolet cameras on the Imager for Magnetopause-to-Aurora Global Exploration(IMAGE)mission.The method can be applied to historical missions and is relevant for upcoming missions,such as the Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)mission.

Surrogate modeling for unsaturated infiltration via the physics and equality-constrained artificial neural networks

Machine learning(ML)provides a new surrogate method for investigating groundwater flow dynamics in unsaturated soils.Traditional pure data-driven methods(e.g.deep neural network,DNN)can provide rapid predictions,but they do require sufficient on-site data for accurate training,and lack interpretability to the physical processes within the data.In this paper,we provide a physics and equalityconstrained artificial neural network(PECANN),to derive unsaturated infiltration solutions with a small amount of initial and boundary data.PECANN takes the physics-informed neural network(PINN)as a foundation,encodes the unsaturated infiltration physical laws(i.e.Richards equation,RE)into the loss function,and uses the augmented Lagrangian method to constrain the learning process of the solutions of RE by adding stronger penalty for the initial and boundary conditions.Four unsaturated infiltration cases are designed to test the training performance of PECANN,i.e.one-dimensional(1D)steady-state unsaturated infiltration,1D transient-state infiltration,two-dimensional(2D)transient-state infiltration,and 1D coupled unsaturated infiltration and deformation.The predicted results of PECANN are compared with the finite difference solutions or analytical solutions.The results indicate that PECANN can accurately capture the variations of pressure head during the unsaturated infiltration,and present higher precision and robustness than DNN and PINN.It is also revealed that PECANN can achieve the same accuracy as the finite difference method with fewer initial and boundary training data.Additionally,we investigate the effect of the hyperparameters of PECANN on solving RE problem.PECANN provides an effective tool for simulating unsaturated infiltration.

A methodology for damage evaluation of underground tunnels subjected to static loading using numerical modeling

We have proposed a methodology to assess the robustness of underground tunnels against potential failure.This involves developing vulnerability functions for various qualities of rock mass and static loading intensities.To account for these variations,we utilized a Monte Carlo Simulation(MCS)technique coupled with the finite difference code FLAC^(3D),to conduct two thousand seven hundred numerical simulations of a horseshoe tunnel located within a rock mass with different geological strength index system(GSIs)and subjected to different states of static loading.To quantify the severity of damage within the rock mass,we selected one stress-based(brittle shear ratio(BSR))and one strain-based failure criterion(plastic damage index(PDI)).Based on these criteria,we then developed fragility curves.Additionally,we used mathematical approximation techniques to produce vulnerability functions that relate the probabilities of various damage states to loading intensities for different quality classes of blocky rock mass.The results indicated that the fragility curves we obtained could accurately depict the evolution of the inner and outer shell damage around the tunnel.Therefore,we have provided engineers with a tool that can predict levels of damages associated with different failure mechanisms based on variations in rock mass quality and in situ stress state.Our method is a numerically developed,multi-variate approach that can aid engineers in making informed decisions about the robustness of underground tunnels.

Artificial intelligence-driven radiomics study in cancer:the role of feature engineering and modeling

Modern medicine is reliant on various medical imaging technologies for non-invasively observing patients’anatomy.However,the interpretation of medical images can be highly subjective and dependent on the expertise of clinicians.Moreover,some potentially useful quantitative information in medical images,especially that which is not visible to the naked eye,is often ignored during clinical practice.In contrast,radiomics performs high-throughput feature extraction from medical images,which enables quantitative analysis of medical images and prediction of various clinical endpoints.Studies have reported that radiomics exhibits promising performance in diagnosis and predicting treatment responses and prognosis,demonstrating its potential to be a non-invasive auxiliary tool for personalized medicine.However,radiomics remains in a developmental phase as numerous technical challenges have yet to be solved,especially in feature engineering and statistical modeling.In this review,we introduce the current utility of radiomics by summarizing research on its application in the diagnosis,prognosis,and prediction of treatment responses in patients with cancer.We focus on machine learning approaches,for feature extraction and selection during feature engineering and for imbalanced datasets and multi-modality fusion during statistical modeling.Furthermore,we introduce the stability,reproducibility,and interpretability of features,and the generalizability and interpretability of models.Finally,we offer possible solutions to current challenges in radiomics research.

Discontinuity development patterns and the challenges for 3D discrete fracture network modeling on complicated exposed rock surfaces

Natural slopes usually display complicated exposed rock surfaces that are characterized by complex and substantial terrain undulation and ubiquitous undesirable phenomena such as vegetation cover and rockfalls.This study presents a systematic outcrop research of fracture pattern variations in a complicated rock slope,and the qualitative and quantitative study of the complex phenomena impact on threedimensional(3D)discrete fracture network(DFN)modeling.As the studies of the outcrop fracture pattern have been so far focused on local variations,thus,we put forward a statistical analysis of global variations.The entire outcrop is partitioned into several subzones,and the subzone-scale variability of fracture geometric properties is analyzed(including the orientation,the density,and the trace length).The results reveal significant variations in fracture characteristics(such as the concentrative degree,the average orientation,the density,and the trace length)among different subzones.Moreover,the density of fracture sets,which is approximately parallel to the slope surface,exhibits a notably higher value compared to other fracture sets across all subzones.To improve the accuracy of the DFN modeling,the effects of three common phenomena resulting from vegetation and rockfalls are qualitatively analyzed and the corresponding quantitative data processing solutions are proposed.Subsequently,the 3D fracture geometric parameters are determined for different areas of the high-steep rock slope in terms of the subzone dimensions.The results show significant variations in the same set of 3D fracture parameters across different regions with density differing by up to tenfold and mean trace length exhibiting differences of 3e4 times.The study results present precise geological structural information,improve modeling accuracy,and provide practical solutions for addressing complex outcrop issues.

Hybrid modeling for carbon monoxide gas-phase catalytic coupling to synthesize dimethyl oxalate process

Ethylene glycol(EG)plays a pivotal role as a primary raw material in the polyester industry,and the syngas-to-EG route has become a significant technical route in production.The carbon monoxide(CO)gas-phase catalytic coupling to synthesize dimethyl oxalate(DMO)is a crucial process in the syngas-to-EG route,whereby the composition of the reactor outlet exerts influence on the ultimate quality of the EG product and the energy consumption during the subsequent separation process.However,measuring product quality in real time or establishing accurate dynamic mechanism models is challenging.To effectively model the DMO synthesis process,this study proposes a hybrid modeling strategy that integrates process mechanisms and data-driven approaches.The CO gas-phase catalytic coupling mechanism model is developed based on intrinsic kinetics and material balance,while a long short-term memory(LSTM)neural network is employed to predict the macroscopic reaction rate by leveraging temporal relationships derived from archived measurements.The proposed model is trained semi-supervised to accommodate limited-label data scenarios,leveraging historical data.By integrating these predictions with the mechanism model,the hybrid modeling approach provides reliable and interpretable forecasts of mass fractions.Empirical investigations unequivocally validate the superiority of the proposed hybrid modeling approach over conventional data-driven models(DDMs)and other hybrid modeling techniques.

Exploring Capillary Fringe Flow:Quasilinear Modeling with Kirchhoff Transforms and Gardner Model

Recent studies have underscored the significance of the capillary fringe in hydrological and biochemical processes.Moreover,its role in shallow waters is expected to be considerable.Traditionally,the study of groundwater flow has centered on unsaturated-saturated zones,often overlooking the impact of the capillary fringe.In this study,we introduce a steady-state two-dimensional model that integrates the capillary fringe into a 2-D numerical solution.Our novel approach employs the potential form of the Richards equation,facilitating the determination of boundaries,pressures,and velocities across different ground surface zones.We utilized a two-dimensional Freefem++finite element model to compute the stationary solution.The validation of the model was conducted using experimental data.We employed the OFAT(One_Factor-At-Time)method to identify the most sensitive soil parameters and understand how changes in these parameters may affect the behavior and water dynamics of the capillary fringe.The results emphasize the role of hydraulic conductivity as a key parameter influencing capillary fringe shape and dynamics.Velocity values within the capillary fringe suggest the prevalence of horizontal flow.By variation of the water table level and the incoming flow q0,we have shown the correlation between water table elevation and the upper limit of the capillary fringe.

Effects of Andrade and Burgers rheologies on glacial isostatic adjustment modeling in Antarctica

Variations in ice mass deform the Earth and modify its gravity field,a process known as Glacial Isostatic Adjustment(GIA).GIA in Antarctica remains poorly constrained due to the cumulative effect of past and present ice-mass changes,the unknown history of the past ice-mass change,and the uncertainties on the mechanical properties of the Earth.This paper investigates the effect of using Andrade and Burgers viscoelastic rheologies,rather than the commonly used Maxwell rheology,to model GIA-induced deformation in Antarctica.The Love number and Green's function formalism are used to compute the radial surface displacements and the gravity changes induced by the past and present ice-mass changes.We consider an Earth model whose elastic properties and radial structure are averaged from the Preliminary Reference Earth Model and two viscosity profiles to account for the recently published results on the present ice-mass changes.Using the three rheological laws affects the temporal response of the Earth differently,leading to smaller discrepancies than those induced by the two viscosity structures.The differences are the largest between Maxwell and Burgers rheologies during the 100-1000 years following the beginning of the surface-mass change.Results show that using the Andrade and Burgers rheologies allows the Earth to respond on decennial to centennial time scales,up to 10 m more than Maxwell.Considering only the recent ice-mass changes,the deformation rates derived from Burgers and Andrade rheologies are several times larger than those estimated by Maxwell rheology.

Integrated Hydrological Modeling of the Godavari River Basin in Maharashtra Using the SWAT Model: Streamflow Simulation and Analysis

Hydrological modeling plays a crucial role in efficiently managing water resources and understanding the hydrologic behavior of watersheds. This study aims to simulate daily streamflow in the Godavari River Basin in Maharashtra using the Soil and Water Assessment Tool (SWAT). SWAT is a process-based hydrological model used to predict water balance components, sediment levels, and nutrient contamination. In this research, we used integrated remote sensing and GIS data, including Digital Elevation Models (DEM), land use and land cover (LULC) maps, soil maps, and observed precipitation and temperature data, as input for developing the SWAT model to assess surface runoff in this large river basin. The Godavari River Basin under study was divided into 25 sub-basins, comprising 151 hydrological response units categorized by unique land cover, soil, and slope characteristics using the SWAT model. The model was calibrated and validated against observed runoff data for two time periods: 2003-2006 and 2007-2010 respectively. Model performance was assessed using the Nash-Sutcliffe efficiency (NSE) and the coefficient of determination (R2). The results show the effectiveness of the SWAT2012 model, with R2 value of 0.84 during calibration and 0.86 during validation. NSE values also ranged from 0.84 during calibration to 0.85 during validation. These findings enhance our understanding of surface runoff dynamics in the Godavari River Basin under study and highlight the suit-ability of the SWAT model for this region.

Temperature-dependent solubility of Rebaudioside A in methanol/ethanol and ethyl acetate mixtures:Experimental measurements and thermodynamic modeling

The equilibrium solubility of Rebaudioside A(Reb A)FormⅡin binary mixtures of methanol/ethanol and ethyl acetate was quantitatively determined within the temperature range of 283.15—328.15 K at ambient pressure.The experimental findings indicate a positive correlation between the solubility of Reb A(FormⅡ)and both the temperature and the methanol/ethanol content in the solvent system.To describe the solubility data,six distinct models were employed:the modified Apelblat equation,theλh model,the combined nearly ideal binary solvent/Redlich—Kister(CNIBS/R—K)model,the van't HoffJouyban-Acree(VJA)model,the Apelblat-Jouyban-Acree(AJA)model,and the non-random two-liquid(NRTL)model.The combined nearly ideal binary solvent/Redlich—Kister model exhibited the most precise fit for solubility in methanol+ethyl acetate mixtures,reflected by an average relative deviation(ARD)of 0.0011 and a root mean square deviation(RMSD)of 12×10^(-7).Conversely,for ethanol+ethyl acetate mixtures,the modified Apelblat equation provided a superior correlation(ARD=0.0014,RMSD=4×10^(-7)).Furthermore,thermodynamic parameters associated with the dissolution of Reb A(FormⅡ),including enthalpy,entropy,and the Gibbs energy change,were inferred from the data.The findings underscore that the dissolution process is predominantly endothermic across the solvent systems examined.Notably,the entropy changes appear to have a significant influence on the Gibbs free energy associated with the dissolution of Reb A(FormⅡ),suggesting that entropic factors may play a pivotal role in the studied systems.

Modeling and Validation of Diamagnetic Rotor Levitated by Permanent Magnetics

As an innovative,low-power consuming,and low-stiffness suspension approach,the diamagnetic levitation technique has attracted considerable interest because of its potential applicability in miniaturized mechanical systems.The foundation of a diamagnetic levitation system is mathematical modeling,which is essential for operating performance optimization and stability prediction.However,few studies on systematic mathematical modeling have been reported.In this study,a systematic mathematical model for a disc-shaped diamagnetically levitated rotor on a permanent magnet array is proposed.Based on the proposed model,the magnetic field distribution characteristics,diamagnetic levitation force characteristics(i.e.,levitation height and stiffness),and optimized theoretical conditions for realizing stable levitation are determined.Experiments are conducted to verify the feasibility of the proposed mathematical model.Theoretical predictions and experimental results indicate that increasing the levitation height enlarges the stable region.Moreover,with a further increase in the rotor radius,the stable regions of the rotor gradually diminish and even vanish.Thus,when the levitation height is fixed,a moderate rotor radius permits stable levitation.This study proposes a mathematical modeling method for a diamagnetic levitation system that has potential applications in miniaturized mechanical systems.

Litterfall production modeling based on climatic variables and nutrient return from stands of Eucalyptus grandis Hill ex Maiden and Pinus taeda L.

Native grasslands in the Pampas of South America are increasingly being replaced by Eucalyptus and Pinus stands.The short rotation regimes used for the stands require high nutrient levels,with litterfall being a major source of nutrient return.To model the litterfall production using climatic variables and assess the nutrient return in 14-year-old Eucalyptus grandis and Pinus taeda stands,we measured litter production over 2 years,using conical litter traps,and monitored climatic variables.Mean temperature,accumulated precipitation,and mean maximum vapor pres-sure deficit at the seasonal level influenced litterfall produc-tion by E.grandis;seasonal accumulated precipitation and mean maximum temperature affected litterfall by P.taeda.The regression tree modeling based on these climatic vari-ables had great accuracy and predictive power for E.grandis(N=33;MAE(mean absolute error)=0.65;RMSE(root mean square error)=0.91;R^(2)=0.71)and P.taeda(N=108;MAE=1.50;RMSE=1.59;R^(2)=0.72).The nutrient return followed a similar pattern to litterfall deposition,as well as the order of importance of macronutrients(E.grandis:Ca>N>K>Mg>P;P.taeda:N>Ca>K>Mg>P)and micronutrients(E.grandis and P.taeda:Mn>Fe>Zn>Cu)in both species.This study constitutes a first approximation of factors that affect litterfall and nutrient return in these systems.

WT-FCTGN:A wavelet-enhanced fully connected time-gated neural network for complex noisy traffic flow modeling

Accurate forecasting of traffic flow provides a powerful traffic decision-making basis for an intelligent transportation system. However, the traffic data's complexity and fluctuation, as well as the noise produced during collecting information and summarizing original data of traffic flow, cause large errors in the traffic flow forecasting results. This article suggests a solution to the above mentioned issues and proposes a fully connected time-gated neural network based on wavelet reconstruction(WT-FCTGN). To eliminate the potential noise and strengthen the potential traffic trend in the data, we adopt the methods of wavelet reconstruction and periodic data introduction to preprocess the data. The model introduces fully connected time-series blocks to model all the information including time sequence information and fluctuation information in the flow of traffic, and establishes the time gate block to comprehend the periodic characteristics of the flow of traffic and predict its flow. The performance of the WT-FCTGN model is validated on the public Pe MS data set. The experimental results show that the WT-FCTGN model has higher accuracy, and its mean absolute error(MAE), mean absolute percentage error(MAPE) and root mean square error(RMSE) are obviously lower than those of the other algorithms. The robust experimental results prove that the WT-FCTGN model has good anti-noise ability.

Displacement field reconstruction in landslide physical modeling by using a terrain laser scanner e Part 2:Application and large strain/displacement and water effect analysis

Deformation analysis is fundamental in geotechnical modeling.Nevertheless,there is still a lack of an effective method to obtain the deformation field under various experimental conditions.In this study,we introduce a processebased physical modeling of a pileereinforced reservoir landslide and present an improved deformation analysis involving large strains and water effects.We collect multieperiod point clouds using a terrain laser scanner and reconstruct its deformation field through a point cloud processing workflow.The results show that this method can accurately describe the landslide surface deformation at any time and area by both scalar and vector fields.The deformation fields in different profiles of the physical model and different stages of the evolutionary process provide adequate and detailed landslide information.We analyze the large strain upstream of the pile caused by the pile installation and the consequent violent deformation during the evolutionary process.Furthermore,our method effectively overcomes the challenges of identifying targets commonly encountered in geotechnical modeling where water effects are considered and targets are polluted,which facilitates the deformation analysis at the wading area in a reservoir landslide.Eventually,combining subsurface deformation as well as numerical modeling,we comprehensively analyze the kinematics and failure mechanisms of this complicated object involving landslides and pile foundations as well as water effects.This method is of great significance for any geotechnical modeling concerning large-strain analysis and water effects.

Combining field data and modeling to better understand maize growth response to phosphorus(P) fertilizer application and soil P dynamics in calcareous soils

We used field experimental data to evaluate the ability of the agricultural production system model (APSIM) to simulate soil P availability,maize biomass and grain yield in response to P fertilizer applications on a fluvo-aquic soil in the North China Plain.Crop and soil data from a 2-year experiment with three P fertilizer application rates(0,75 and 300 kg P_(2)O_(5) ha^(–1)) were used to calibrate the model.Sensitivity analysis was carried out to investigate the influence of APSIM SoilP parameters on the simulated P availability in soil and maize growth.Crop and soil P parameters were then derived by matching or relating the simulation results to observed crop biomass,yield,P uptake and Olsen-P in soil.The re-parameterized model was further validated against 2 years of independent data at the same sites.The re-parameterized model enabled good simulation of the maize leaf area index (LAI),biomass,grain yield,P uptake,and grain P content in response to different levels of P additions against both the calibration and validation datasets.Our results showed that APSIM needs to be re-parameterized for simulation of maize LAI dynamics through modification of leaf size curve and a reduction in the rate of leaf senescence for modern staygreen maize cultivars in China.The P concentration limits (maximum and minimum P concentrations in organs)at different stages also need to be adjusted.Our results further showed a curvilinear relationship between the measured Olsen-P concentration and simulated labile P content,which could facilitate the initialization of APSIM P pools in the NCP with Olsen-P measurements in future studies.It remains difficult to parameterize the APSIM SoilP module due to the conceptual nature of the pools and simplified conceptualization of key P transformation processes.A fundamental understanding still needs to be developed for modelling and predicting the fate of applied P fertilizers in soils with contrasting physical and chemical characteristics.

A Novel Framework to Construct S-Box Quantum Circuits Using System Modeling: Application to 4-Bit S-Boxes

Quantum computers accelerate many algorithms based on the superposition principle of quantum mechanics.The Grover algorithm provides significant performance to malicious users attacking symmetric key systems.Since the performance of attacks using quantum computers depends on the efficiency of the quantum circuit of the encryption algorithms,research research on the implementation of quantum circuits is essential.This paper presents a new framework to construct quantum circuits of substitution boxes(S-boxes)using system modeling.We model the quantum circuits of S-boxes using two layers:Toffoli and linear layers.We generate vector spaces based on the values of qubits used in the linear layers and apply them to find quantum circuits.The framework finds the circuit bymatching elements of vector spaces generated fromthe input and output of a given S-box,using the forward search or themeet-in-the-middle strategy.We developed a tool to apply this framework to 4-bit S-boxes.While the 4-bit S-box quantum circuit construction tool LIGHTER-R only finds circuits that can be implemented with four qubits,the proposed tool achieves the circuits with five qubits.The proposed tool can find quantum circuits of 4-bit odd permutations based on the controlled NOT,NOT,and Toffoli gates,whereas LIGHTER-R is unable to perform this task in the same environment.We expect this technique to become a critical step toward optimizing S-box quantum circuits.