Next Article in Journal
An Evolutionary Game-Based Regulatory Path for Algorithmic Price Discrimination in E-Commerce Platforms
Previous Article in Journal
The Discovery of Truncated M-Fractional Exact Solitons and a Qualitative Analysis of the Generalized Bretherton Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

FMEA-TSTM-NNGA: A Novel Optimization Framework Integrating Failure Mode and Effect Analysis, the Taguchi Method, a Neural Network, and a Genetic Algorithm for Improving the Resistance in Dynamic Random Access Memory Components

by
Chia-Ming Lin
and
Shang-Liang Chen
*
Institute of Manufacturing Information and Systems, National Cheng Kung University, Tainan 70101, Taiwan
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(17), 2773; https://doi.org/10.3390/math12172773
Submission received: 13 August 2024 / Revised: 1 September 2024 / Accepted: 5 September 2024 / Published: 7 September 2024

Abstract

:
Dynamic random access memory (DRAM) serves as a critical component in medical equipment. Given the exacting standards demanded by medical equipment products, manufacturers face pressure to improve their product quality. The electrical characteristics of these products are based on the resistance value of the DRAM components. Hence, the purpose of this study is to optimize the resistance value of DRAM components in medical equipment. We proposed a novel FMEA-TSTM-NNGA framework that integrates failure mode and effect analysis (FMEA), the two-stage Taguchi method (TSTM), neural networks (NN), and genetic algorithms (GA) to optimize the manufacturing process. Moreover, the proposed FMEA-TSTM-NNGA framework achieved a substantial reduction in experimental trials, cutting the required number by a factor of 85.3 when compared to the grid search method. Our framework successfully identified optimal manufacturing condition settings for the resistance values of DRAM components: Depo time = 27 s, Depo O2 flow = 151 sccm, ARC-LTO etch time = 43 s, ARC-LTO etch pressure = 97 mTorr, Ox-SiCO etch time = 91 s, Ox-SiCO gas ratio = 22%, and Polish time = 84 s. The results helped the case company improve the resistance value of DRAM components from 191.1 × 10−3 Ohm to 176.84 × 10−3 Ohm, which is closer to the target value of 176.5 × 10−3 Ohm. The proposed FMEA-TSTM-NNGA framework is designed to operate efficiently on resource-constrained, facilitating real-time adjustments to production attributes. This capability enables DRAM manufacturers to swiftly optimize product quality.

1. Introduction

The overall usage of DRAM has been greatly improved, particularly in applications like telemedicine, as the health sector increasingly relies on digital solutions. This trend is expected to continue, with the pandemic accelerating the intensified use of technology in healthcare settings, thereby raising the demand for DRAM components in medical applications [1]. As the demand for high-performing memory tools capable of handling large volumes of data and performing tasks quickly and reliably grows, medical devices are increasingly utilizing DRAM components. As medical technologies advance, DRAM plays a significant role in optimizing devices used for diagnosis, sensing the presence of illnesses, and various other healthcare purposes [2].
Healthcare equipment and appliances require enhanced precision and reliability, and the same is true for the DRAM components used in these devices. This underscores the need for DRAM to meet the stringent demands of complex medical applications that require high performance and data authenticity. The evolution of medical technology requirements indicates a growing need for DRAM with greater stability and performance characteristics [3,4].
These devices are critical instruments in diagnosis and in the monitoring of various physiological parameters, which are crucial for clinical assessment. For example, information from an accurate analysis of olfactory evoked potentials may greatly aid in understanding neurological functions and diagnosing related disorders [5]. Similarly, advancements in deep learning algorithms, such as ResUNet for liver and tumor segmentation in CT images, enhance the level of detail in medical imaging, thereby improving the diagnosis of medical conditions [6]. These technologies ensure that accurate and efficient medical equipment relay information that facilitates informed decision making by healthcare practitioners in patient management.
Regarding a sustainable smart healthcare system, regulations focus on efficacy, dependability, and risk assessment for smart healthcare technology. These stringent requirements arise because healthcare applications are typically mission-critical, where errors can be costly. It is crucial to ensure that all medical devices and systems meet these standards to safeguard patient well-being and deliver dependable healthcare solutions. The COVID-19 pandemic has amplified the need for these criteria, rendering durable and reliable health IT essential to the system [7].
DRAM components are embedded in various types of medical analysis and application devices, including ventilators, medical tablets, medical computers, and physiological monitoring systems. These components provide the memory storage and speed needed to handle substantial amounts of critical data processed in real time, enabling the smooth and effective functioning of medical devices. Integrating DRAM into these systems is crucial for meeting the performance and stability requirements of medical applications [8,9].
Therefore, the difference between the theoretical resistance of DRAM components and the intended value should be as small as possible to ensure optimal results. Significant resistance deviations in DRAM components can lead to disruptions in electronic signals and a reduction in the quality and yield of medical products. This study is set in a semiconductor firm in the Hsinchu Science Park in Taiwan, which specializes in producing DRAM components. During production, the firm encountered issues with electronic signals in DRAM parts, and a further examination revealed that the resistance fluctuations in the DRAMs with signal problems were beyond tolerance. To facilitate timely shipping and meet customer demand, this study will assist the company in analyzing and enhancing these aspects.
During production, if the yield rate is too low, production engineers need to constantly monitor the production situation. Low production yield rates can affect the shipping schedule, so customers often request rapid improvements in production yield rates. To improve product yield rates, production engineers often use previous production data for predictive improvement. Although there is a lot of production data, most of it falls within the specification tolerance range, which can only be predicted within that region, not beyond it. Because of these limitations, the results of yield rate improvements are often not very significant. In addition, when production engineers conduct product improvement experiments, they often use trial-and-error methods, which consume a large amount of experimental resources. Moreover, these trial-and-error experiments often result in imbalanced data, which increases the prediction error.
In order to improve the resistance values of DRAM components in medical equipment, it is necessary to first collect a high-quality dataset for predicting the resistance values of DRAM components. Therefore, there is a need to introduce a fast and effective data collection method. Before collecting data, it is necessary to identify the important production attributes on the production line. Thus, this study will use the FMEA and Taguchi method to identify the important production attributes with the minimum number of experiments and help the case company obtain the most valuable production data. Then, the Taguchi experimental data will be integrated with artificial neural network for modeling, in order to find the best attribute setting of production attribute to assist the case study in improving the resistance values of DRAM components.

2. Motivation for Problem Description

The resistance values of DRAM components in medical equipment can vary due to different production attributes settings. As the manufacturing process becomes more complex, improving the resistance value of DRAM components requires considering multiple possible production attributes. During production, production engineers often rely on past engineering experience to make improvements, which can be time-consuming and delay product shipping schedules. Companies often find that their equipment schedules are overloaded, customers require tight improvement times, and experiment resources are limited, among other constraints. Therefore, in order to meet shipping schedule as quickly as possible, the best production attribute level setting must be determined within a short period of time.
The semiconductor manufacturing processes affecting the resistance of DRAM components include deposition, photolithography, etching, and polishing [10,11,12,13] and are described below:
(1)
Thin film deposition operates by building up structures on a wafer using different types of materials in a layered manner. The first deposition method is chemical vapor deposition (CVD), which involves the use of corrosive gasses and chemicals which provide a thin covering to the wafers upon a chemical reaction. First of all, the gas containing target material is injected into the reaction chamber, and the wafer is located in this chamber. After that, by adjusting parameters such as the temperature and pressure, the gas molecules are prompted to chemically react with the aim of depositing the correct thin film material. These materials are deposited at the surface of the wafer, and after some time, they stack up and develop a film-like structure. The second deposition technique is physical vapor deposition (PVD), which is an advanced technology that transfers physical state materials into a gas or plasma state and deposits them on the surface of the wafer. Some PVD technologies are sputter deposition and evaporation deposition. Often, the material that is going to be deposited in any PVD technique becomes heated or is bombed first so that it can evaporate in the form of gas or plasma. These high-energy particles are then directed to the surface of the wafer, and in the process of collision, it is transformed into a thin film structure.
(2)
A microfabrication technique that is closely related to the manufacturing of a semiconductor is photolithography, where fine semiconductor circuit patterns are printed. The four sub-process operations are cleaning, photoresist application, exposure, and development, which the wafers undergo in the main process. Before photolithography is carried out, there is always a need to cleanse the wafer to remove any particles which may hinder the photolithography process. The cleaning methods which are employed for the wafers are solution washing, high-pressure water cleaning, and ultrasonic cleaning. Then, a layer of material known as the photoresist, a polymer that can be dissolved by exposure to light of a specific wavelength, is formed on the wafer’s surface. The photoresist is typically used and can be dispensed uniformly on the surface of the wafer by the spin-coating method. Incidentally, the photomask may be exposed to a light source of a given wavelength with the intention of projecting some patterns on the surface of the wafer. These light rays will only pass through the parts of the photo res dose where the thickness of the photo res dose will be consistent with that of the underlying layer, and they will not pass through areas where the layers of the photo res dose will have to be stripped off. Finally, the wafer is put into a developer solution where the other parts of the photoresist layer exposed to light are dissolved. This method is called positive photoresist development, where the exposed layer of the photoresist becomes the pattern in the shape and size of the required design.
(3)
In the etching process, the un-masked areas in the wafer are chemically or physically etched to make the circuit structure of the required semiconductor materials by erasing the unwanted areas. Dry etching can be described as the technique of taking material on a wafer’s surface by means of gas plasma or reactive plasma. In the course of the dry etching process, the wafer is located in the vacuum chamber, and a gas flow containing the reactive gas is used. Finally, the gas experiences a reaction under plasma excitation, which allows the removal of material deposited on the surface of the wafer. Wet etching also involves the process of removing material on the wafer’s surface by an actual chemical process which occurs in a liquid medium. In the wet etching process, the wafer is completely submerged in a bath which contains an etching solution. The etching solution decomposes chemically with the materials on the wafer, and this results in their removal.
(4)
The polishing process is employed in the manufacturing process to planarize the wafer. The steps include coarse polishing, fine polishing, cleaning, and inspection. First, the wafer is placed on the rotating polishing disk, and then, with the polishing particles, a preliminary polish is performed to minimize surface roughness, protrusions, and depressions on the wafer. However, an additional step of planarization is required to further smooth the wafer’s surface. This step uses the smallest polishing particles and the highest pressure to achieve the desired flatness. During the polishing process of ion-implanted silicon wafers, polish debris and dirt are extensively produced. Therefore, it is necessary to perform chemical or physical cleaning on the wafer’s surface to meet the required cleanliness and smoothness conditions, as the presence of high amounts of polishing powder and dirt may hinder this. After polishing and cleaning, there may still be slight roughness or scratches on the wafer’s surface in terms of flatness. To further enhance the flatness and surface quality of the wafer, other methods like chemical mechanical polishing or even conventional mechanical polishing can be employed. Lastly, the wafers are required to undergo a final cleaning step to ensure that the entire surface is free of any dirt and polishing residue and to verify that the wafer’s flatness meets the required standards.

3. Methodology

3.1. Taguchi Method

The Taguchi method (TM) is an experimental design method which can be used to improve the product quality of manufacturing processes [14,15]. One of the key characteristics of the Taguchi method is the use of orthogonal arrays (OA). OA can make the experiment more efficient [16]. Table 1 shows the Taguchi orthogonal array L12(211). The letters A~K in the OA represent different production attributes, while the numbers 1 and 2 represent the two different levels for each production attribute. Each column in OA represents a specific production attribute. L12(211) has a total of 11 columns, which means that up to 11 production attributes can be arranged. The number of columns in OA corresponds to the number of experiments, so the number of experiments in the L12(211) orthogonal array is 12.
Genichi Taguchi developed SN (signal-to-noise) to assess the impact of manufacturing process attributes on product quality. SN is a measure of the quality of a signal relative to the background noise. In manufacturing, SN can be used to identify key attributes that significantly affect product quality. A higher SN indicates greater process stability and higher product quality. Conversely, a lower SN indicates process instability, which may require corrective measures. The SN can be categorized into three types based on quality requirements [17]: Nominal The Better (NTB), Smaller The Better (STB), and Larger The Better (LTB), each with a corresponding Formula (1)–(3). The SN has become an essential tool in quality control and improvement, widely applied by industries worldwide.
In this context, y ¯ denotes the average value of a group of experiments; while m represents the target quality value; S2 represents the variance of each group of experiments. It is calculated as the sum of the squared deviations from the mean, divided by the degrees of freedom; The value of the ith experiment is denoted by yi, and n represents the total number of experiments conducted.
S N N T B = 10 log i = 1 n ( y i m ) 2 n = 10 log S 2 y ¯ 2
S N S T B = 10 log i = 1 n y i 2 n = 10 log ( y ¯ 2 + S 2 )
S N L T B = 10 log i = 1 n 1 y i 2 n
The algorithmic process of the Taguchi method is as follows:
(1)
Problem Definition and Objective Setting: The output attributes and the inputs related to the performance, for instance, the quality of the product, are determined. When writing down the specifications, it is necessary to state the target performance issues of the system as well as define its target values or response curves if any.
(2)
Experimental Design: The control factors that define the system’s performance are chosen and their possible levels are defined. As noted in the Taguchi method, it is possible to test more than one input factor. It is recommended to apply Taguchi’s orthogonal array to build an experimental matrix that would cover all possible combinations of the factor level settings.
(3)
Experimentation and Data Collection: Experiments are performed on the basis of designed specifications of the input and output characteristics. For each run of the experiment, the response values (the output characteristics) are documented.
(4)
Data Analysis: Regarding quality characteristics, the SN is used in order to determine the stability of the system. On that basis, the key factors that play a crucial role in shaping the output characteristics are defined, and their importance is assessed.
(5)
Optimization and Validation: Multiple comparison is performed using the post hoc test to define the best level of the factor on the basis of the mean value of the SN ratio and the ANOVA. Using a practical test or more precise factor confirmation experiments, the advantages of the selected factor setting are confirmed.
The pseudo code of the Taguchi method is determined as follows:
(1)
Define the Problem
To ensure clarity on what the experiment would set out to achieve, the following is carried out:
The output (response) characteristics, which should be optimized, are identified.
The relations between the input factors (variables) and output are determined.
(2)
Select Factors and Levels
The research objectives and questions or hypotheses are formulated for the control factors for the study.
The method used to measure each factor is described.
(3)
Design the Experiment
An appropriate orthogonal array is chosen depending on the factor number and level number.
The experimental matrix is developed through the choice of orthogonal array, chosen above.
(4)
Conduct the Experiments
The experiments are carried out as per the design matrix laid out above.
A record of the final output values obtained at the end of each of the experimental runs is made.
(5)
Analyze the Data
The signal-to-noise (SN) ratio values of each experimental runs are determined.
The trends in the SN ratios could be analyzed in order to obtain an idea about the performance of the system.
An analysis of variance (ANOVA) test is performed to determine the contribution of each factor.
(6)
Optimize the Process
The combination of factor levels is determined, which could be considered as optimal according to the analysis.
All of the process parameters are changed to the best possible values.

3.2. Artificial Neural Network

Artificial neural networks (ANNs) are among the most popular tools used for data analysis. This approach involves the creation of a model through the use of statistics from a given dataset of input attributes and output variables. The model identifies the interdependence of variables and can be used for prediction, decision making, and diagnosing engineering problems, as pointed out in [18].
In fact, to train a neural network, the weights of the links between neurons must be changed iteratively. An important fact is that weights determine the probabilities of activation in connected neurons and affect the output variable. Neural networks are formed by integrating several neurons; this can be explained by the following reasons [19]. In Figure 1 below, we provide a three-layer neural network architecture. The training process in the network has to be performed multiple times so that the input attributes are mapped correctly to the right output variable. The primary goal of training is to minimize error and ensure the output value of the neural network is as close as possible to the engineering target. This means that training is considered complete once the error rate no longer fluctuates, thus suggesting that the network has completely converged. It then requires evaluation on a set of samples that were never used in training to ensure that the predicted value matches the output variable in the sample set [20]. With respect to training neural networks, the learning rate is generally a hyperparameter that determines the speed at which the network converges. The advantage of a higher learning rate is that it results in convergence at a much faster rate than a lower learning rate [21,22].
The pseudo code of an artificial neural network is determined as follows:
(1)
Initialize the Network
The number of neurons in the input layer, hidden layer, and output layer is decided.
Initial values of the weights and biases are assigned with respect to each neuron (these may preferably be assigned randomly).
(2)
Forward Propagation
The following is carried out for each input sample:
  • Input Layer:
    Let the input layer take in the input data.
  • Hidden Layer:
    The following is carried out for each hidden layer neuron:
    The weighted sum of inputs is found, and bias is added to them.
    An activation function, such as ReLU, sigmoid, or other, is applied to the sum.
  • Output Layer:
    The following is carried out for each output neuron:
    The sum of the products of the weight in the connection between two neurons and the hidden layer’s output plus the bias term are found.
    An activation function is applied to the sum to come up with the final output.
(3)
Calculate Loss
The difference between the predicted output and the target output is found by using a loss function, e.g., mean squared error or cross entropy.
(4)
Backpropagation
The backpropagation algorithm is used to determine the gradient with respect to all of the weights and biases of the network.
The following is carried out for each layer:
  • The derivative of the activation function is determined.
  • The first-order derivative of the loss with the weights and the biases is calculated.
  • The weights and the biases are updated according to the gradients calculated at the beginning of the current iteration.
(5)
Iterate
Steps 2–4 are repeated for a certain number of epochs or until the loss function falls within an acceptable range.
(6)
Output the Results
Once training is complete, the network can be used to predict the outputs of other inputs.
The effectiveness of the network is assessed using a testing dataset.

3.3. Genetic Algorithm

Holland J. developed the genetic algorithm (GA) in 1975. Since then, the GA has undergone extensive development. The GA has been shown to be an effective search algorithm, especially in the proposed model, as it can locate the global optimal solution rather than just a local one. Through natural selection, the GA identifies the fittest set of individuals to survive, forming an adaptive optimization technology capable of searching the high-dimensional solution space.
The GA is widely used in work related to hyperparameter tuning in machine learning models and engineering parameter optimization [23]. The evolutionary process of the GA consists of three main steps: chromosome replication or duplication, crossover, and mutation. The GA operates in the solution space of engineering problems, and after multiple iterations, it finds the most adaptable solution that satisfies all constraint conditions. This is the chromosome with the highest fitness function value, which is the global optimal solution [24]. The algorithmic process of the GA is as follows:
(1)
Initialize population: First, a set of potential solutions is randomly created, which forms the first generation of the population.
(2)
Evaluate population: The performance value of each possible solution is measured, and it is referred to as the fitness function value.
(3)
Select parents: Some solutions with high fitness function values are chosen as parents. In most cases, roulette wheel selection or the elite preservation strategy is employed for the selection of parents.
(4)
Crossover: Two parents are randomly chosen to crossover in order to create two new chromosomes.
(5)
Mutation: Sometimes, it is necessary to mutate the offspring with some certain probability in order to avoid falling into local optima.
(6)
Select next generation: Numbers of offspring and parents are chosen as the next generation of the population.
(7)
Iteration steps 2 to 6 are performed over and over until either the condition for convergence is reached or the maximum number of iterations has been exhausted.
(8)
Finally, the candidate with the highest objective value is chosen as the global optimal solution to the presented mathematical model.
The pseudo code of the genetic algorithm is determined as follows:
(1)
Initialize Population
The first generation of the candidate solutions, commonly known as the chromosomes, is produced.
Each chromosome is a candidate solution, and most commonly, these chromosomes are stored as binary strings of numbers and real numbers, amongst other formats.
(2)
Evaluate Fitness
The following is carried out for each chromosome in the population:
  • The chromosome is translated into the solution on the chromosome it represents.
  • Another measure of the solution obtained is determined by considering the goal or objective function of the problem.
(3)
Selection
Parent chromosomes in the current population are chosen using the fitness measure defined.
Higher fitness enhances the likelihood of selection, for example, by means of roulette wheel and tournament selections.
(4)
Crossover
The following is carried out for each pair of selected parents:
  • The crossover point(s) is selected as a victim.
  • Equal segments of the parent’s chromosomes are exchanged at crossover point(s) to form the offspring commonly referred to as children.
With a certain crossover probability, a crossover operation should be performed in order to generate new offsprings.
(5)
Mutation
The following is carried out for each offspring:
  • A possibility of mutation is created that will cause some arbitrary change to some genes in the chromosome.
  • It must be ensured that mutation contributes to the creation of diversity but does not interfere with searching.
(6)
Create New Population
The waiting list state is populated by replacing the old population with the new population, which is the offspring.
Optionally, some of the chromosomes with good quality are transferred to the old population in order to guarantee the existence of good solutions.
(7)
Iterate
Steps 2–6 are repeated for a specified number of generations or until a termination condition is met (e.g., the fitness threshold or maximum number of generations).
(8)
Output the Best Solution
After the final generation, the chromosome with the highest fitness is identified and output as the best solution found by the algorithm.

4. The Proposed Methodology

To improve the quality of DRAM components in medical equipment within a limited timeframe, this study employed a novel FMEA-TSTM-NNGA framework, which integrates Failure mode and effect analysis (FMEA), two stage Taguchi method (TSTM), artificial neural network (NN), and genetic algorithm (GA). The primary objective is to aid the case company in enhancing product yield despite resource and time limitations. Using these methods, the research aims to improve the resistance value of DRAM components, thereby achieving a higher yield rate.
To ensure high-quality data, Taguchi method can be a beneficial approach. Unlike trial-and-error or one-factor-at-a-time experiments, Taguchi method allows for the simultaneous consideration of multiple production attributes on the resistance value of DRAM components under different experimental level settings. This Taguchi method also helps balance data, using the characteristics of orthogonal arrays. Moreover, Taguchi method can aid in identifying important production attributes that affect the resistance of DRAM components. Traditionally, companies have relied on work experience to select important production attributes for experiments. However, this approach can result in the testing of insignificant effects, wasting valuable experimental resources. By determining the important attributes of the product before initiating quality improvement efforts and then optimizing and adjusting them, companies can more efficiently achieve the best quality performance of the product.
Taguchi method is a useful for finding an optimal production attribute setting under a specific experimental setting level. However, Taguchi method may only identify a local optimal solution rather than the best global production attribute setting for the entire experimental region. To address this limitation, artificial intelligence algorithms such as artificial neural network and genetic algorithms can be used to further optimize product quality.
Artificial neural networks are particularly well-suited for performing nonlinear modeling on the complex relationship between input production attributes and the resistance of DRAM components. The predictive model can then be used to help identify the best production attribute setting across the entire experiment region. By incorporating these AI techniques into the design and optimization process, companies can achieve even better product quality while minimizing the amount of time and resources spent on experimentation.
The proposed methodology of this study can be divided into seven steps, as shown in Figure 2. A schematic diagram of the methodology is provided in Figure 3.
  • Defining the problem: In this step, the existing abnormal phenomena and their yield loss rate of defects in the case study are described. The manufacturing process that may be causing the abnormalities is identified based on the adverse phenomena observed.
  • Using FMEA to assess production attributes: FMEA was used to assess the production attributes affecting memory component resistance in thin film deposition, etching, plating, and polish processes. Relevant production attribute setting is then identified to improve the average resistance of DRAM components in the medical equipment through the optimization of these production attributes.
  • Applying first stage Taguchi method to find important attributes: In this step, we will use the Taguchi method L12(211) orthogonal array to collect experimental data. We will then use attribute response tables, attribute response charts, and ANOVA (analysis of variance) tables to analyze the data and identify important production attributes that affect the average resistance of DRAM components. At present, many researchers have used Taguchi method to assist in the optimization and improvement of engineering problems. Gao and Zhou (2024) employs the Taguchi Method to conduct sensitivity analysis on operating parameters for a proton exchange membrane fuel cell (PEMFC) combined heat and power (CHP) cogeneration system in variable climate regions, proposing cost-effective strategies for electricity- and heat-dominated outputs [25]. Xie et al. (2024) used the Taguchi method to explore key parameters in pure waterjet surface treatment of Ti6Al4V specimens, emphasizing the significance of operation pressure for biomaterial surface fine-tuning [26]. Tanürün et al. (2024) employed the Taguchi method to optimize Vertical Axis Wind Turbine (VAWT) performance with an adaptive flap design, achieving a 74.01% increase in power coefficient (CP) compared to conventional VAWTs, with flap position identified as the most influential factor [27]. Lu et al. (2024) used the Taguchi method to optimize CO2 laser polishing of fused silica, revealing laser beam scanning speed as a crucial factor, resulting in a substantial reduction in surface roughness from Ra = 0.157 μm to 0.005 μm [28]. Yang et al. (2024) utilized the Taguchi method and Gray relation analysis to identify the optimal design parameters (perforation length, height, and tilt angle) for slit fins with lateral perforations, achieving the highest j-factor and lowest f-factor [29]. Based on the empirical evidence presented in the literature, it is clear that Taguchi’s experimental design can effectively assist researchers in identifying important parameters across various fields and finding optimal parameter settings for different engineering problems.
  • Applying second stage Taguchi method to recollect data and find out the optimal attribute setting: To collect experimental data, we used the Taguchi method L18(21 × 37) orthogonal array. We conducted a second Taguchi method in the production attribute region, with the aim of recollecting experimental data that could potentially contain optimal solutions using the L18(21 × 37) orthogonal array. We analyzed the collected data using attribute response tables, attribute response charts, and ANOVA table to determine the optimal attribute setting in the region.
  • Using artificial neural networks to build a model: To construct the relationship between the important production attributes and the resistance value of DRAM components of medical equipment in the L18(21 × 37) orthogonal array, an artificial neural network can be used. Recent studies frequently use neural networks to establish relationships between production attributes and response variables, and subsequently use these models to determine the optimal attribute settings and enhance product quality. Acharjee et al. (2024) developed two artificial neural network based models (W-ANN and H-ANN) for predicting dynamic modulus in Colombian hot-mix asphalt mixtures, outperforming previous models and providing practical tools for pavement design with reduced testing requirements [30]. Liu (2024) presents a sensor array based on noble metal-doped In2O3, employing a back propagation neural network (BPNN) integrated with the whale optimization algorithm (WOA) for anti-interference detection of mixed NOX, achieving quantitative prediction of components in the presence of cross interference [31]. Mudawar et al. (2024) employed artificial neural networks to predict heat transfer and critical heat flux (CHF) in both microgravity and Earth gravity [32]. Sayed et al. (2024) applied artificial intelligence (ANN and FBI) to optimize yeast and wastewater concentrations in a microalgae microbial fuel cell, achieving significant improvements in power density and COD removal [33]. Sanni et al. (2024) employed an adaptive neuro-fuzzy inference system (ANFIS) to model the corrosion rate of AISI 316 stainless steel under various inhibitor dosages and schedules, achieving superior prediction accuracy [34].
  • Adopting genetic algorithm to identify the best production attribute setting: The genetic algorithm will be used to identify the best setting of production attribute levels across the entire production attribute region. Many researchers have also confirmed the feasibility of the genetic algorithm. Park et al. (2024) integrated genetic algorithms and deep learning to enhance the I-V modeling of DRAM transistors, addressing limitations of the BSIM model and accurately modeling devices with hot-carrier degradation effects [35]. Mukhanov et al. (2020) introduced DStress, a framework using genetic algorithms to identify worst-case DRAM error patterns, enhancing testing by detecting critical reliability issues and improving testing mechanisms [36]. Leu et al. (2021) used genetic algorithms to reduce thickness deviation in DRAM manufacturing. The approach successfully decreased deviation from 45.0 to 12.9 Å [37]. Babu et al. (2024) proposed a novel IoT network optimization strategy using a genetic algorithm (GA) for energy efficiency and mixed integer linear programming (MILP) for node deployment, with blockchain enhancing data privacy; outperforming existing models in network lifetime and throughput [38]. Nigam et al. (2024) accelerated molecular design using a genetic algorithm and an artificial neural network, identifying over 10,000 potential organic emitters with inverted singlet-triplet gaps (INVEST) and appreciable fluorescence rates for potential use in new-generation organic light-emitting diodes [39]. Wang et al. (2024) applied an improved genetic algorithm for optimal scheduling in hybrid energy ship power systems, reducing costs and environmental impacts [40]. Bhat et al. (2024) enhanced aluminum alloy design using machine learning, incorporating data-driven classes and genetic algorithms for feature optimization, achieving improved tensile strength and elongation predictions [41]. Ildarabadi et al. (2024) applied a genetic algorithm with elitism mechanism to optimize the cost of restoring power to a distribution network after disasters by building Tie-lines between damaged and healthy sections [42]. GA is an optimization technique that can be applied to a wide range of engineering problem areas. Therefore, it is favored by researchers who want to find optimal solutions to complex problems.
  • Performing confirmation experiments and comparisons: After identifying the optimal settings for production attributes, confirmation experiments are carried out to confirm the efficacy of the suggested methodology., and the results are compared with the average resistance of DRAM components before improvement to confirm the effectiveness of the improvements.
The pseudo code of the proposed methodology is determined as follows:
(1)
Define the problem.
The various abnormal phenomena and the yield loss rate of defects that is associated with it are explained.
The manufacturing processes that is responsible for the formation of such abnormalities is determined.
(2)
Examine the possibility of failure in the future and identify product attributes in the FMEA.
The FMEA is used to assess the production characteristics that impact memory component resistance.
The critical production attributes of thin film deposition, etching and plating, and polishing are determined.
The optimal settings of the production attributes that will cause the DRAM components to have an average resistance greater than past records are determined.
(3)
Use the first stage of the Taguchi method.
The L12(211) orthogonal array is chosen to plan the experiments.
Research is conducted and information is obtained.
The analysis of attribute data can be carried out with the help of attribute response tables, response chars, and an ANOVA in order to determine important production attributes influencing DRAM component resistance.
(4)
Use the second stage of the Taguchi method.
Further experimentation should be carried out using the L18(21 × 37) orthogonal array.
More tests should be conducted in order to gather information on the area of important attributes.
This means that one has to perform the analysis of the obtained data to identify the appropriate values of production attributes.
(5)
Develop a model employing artificial neural networks.
An artificial neural network is performed with the dependent variable of resistance values of DRAM components, with the independent variables being the important production attributes.
The neural network model should be trained and validated using the experimental data collected from the second stage of the Taguchi method.
(6)
Optimize the production attribute settings with the genetic algorithm.
A genetic algorithm search is performed regarding the overall attribute region of the production attribute in order to obtain the global optima of the attribute settings.
The suitability of several contexts is analyzed and the most appropriate one is determined.
(7)
Confirmation experiments and comparison.
Confirmation experiments are carried out using the best production attribute settings which were found.
The findings are assessed with pre-change data to ensure that the above proposed methodology has a powerful impact.
Details, especially the outline of the results of their work, are determined with an emphasis on increases in the average resistance of DRAM parts.

5. Case Study

Details Content of the Proposed Methodology

There are seven steps involved in the proposed methodology:
(1)
Problem Definition
This case study is based on a DRAM manufacturing company located in Hsinchu Science Park, Taiwan. In 2023, it was discovered that abnormal resistance values occurred in the DRAM manufacturing processes. By the second half of 2023, the number of abnormal wafers had reached 360, and the yield loss rate had increased to 1.89%, significantly impacting the DRAM quality. To address the abnormality of the DRAM components, the company analyzed the production data from 2023 and created a trend chart of weekly average resistance values, as shown in Figure 4. The chart indicates a gradual increase in the resistance values of the DRAM components, with the abnormal position occurring at RM5_01 in a memory component, as shown in Figure 5. The primary reason for the abnormal resistance value was not optimized production attributes in the deposition, etching, plating, and polishing processes. This study analyzes the production data of DRAM components from 2023 and finds that the average resistance value is 191.1 × 10−3 Ohm, which exceeds the target value of 176.5 × 10−3 Ohm specified by the customer. These findings are shown in Figure 6. Therefore, the purpose of this study is to optimize the production attributes of these DRAM components to improve their average resistance value and meet the customer’s specification target in a short period.
(2)
Using FMEA to assess production attributes
The semiconductor manufacturing processes that impact the resistance of DRAM components in medical equipment primarily include thin film deposition, etching, plating, and polishing processes. These processes have specific names, namely SiCO TF depo, ARC-LTO etch, Oxide-SiCO etch, Cu ECP, and Cu polish, as shown in Figure 7. The abbreviations used in Figure 7 are SiCO for Silicon-OxyCarbide, ARC for Anti-Reflective Coating, LTO for Low Temperature Oxide, and ECP for Electro-Chemical Plating.
To assess the production attributes that may affect the resistance of DRAM components in the thin film deposition, etching, plating, and polish processes, the case company formed a group of production engineers to improve the relevant production attributes for these processes. The production engineers evaluate three indicators, Severity (S), Occurrence (O), and Detection (D), based on the criteria of failure mode and effect analysis (FMEA). Severity (S) presents the magnitude of the impact of a failure. Severity is assessed on a scale of 1 to 10, where 1 indicates a minor impact and 10 indicates a very severe impact. Occurrence (O) shows the probability of a failure occurring. The probability is evaluated on a scale of 1 to 10, where 1 indicates a very low probability and 10 indicates a very high probability. Detection (D) presents how easily a failure can be detected when it occurs. Detection is assessed on a scale of 1 to 10, where 1 indicates detection is easy and 10 means detection is very difficult.
The corresponding S, O, and D values are assessed for each production attribute, and the risk priority number (RPN) is calculated as the process of S, O, and D. A higher RPN value indicates a higher risk, requiring the implementation of more control measures to mitigate the risk. The results of this evaluation are shown in Table 2. Based on the team’s professional experience, they selected ten production attributes that may affect the DRAM components in medical equipment. These attributes include Depo time (s), Depo He flow (sccm), Depo O2 flow (sccm), ARC-LTO etch time (s), ARC-LTO etch pressure (mTorr), Ox-SiCO etch time (s), Ox-SiCO gas ratio (%), Ox-SiCO CHF3 flow (sccm), Polish time (s), and Polish pressure (psi).
(3)
Applying first stage Taguchi method to find important attributes
In the first stage of the Taguchi method, ten production attributes potentially affecting DRAM components in medical equipment were considered, and the levels of these attributes were determined based on the engineer’s experience with orthogonal arrays. The experimental setting levels for the production attributes are shown in Table 3, where Level 1 shows a low setting level and Level 2 shows a high setting level. Using the Taguchi method, the L12(211) orthogonal array was selected as the screening experiment to identify the important production attributes. The production attributes (A~J) are, respectively, arranged in columns 1~10 of the L12(211) orthogonal array. In this phase, if we use a full factorial design to identify better setting for ten production attributes, each possessing two levels, would demand a staggering 1024 experiments (210 = 1024). However, employing the L12(211) orthogonal array, this investigation efficiently conducted a mere 12 experiments, encompassing diverse production attribute settings.
Manufacturing engineers were assigned to conduct a Taguchi method for the identified production attributes. In this design, a total of 12 experiments were carried out with different setting levels for each production attribute. Table 4 shows the measured resistance values of five DRAM components for each experiment of the Taguchi orthogonal array L12(211). Table 5 shows the outcomes of the attribute response table regarding the average resistance value of DRAM components, while Figure 8 shows the results of the attribute response chart.
According to Table 5, the impact of each production attribute on the average resistance of the memory component is ranked as: F(8.8) > E(7.5) > I(4.2) > A(3.6) > H(2.6) > J(2.1) > C(2.0) > G(1.5) > B(1.0) > D(0.8). Based on Figure 8, the better settings of production attributes can be set as A = 29 s, E = 85 mTorr, F = 80 s, and I = 90 s for average resistance.
Furthermore, as the focus of this study is on improving the resistance value of DRAM components, which is a nominal-the-better (NTB) quality characteristic, the signal-to-noise ratio (SN) of each experiment was calculated using Formula (1) and the results are presented in the last column of Table 4. Table 6 and Figure 9 exhibit the attribute response table and attribute response chart for the SN of DRAM components, respectively.
According to the information presented in Table 6, the production attributes’ effects on the SN of DRAM components can be ranked in descending order: E(2.3) > C(2.2) > G(1.7) > I(1.4) > D(1.4) > A(0.8) > J(0.4) > H(0.3) > F(0.2) > B(0.1). Based on these findings, the better setting of production properties can be determined as A = 33 s, C = 160 sccm, D = 45 s, E = 55 mTorr, G = 30%, and I = 90 s for SN, as shown in Figure 9.
To further investigate the optimal production attribute settings for improving the resistance value of the DRAM components, a new Taguchi method L18(21 × 37) will be conducted. The study will use the analysis results of the average resistance value and SN of the DRAM components from the previous experiment, focusing on seven important production attributes: A (Depo time), C (Depo O2 flow), D (ARC-LTO etch time), E (ARC-LTO etch pressure), F (Ox-SiCO etch time), G (Ox-SiCO gas ratio), and I (Polish time). The new experiment aims to recollect relevant data within the region of experiments that may include the optimal solution and determine the optimal setting of production attribute levels suggested by Taguchi method.
(4)
Applying second stage Taguchi method to recollect data and find out the optimal attribute setting
In this step, this study adopts the Taguchi method once again to gather relevant data for the experimental region that is likely to contain the optimal solution. Based on the average resistance value of the DRAM components obtained in the previous step, the better attribute setting for average resistance is A = 29 s, E = 85 mTorr, F = 80 s, and I = 90 s. On the other hand, the best attribute setting for SN is A = 33 s, C = 160 sccm, D = 45 s, E = 55 mTorr, G = 30%, and I = 90 s. However, since there is a conflict between the levels of production attributes A and E, the production engineers decided to focus on the better attribute setting of the average resistance and finalized the attribute combination as A = 29 s, C = 160 sccm, D = 45 s, E = 85 mTorr, F = 80 s, G = 30%, and I = 90 s.
In the second stage, we adjusted the levels of the selected attributes based on the better attribute setting, as shown in Table 7. For example, attribute A (Depo time) was initially set to 29 and 33 s in the first stage. However, the experiments revealed that A = 29 s yielded a better average resistance value. Consequently, in the second stage of the Taguchi method, we refined the attribute levels for A, adjusting them to 25, 27, and 29 s to further optimize the settings and continue identifying the best settings.
To rearrange new experimental levels, the experimental design at this step will use the L18(21 × 37) orthogonal array, where the seven important production attributes (A, C, D, E, F, G, and I) will be arranged in the L18(21 × 37) columns 2 to 8 of the orthogonal array, as shown in Table 8. During this phase, employing a full factorial design to pinpoint optimal setting for seven production attributes, each featuring three levels, would necessitate a substantial 2187 experiments (37 = 2187). However, by using the L18(21 × 37) orthogonal array, this study conducted a streamlined set of only 18 experiments, effectively exploring a wide range of production attribute settings.
Table 8 shows the experimental results for the resistance value of DRAM components and the SN value of Taguchi’s experimental design L18(21 × 37) at different experimental level settings. Table 9 and Figure 10 present the attribute response table and attribute response chart, respectively, for the average resistance analysis results of DRAM components. The effect of each production attribute on the average resistance value of DRAM components in descending order is shown in Table 9 as follows: F(7.57) > A(6.73) > I(6.56) > E(5.68) > G(0.87) > D(0.54) > C(0.37). As seen in Figure 10, the optimal settings of production attributes are set as A = 25 s, E = 95 mTorr, F = 100 s, and I = 85 s for average resistance.
Table 10 and Figure 11 are the attribute response table and attribute response chart, respectively, for the SN analysis results of DRAM components. From Table 10, it can be inferred that the production attribute D has the highest contribution to the SN of DRAM components, followed by attributes A, G, C, I, E, and F in descending order of importance. In Figure 11, the optimal settings of production attributes can be set as A = 25 s, C = 150 sccm, D = 35 s, and G = 30% for SN.
Based on these analysis results of the average resistance and SN, the optimal settings of attributes for the Taguchi method are A = 25 s, C = 150 sccm, D = 35 s, E = 95 mTorr, F = 100 s, G = 30%, and I = 85 s.
To confirm the optimal production attribute setting (A = 25 s, C = 150 sccm, D = 35 s, E = 95 mTorr, F = 100 s, G = 30%, and I = 85 s) determined by the Taguchi method, the case company conducted three batches on the production line. Each batch sampled one wafer, and the resistance of five DRAM components from each wafer was measured to confirm the experimental data. Table 11 shows the measurement results. Initially, the average resistance value of the memory component was 191.1 × 10−3 Ohm, and the deviation from the design target value of 176.5 × 10−3 Ohm was 14.6 × 10−3 Ohm. However, after applying the Taguchi method, the average resistance value was reduced to 172.1 × 10−3 Ohm, and the deviation from the design target value was only 4.4 × 10−3 Ohm, which represents an improvement of 69.89%. Table 12 compares the results before and after the implementation of the Taguchi method, confirming that the local optimal setting of production attribute levels identified by the Taguchi method is suitable for mass production. Nevertheless, to further reduce the deviation from the design target value of 176.5 × 10−3 Ohm, this study plans to use an artificial neural network to construct the relationship between the production attributes and the resistance value of DRAM components in the next stage.
(5)
Using artificial neural networks to build a model
Based on results of ANOVA, we can identify the important production attributes as A (Depo time), C (Depo O2 flow), D (ARC-LTO etch time), E (ARC-LTO etch pressure), F (Ox-SiCO etch time), G (Ox-SiCO gas ratio), and I (Polish time). These attributes will be used as input variables for the artificial neural network, while the resistance value of the memory component will serve as the output variable.
To create an artificial neural network model with a rich and varied dataset, this research utilized the Monte Carlo technique to enhance the resistance data sourced from Table 8. The Monte Carlo method introduces randomness or stochastic components into its forecasts, offering varying results with each simulation. By employing a probabilistic model like Monte Carlo, diverse datasets can be generated, leading to a broader scope of outcomes [43]. In this investigation, a normal distribution probabilistic model was adopted for augmenting the data using the Monte Carlo approach. The outcomes of this data augmentation process are presented in Table 13, specifically in rows P1 to P10.
To train the network, 80% of the dataset is randomly selected as the training dataset, and the remaining 20% is used as the testing dataset. Finally, the optimal neural network architecture was identified as 7-2-6-5-1, as shown in Figure 12.
(6)
Adopting genetic algorithm to identify the best production attribute setting
In this step, the objective function for the genetic algorithm is defined as the inverse of the deviation between the predicted resistance values from the previously constructed 7-2-6-5-1 neural network architecture and the target design values.
Finally, the genetic algorithm finds that the minimum resistance deviation of the memory component is 1.2 × 10−4 Ohm, and the global optimal production attributes setting are set as follows: A (Depo time) = 27 s, C (Depo O2 flow) = 151 sccm, D (ARC-LTO etch time) = 43 s, E (ARC-LTO etch pressure) = 97 mTorr, F (Ox-SiCO etch time) = 91 s, G (Ox-SiCO gas ratio) = 22%, and I (Polish time) = 84 s.
(7)
Performing confirmation experiments and comparisons
To verify the proposed framework FMEA-TSTM-NNGA in this study, a small trial run was conducted by the case company on their production line to identify the optimal production attributes for improving the resistance value of DRAM components. The trial run involved producing 3 wafers, and 5 DRAM components were selected from each wafer. The results of the trial run are shown in Table 14.
Using artificial neural networks and genetic algorithms, the DRAM components of the medical equipment produced by the case company were improved. As a result, there was a significant improvement in the resistance value of the DRAM components. Specifically, the average resistance deviation of the memory element was reduced by 97.67%, from 14.6 × 10−3 Ohm to 0.34 × 10−3 Ohm, as shown in Table 15.
Based on the results of the trial run and the optimization process, it can be concluded that the best setting of the production attributes for DRAM components is a mass-production setting.
Additionally, the execution results of the proposed methodology can save more experimental resources compared to other existing optimization methods. The first stage of the Taguchi method employs the L12(211) orthogonal array, significantly minimizing the number of experimental trials needed to identify critical production attributes that influence memory component resistance. Typically, a full factorial experiment would require 1024 tests to cover all combinations of ten production attributes across two levels. However, using the L12(211) array, the number of necessary experiments is reduced to just 12, resulting in an 85.3-fold decrease in the experimental workload. The second stage of the Taguchi method applies the L18(21 × 37) orthogonal array. Normally, a complete factorial design would necessitate 2187 experiments to explore all combinations of seven production attributes at three levels. The L18(21 × 37) array reduces the number of required experiments to 18, thereby achieving a 121.5-fold reduction in experimental effort compared to the full factorial experiment. In the genetic algorithm (GA) phase, 20 initial population tests were conducted, followed by 300 iterations, resulting in a total of 6000 runs. In comparison, a grid search approach would require approximately 2.4 × 10⁷ runs (calculated as 4 × 20 × 10 × 15 × 20 × 10 × 10), demonstrating a substantial reduction in the number of model evaluations.

6. Numerical Simulations and Discussion

(a)
Analysis of variance (ANOVA)
When researchers aim to analyze multiple categorical independent variables and compare the means across different groups, an ANOVA becomes essential [44]. Specifically, when the independent variables are categorical and the dependent variable is continuous, a statistical analysis is necessary to evaluate the relationships between the groups’ means. This analysis helps determine how variations in the dependent variable are affected by the different levels of the independent variables.
For the first stage of the Taguchi method, Table 16 provides the analysis of variance (ANOVA) of the average resistance values of DRAM components for the orthogonal array L12(211). Based on Table 16, it can be observed that the p-values for production attributes A, E, F, and I are below 0.1, indicating that these four production attributes have a significant contribution to the average resistance of DRAM components and are important production attributes that affect resistance. Table 17 provides the ANOVA results of the SN of DRAM components for the orthogonal array L12(211). Table 17 shows that the p-values for production attributes A, C, D, E, G, and I are less than 0.1. Therefore, these six production attributes have a significant contribution to the SN of DRAM components and are important production attributes that affect SN.
Regarding the second stage of the Taguchi method, Table 18 presents the ANOVA results of the average resistance analysis of DRAM components for the Taguchi’s experimental design L18(21 × 37). Based on the information provided in Table 18, it can be observed that the p-values for production attributes A, E, F, and I are below 0.1, indicating that these four production attributes significantly contribute to the average resistance of DRAM components and are important production attributes that affect the resistance of DRAM components. Table 19 shows the ANOVA results of the SN analysis of DRAM components for the orthogonal array L18(21 × 37). Table 19 indicates that the p-values for production attributes A, C, D, and G are less than 0.1, signifying that these attributes significantly contribute to the SN of DRAM components and are important factors that affect SN.
(b)
Neural networks (NN)
This study proposes an artificial neural network architecture with three hidden layers. To train the network, 80% of the dataset (Table 13) is randomly selected as the training dataset, and the remaining 20% is used as the testing dataset. The neural network model is configured with the ReLU activation function and the Adam optimizer. A grid search method is used to find the optimal settings for hyperparameters, including the learning rate, momentum rate, and number of nodes in each hidden layer. The learning rate is tested with values of [0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4], while the momentum rate is tested with values of [0.7, 0.75, 0.8, 0.85, 0.9, 0.95]. The number of nodes in each hidden layer is set as [1:10]. The initial conditions of the neural network are shown in Figure 13.
The artificial neural network is trained on 1000 iterations, and the minimum root mean square error (RMSE) of the testing dataset is used as the criterion to select the final weights. After conducting the grid search, the optimal hyperparameters settings are determined as follows: the learning rate is set to 0.1, the momentum rate is set to 0.95, and the number of nodes in the three hidden layers are set to 2, 6, and 5, respectively. Therefore, the optimal neural network architecture is 7-2-6-5-1.
To verify the convergence of the optimal hyperparameters identified by the grid search, we compare the RMSE of the training and testing datasets for the ANN model with the best hyperparameters. Figure 14 shows the convergence curve, which plots RMSE values against the number of iterations, offering essential insights into a model’s training dynamics. In the early learning phase, both the training and testing RMSE values are usually elevated, reflecting the model’s initial state with weights that are far from optimal. As training progresses, these RMSE values typically decline, signaling enhanced model performance. Around the 60th iteration, the RMSE values for both the training and testing datasets converge closely, indicating consistent learning. However, if the gap between the training and testing RMSE values begins to widen, it could indicate overfitting, where the model becomes overly adapted to the training data, thereby reducing its generalization capability. In the final stages of training, the RMSE values tend to stabilize. If the testing RMSE values stabilize near the training RMSE values, this strongly suggests that the model has achieved robustness.
(c)
Genetic algorithm (GA)
The genetic algorithm (GA) operates by first standardizing the input variables of artificial neural networks. Once standardized, the values are converted into binary code, which is then combined into a string and used in subsequent genetic processes, including reproduction, crossover, and mutation. In this study, the fitness function is defined as the reciprocal of the deviation between the memory component and the target value, as shown in Formula (4):
F i t n e s s   f u n c t i o n = 1 R e s i s t a n c e T a r g e t   v a l u e
In this study, the initial operation of the genetic algorithm begins with the random generation of 20 populations for reproduction, Using the roulette method. The roulette wheel’s area size is determined by the fitness function value, where a larger area on the wheel indicates a higher fitness function value and easier reproduction of the next generation. Based on past experience, the crossover rate and mutation rate for the genetic algorithm have been established as 0.85 and 0.05, respectively. The stopping condition of the genetic algorithm is set to 300 iterations.
Table 20 shows the production attribute setting levels and the resistance value deviations of the DRAM components after each genetic algorithm. The genetic algorithm finds that the minimum resistance deviation of the memory component is 1.2 × 10−4 Ohm, and the global optimal production attribute settings are set as follows: A (Depo time) = 27 s, C (Depo O2 flow) = 151 sccm, D (ARC-LTO etch time) = 43 s, E (ARC-LTO etch pressure) = 97 mTorr, F (Ox-SiCO etch time) = 91 s, G (Ox-SiCO gas ratio) = 22%, and I (Polish time) = 84 s.

7. Conclusions

In this case study, we used FMEA to identify ten production attributes in thin film deposition, etching, plating, and polish processes. We then conducted experiments using the Taguchi method L12(211) with a limited number of trials to identify key production attributes that affect memory component resistance. Using the factor response table and an ANOVA, we identified seven critical attributes: A (Depo time), C (Depo O2 flow), D (ARC-LTO etch time), E (ARC-LTO etch pressure), F (Ox-SiCO etch time), G (Ox-SiCO gas ratio), and I (Polish time). The dataset, derived from seven key attributes, was collected using the Taguchi method L18(21 × 37). This dataset was subsequently utilized as both the training and testing dataset for an artificial neural network in the following stage.
Next, we integrated the Monte Carlo method with a normal distribution probabilistic model into the L18(21 × 37) to expand the modeling dataset. During the artificial neural network modeling process, the artificial neural network in this study was fed with seven important attributes as input variables. The optimal hyperparameter setting of the artificial neural network were then found through the grid search method, with a learning rate of 0.1, momentum rate of 0.95, and three hidden layers with 2, 6 and 5 nodes, respectively. Finally, the resulting architecture was 7-2-6-5-1.
Subsequently, this study applied a genetic algorithm to determine the best production attribute setting across the entire region, resulting in the following settings: A (Depo time) = 27 s, C (Depo O2 flow) = 151 sccm, D (ARC-LTO etch time) = 43 s, E (ARC-LTO etch pressure) = 97 mTorr, F (Ox-SiCO etch time) = 91 s, G (Ox-SiCO gas ratio) = 22%, and I (Polish time) = 84 s.
After verifying the experimental results, this study found that the average resistance of the DRAM components had improved from 191.1 × 10−3 Ohm to 176.84 × 10−3 Ohm, which is very close to the target value of 176.5 × 10−3 Ohm. The average resistance deviation was also reduced from 14.6 × 10−3 Ohm to 0.34 × 10−3 Ohm, representing an improvement of 97.67%. The FMEA-TSTM-NNGA framework presented in this paper enables comprehensive analysis of the process and helps identify high-risk factors within it, ensuring that these key considerations are accounted for during improvement efforts.
This framework demonstrates potential for scalability and adaptability across various semiconductor manufacturing processes. However, it may require customization, including adjustments to the input and response variables, to align with the specific needs of different processes. Additionally, the training data for the ANN and the GA’s objective functions may need to be reconfigured and calibrated to match new process objectives, ensuring that the optimization outcomes remain accurate and applicable.
The key limitation of this research is the reliance on a single dataset, which may constrain the generalizability of the findings. To assess the broader applicability of the proposed FMEA-TSTM-NNGA framework, future research could investigate its effectiveness in other semiconductor fields and validate its generalizability using independent resistance datasets from multiple companies.
Our proposed framework can be seamlessly adapted to address image-based semiconductor product defect detection tasks using deep learning techniques. Notably, its design prioritizes efficient training and accurate inference, making it particularly well suited for deployment on resource-constrained edge AI devices within semiconductor manufacturing environments.

Author Contributions

S.-L.C. initiated the study. S.-L.C. and C.-M.L. designed the experiments and algorithm. C.-M.L. wrote the code for the algorithm and performed the experiments. S.-L.C. and C.-M.L. analyzed the experimental results. C.-M.L. and S.-L.C. wrote and edited the manuscript. All authors contributed to the preparation of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mathews, K. COVID-19’s Impact on the Semiconductor Industry. 2021. Available online: https://kidder.com/trend-articles/covid-19s-impact-on-the-semiconductor-industry/ (accessed on 30 August 2024).
  2. Kim, K. The smallest engine transforming humanity: The past, present, and future. In Proceedings of the 2021 IEEE International Electron Devices Meeting (IEDM), San Francisco, CA, USA, 11–16 December 2021. [Google Scholar] [CrossRef]
  3. Cohen, S.; Hershcovitch, M.; Taraz, M.; Kißig, O.; Wood, A.; Waddington, D.; Chin, P.; Friedrich, T. Drug repurposing using link prediction on knowledge graphs with applications to non-volatile memory. In Complex Networks & Their Applications X: Volume 2, Proceedings of the Tenth International Conference on Complex Networks and Their Applications COMPLEX NETWORKS 2021; Springer International Publishing: Cham, Switzerland, 2022; Volume 1073, pp. 742–753. [Google Scholar] [CrossRef]
  4. Shiratake, S. Scaling and performance challenges of future DRAM. In Proceedings of the 2020 IEEE International Memory Workshop (IMW), Dresden, Germany, 17–20 May 2020. [Google Scholar]
  5. Arpaia, P.; Cataldo, A.; Criscuolo, S.; De Benedetto, E.; Masciullo, A.; Schiavoni, R. Assessment and Scientific Progresses in the Analysis of Olfactory Evoked Potentials. Bioengineering 2022, 9, 252. [Google Scholar] [CrossRef]
  6. Rahman, H.; Bukht, T.F.N.; Imran, A.; Tariq, J.; Tu, S.; Alzahrani, A. A Deep Learning Approach for Liver and Tumor Segmentation in CT Images Using ResUNet. Bioengineering 2022, 9, 368. [Google Scholar] [CrossRef] [PubMed]
  7. Chen, T.C.T. Sustainable smart healthcare applications: Lessons learned from the COVID-19 pandemic. In Sustainable Smart Healthcare: Lessons Learned from the COVID-19 Pandemic; Springer Nature: Cham, Switzerland, 2023; pp. 65–92. [Google Scholar] [CrossRef]
  8. Kim, D.H.; Kwon, H.J.; Bae, S.J. DRAM circuit and process technology. In Semiconductor Memories and Systems; Woodhead Publishing: Sawston, UK, 2022; pp. 87–117. [Google Scholar] [CrossRef]
  9. Kim, S.S.; Yong, S.K.; Kim, W.; Kang, S.; Park, H.W.; Yoon, K.J.; Sheen, D.S.; Lee, S.; Hwang, C.S. Review of semiconductor flash memory devices for material and process issues. Adv. Mater. 2023, 35, 2200659. [Google Scholar] [CrossRef]
  10. Linliu, K. DRAM (Dynamic Random Access Memory) Process Flow; Independently Published: 2020. Available online: https://www.amazon.com/Dynamic-Random-Access-Memory-Process/dp/B084DLVSSP (accessed on 12 August 2024).
  11. Xiao, H. 3D IC Devices, Technologies, and Manufacturing; SPIE Press Monographs: Washington, WA, USA, 2016. [Google Scholar]
  12. Prince, B. Semiconductor Memories: A Handbook of Design, Manufacture and Application, 2nd ed.; John Wiley & Sons: New York, NY, USA, 1996. [Google Scholar]
  13. Kuo, Y. Thin Film Transistor Technologies VI: Proceedings of the International Symposium; The Electrochemical Society: New Jersey, NJ, USA, 2003. [Google Scholar]
  14. Phadke, M.S. Quality Engineering Using Robust Design; Prentice-Hall: New Jersey, NJ, USA, 1989. [Google Scholar]
  15. Fowlkes, W.Y.; Creveling, C.M. Engineering Methods for Robust Product Design; Addison-Wesley Publishing Company: Boston, MA, USA, 1995. [Google Scholar]
  16. Su, C.T. Quality Engineering: Off-Line Methods and Applications, 1st ed.; CRC Press: Boca Raton, FL, USA, 2013. [Google Scholar]
  17. Ross, P.J. Taguchi Techniques for Quality Engineering: Loss Function, Orthogonal Experiments, Parameter and Tolerance Design, 2nd ed.; McGraw-Hill: New York, NY, USA, 1996. [Google Scholar]
  18. Rosenblatt, F. Principles of Neurodynamics: Perceptrons and the Theory of Brain Mechanisms; Spartan Books: New York, NY, USA, 1962. [Google Scholar]
  19. Stern, H.S. Neural networks in applied statistics. Technometrics 1996, 38, 205–214. [Google Scholar] [CrossRef]
  20. McClelland, J.L.; Rumelhart, D.E. Explorations in Parallel Distributed Processing: A Handbook of Models, Programs, and Exercises, 1st ed.; MIT Press: London, UK, 1989. [Google Scholar]
  21. Fausett, L. Fundamentals of Neural Networks: An Architecture, Algorithms, and Applications; Prentice Hall: New Jersey, NJ, USA, 1994. [Google Scholar]
  22. Hagan, M.T.; Demuth, H.B.; Beale, M. Neural Network Design; PWS: Boston, MA, USA, 1995. [Google Scholar]
  23. Goldberg, D.E. Genetic Algorithm in Search, Optimization and Machine Learning, 1st ed.; Addison-Wesley: New York, NY, USA, 1989. [Google Scholar]
  24. Renders, J.M.; Flasse, S.P. Hybrid methods using genetic algorithms for global optimization. IEEE Trans. Syst. Man Cybern. Part B (Cybern.) 1966, 26, 243–258. [Google Scholar] [CrossRef]
  25. Gao, B.; Zhou, Y. In-field experimental study and multivariable analysis on a PEMFC combined heat and power cogeneration for climate-adaptive buildings with taguchi method. Energy Convers. Manag. 2024, 301, 118003. [Google Scholar] [CrossRef]
  26. Xie, J.; Qiao, Y.; Qi, Y.; Xu, Q.; Shemtov-Yona, K.; Chen, P.; Rittel, D. Application of the Taguchi method to areal roughness-based surface topography control by waterjet treatments. Appl. Surf. Sci. Adv. 2024, 19, 100548. [Google Scholar] [CrossRef]
  27. Tanürün, H.E. Improvement of vertical axis wind turbine performance by using the optimized adaptive flap by the Taguchi method. Energy Sources Part A Recovery Util. Environ. Eff. 2024, 46, 71–90. [Google Scholar] [CrossRef]
  28. Lu, G.; Li, X.; Wang, D.; Wang, K. Process Parameter Optimization for CO2 Laser Polishing of Fused Silica Using the Taguchi Method. Materials 2024, 17, 709. [Google Scholar] [CrossRef]
  29. Yang, P.; Wang, L.; Qing, T.; Liu, Y. Multi-objective optimization of geometrical parameters of laterally perforated on silt fin by gray relation analysis based on Taguchi method. Int. J. Therm. Sci. 2024, 196, 108705. [Google Scholar] [CrossRef]
  30. Acharjee, P.K.; Souliman, M.I.; Freyle, F.; Fuentes, L. Development of dynamic modulus prediction model using artificial neural networks for colombian mixtures. J. Transp. Eng. Part B Pavements 2024, 150, 04023038. [Google Scholar] [CrossRef]
  31. Liu, Y.; Yang, Z.; Huang, L.; Zeng, W.; Zhou, Q. Anti-interference detection of mixed NOX via In2O3-based sensor array combining with neural network model at room temperature. J. Hazard. Mater. 2024, 463, 132857. [Google Scholar] [CrossRef] [PubMed]
  32. Mudawar, I.; Darges, S.J.; Devahdhanush, V.S. Prediction technique for flow boiling heat transfer and critical heat flux in both microgravity and Earth gravity via artificial neural networks (ANNs). Int. J. Heat Mass Transfer. 2024, 220, 124998. [Google Scholar] [CrossRef]
  33. Sayed, E.T.; Rezk, H.; Abdelkareem, M.A.; Olabi, A.G. Artificial neural network based modelling and optimization of microalgae microbial fuel cell. Int. J. Hydrogen Energy 2024, 52, 1015–1025. [Google Scholar] [CrossRef]
  34. Sanni, O.; Adeleke, O.; Ukoba, K.; Ren, J.; Jen, T.C. Prediction of inhibition performance of agro-waste extract in simulated acidizing media via machine learning. Fuel 2024, 356, 129527. [Google Scholar] [CrossRef]
  35. Park, J.H.; Kim, J.N.; Lee, S.; Kim, G.J.; Lee, N.; Baek, R.H.; Kim, D.H.; Kim, C.; Kang, M.; Kim, Y. Current-Voltage Modeling of DRAM Cell Transistor Using Genetic Algorithm and Deep Learning. IEEE Access 2024, 12, 23881–23886. [Google Scholar] [CrossRef]
  36. Mukhanov, L.; Nikolopoulos, D.S.; Karakonstantis, G. Dstress: Automatic synthesis of dram reliability stress viruses using genetic algorithms. In Proceedings of the 2020 53rd Annual IEEE/ACM International Symposium on Microarchitecture (MICRO), Athens, Greece, 17–21 October 2020; pp. 298–312. [Google Scholar] [CrossRef]
  37. Leu, Y.; Lin, C.M.; Yang, W.N. Reducing Thickness Deviation of W-Shaped Structures in Manufacturing DRAM Products Using RSM and ANN_GA. IEEE Trans. Compon. Packag. Manuf. Technol. 2021, 11, 899–910. [Google Scholar] [CrossRef]
  38. Babu, R.M.; Satamraju, K.P.; Gangothri, B.N.; Malarkodi, B.; Suresh, C.V. A hybrid model using genetic algorithm for energy optimization in heterogeneous internet of blockchain things. Telecommun. Radio Eng. 2024, 83, 1–16. [Google Scholar] [CrossRef]
  39. Nigam, A.; Pollice, R.; Friederich, P.; Aspuru-Guzik, A. Artificial design of organic emitters via a genetic algorithm enhanced by a deep neural network. Chem. Sci. 2024, 15, 2618–2639. [Google Scholar] [CrossRef]
  40. Wang, X.; Zhu, H.; Luo, X.; Chang, S.; Guan, X. An energy dispatch optimization for hybrid power ship system based on improved genetic algorithm. Proc. Inst. Mech. Eng. Part A J. Power Energy 2024, 238, 348–361. [Google Scholar] [CrossRef]
  41. Bhat, N.; Barnard, A.S.; Birbilis, N. Inverse Design of Aluminium Alloys Using Genetic Algorithm: A Class-Based Workflow. Metals 2024, 14, 239. [Google Scholar] [CrossRef]
  42. Ildarabadi, R.; Lotfi, H.; Hajiabadi, M.E. Resilience enhancement of distribution grids based on the construction of Tie-lines using a novel genetic algorithm. Energy Syst. 2024, 15, 371–401. [Google Scholar] [CrossRef]
  43. Kroese, D.P.; Brereton, T.; Taimre, T.; Botev, Z.I. Why the Monte Carlo method is so important today. Wiley Interdiscip. Rev. Comput. Stat. 2014, 6, 386–392. [Google Scholar] [CrossRef]
  44. St»Hle, L.; Wold, S. Analysis of variance (ANOVA). Chemom. Intell. Lab. Syst. 1989, 6, 259–272. [Google Scholar] [CrossRef]
Figure 1. Neural network architecture such as three layers.
Figure 1. Neural network architecture such as three layers.
Mathematics 12 02773 g001
Figure 2. Proposed methodology.
Figure 2. Proposed methodology.
Mathematics 12 02773 g002
Figure 3. A schematic diagram for proposed methodology.
Figure 3. A schematic diagram for proposed methodology.
Mathematics 12 02773 g003
Figure 4. Trend chart of weekly average resistance values of DRAM components.
Figure 4. Trend chart of weekly average resistance values of DRAM components.
Mathematics 12 02773 g004
Figure 5. Abnormal position of DRAM components.
Figure 5. Abnormal position of DRAM components.
Mathematics 12 02773 g005
Figure 6. Results of the analysis of resistance values in DRAM components.
Figure 6. Results of the analysis of resistance values in DRAM components.
Mathematics 12 02773 g006
Figure 7. Manufacturing process of DRAM components.
Figure 7. Manufacturing process of DRAM components.
Mathematics 12 02773 g007
Figure 8. Average resistance of memory component for L12(211) in attribute response chart.
Figure 8. Average resistance of memory component for L12(211) in attribute response chart.
Mathematics 12 02773 g008
Figure 9. SN of memory component for L12(211) in attribute response chart.
Figure 9. SN of memory component for L12(211) in attribute response chart.
Mathematics 12 02773 g009
Figure 10. Average resistance of memory component for L18(21 × 37) in attribute response chart.
Figure 10. Average resistance of memory component for L18(21 × 37) in attribute response chart.
Mathematics 12 02773 g010
Figure 11. SN of memory component for L18(21 × 37) in attribute response chart.
Figure 11. SN of memory component for L18(21 × 37) in attribute response chart.
Mathematics 12 02773 g011
Figure 12. Proposed architecture of artificial neural network.
Figure 12. Proposed architecture of artificial neural network.
Mathematics 12 02773 g012
Figure 13. Hyperparameter settings of neural network.
Figure 13. Hyperparameter settings of neural network.
Mathematics 12 02773 g013
Figure 14. Convergence curve of training and testing datasets for the proposed 7-2-6-5-1 neural network.
Figure 14. Convergence curve of training and testing datasets for the proposed 7-2-6-5-1 neural network.
Mathematics 12 02773 g014
Table 1. Orthogonal array L12(211).
Table 1. Orthogonal array L12(211).
NoABCDEFGHIJK
111111111111
211111222222
311222111222
412122122112
512212212121
612221221211
721221122121
821212221112
921122212211
1022211112212
1122121211122
1222112121221
Table 2. Relevant attributes in the production process of DRAM components.
Table 2. Relevant attributes in the production process of DRAM components.
NoProcessAttributesSeverityOccurrenceDetectionRPN
1SiCo TF-depoHtrSpace1/HtrSpace234448
2Pressure (Torr)43336
3RF1Time/RF2Time (s)782112
4HighFreqRF1Pwr/HighFreqRF2Pwr (W)34336
5LowFreqRF1Pwr/LowFreqRF2Pwr (W)54360
6OMCTS Flow Set (mgm)23212
7HE_OMCTS Flow Set (sccm)34448
8O2 Flow Set (sccm)24432
9ARC-LTO etchE.Time (s)34224
10Pressure (mTorr)983216
11RF Power 60 MHz33436
12RF Power 13.56 MHz24432
13HV (-V)32318
14Gas O2 (sccm)53460
15Gas CF4 (sccm)23424
16Gas CHF3 (sccm)42432
17Gas ratio (center%)23424
18Temp Top (°C)44348
19Temp Wall (°C)34336
20Temp Bottom (°C)34336
21Oxide-SiCO etchE.Timt (s)773147
22Pressure (mTorr)35345
23RF Power (W) 60 MHz34448
24RF Power (W) 13.56 MHz33436
25HV (-V)33327
26Gas O2 (sccm)53460
27Gas CF4 (sccm)24432
28Gas CHF3 (sccm)954180
29Gas ratio (center%)844128
30Temp Top (°C)53345
31Temp Wall (°C)23318
32Temp Bottom (°C)33327
33Cu polishHead States: Zone 1–5 (psi)684192
34End Step Criteria: Max Time (s)782112
Table 3. Production attributes and their setting levels of L12(211).
Table 3. Production attributes and their setting levels of L12(211).
FactorDepo Time
(s)
Depo He Flow
(sccm)
Depo O2 Flow
(sccm)
ARC-LTO Etch Time
(s)
ARC-LTO Etch Pressure
(mTorr)
Ox-SiCO Etch Time
(s)
Ox-SiCO Gas Ratio
(%)
Ox-SiCO CHF3 Flow
(sccm)
Polish Time
(s)
Polish Pressure
(psi)
ABCDEFGHIJ
Level 1298501604555603040902.2
Level 23395017055858050601004.2
Table 4. Taguchi method L12(211) and experimental data.
Table 4. Taguchi method L12(211) and experimental data.
EXP.Control FactorsResistance (Ohm × 10−3)Average Resistance (Ohm × 10−3)Standard DeviationSN
ABCDEFGHIJN1N2N3N4N5
11111111111224.44223.55222.37223.08231.02224.8933.50536.15
21111122222219.18227.35218.39229.68225.53224.0275.01433.00
31122211122219.59211.61225.11227.83231.26223.0807.70329.24
41212212211222.06226.81209.83225.08220.27220.8116.64430.43
51221221212205.75217.19213.93217.38206.87212.2245.58431.60
61222122121227.24223.62215.12218.87207.80218.5287.56229.22
72122112212235.70234.86228.51221.27232.19230.5055.87431.88
82121222111214.47215.74204.49210.04201.42209.2326.20430.56
92112221221215.29225.17214.69227.10216.86219.8215.85831.49
102221111221242.10234.50238.00230.50229.20234.8645.32632.89
112212121112219.78222.01227.80222.68228.03224.0613.68035.69
122211212122226.32226.94217.65230.61232.51226.8055.72331.96
222.4045.72332.01
Table 5. Average resistance of memory component for L12(211) in attribute response table.
Table 5. Average resistance of memory component for L12(211) in attribute response table.
FactorABCDEFGHIJ
Level 1220.6221.9223.4222.0226.1226.8223.2221.1220.3221.4
Level 2224.2222.9221.4222.8218.7218.0221.7223.7224.5223.5
Effect3.61.02.00.87.58.81.52.64.22.1
Rank49710218536
Table 6. SN of memory component for L12(211) in attribute response table.
Table 6. SN of memory component for L12(211) in attribute response table.
FactorABCDEFGHIJ
Level 131.6132.0533.1232.6933.1432.0932.8432.1332.7231.79
Level 232.4131.9630.9031.3230.8831.9331.1731.8831.3032.23
Effect0.800.092.221.372.260.161.670.251.420.44
Rank61025193847
Table 7. Important attributes and setting levels of L18(21 × 37).
Table 7. Important attributes and setting levels of L18(21 × 37).
FactorDepo Time
(s)
Depo O2 Flow
(sccm)
ARC-LTO Etch Time
(s)
ARC-LTO Etch Pressure
(mTorr)
Ox-SiCO Etch Time
(s)
Ox-SiCO Gas Ratio
(%)
Polish Time
(s)
ACDEFGI
Level 1251403585802080
Level 2271504095902585
Level 329160451001003090
Table 8. Taguchi method L18(21 × 37) and experimental data.
Table 8. Taguchi method L18(21 × 37) and experimental data.
EXP.ParametersResistance (Ohm × 10−3)Average Resistance (Ohm × 10−3)Standard DeviationSN
ACDEFGIN1N2N3N4N5
11111111178.88179.21176.94175.53175.96177.301.67640.49
21222222172.46170.81176.63172.20172.89173.002.17538.01
31333333166.60165.81171.62168.77175.64169.694.01932.51
42112233184.08177.80179.70180.63180.56180.552.27837.98
52223311166.48168.81173.47166.24173.67169.743.64433.36
62331122179.81190.91182.01185.43188.41185.324.53132.23
73121323182.61179.65186.47183.07186.50183.662.89336.05
83232131187.20181.33178.89178.84181.11181.483.41234.52
93313212174.91184.51178.74179.64178.54179.273.44234.33
101133221163.96167.04168.76167.05172.11167.782.97435.03
111211332176.48175.28175.07174.30174.94175.210.79946.82
121322113178.98186.19184.16177.66189.01183.204.80031.63
132123132178.66175.52181.45183.50179.51179.733.00435.54
142231213179.63184.18189.82188.90178.26184.155.24030.92
152312321169.94171.84168.87173.72170.49170.971.87039.22
163132312185.80176.34176.06174.46174.23177.384.80031.35
173213123185.89183.95189.85183.68182.00185.083.00635.79
183321231186.32174.52177.17180.35180.22179.724.41032.20
177.963.2835.44
Table 9. Average resistance of memory component for L18(21 × 37) in attribute response table.
Table 9. Average resistance of memory component for L18(21 × 37) in attribute response table.
LevelACDEFGI
1174.36177.73178.06180.89182.02178.51174.50
2178.41178.11178.17177.76177.41177.63178.32
3181.10178.03177.63175.21174.44177.73181.06
Effect6.730.370.545.687.570.876.56
Rank2764153
Table 10. SN of memory component for L18(21 × 37) in attribute response table.
Table 10. SN of memory component for L18(21 × 37) in attribute response table.
LevelACDEFGI
137.4236.0739.1136.4535.0333.6835.80
234.8836.5734.4735.4534.7536.0636.38
334.0433.6932.7634.4336.5536.6034.15
Effect3.372.886.352.031.812.912.23
Rank2416735
Table 11. Confirmation results of optimal production attribute setting in the second Taguchi method.
Table 11. Confirmation results of optimal production attribute setting in the second Taguchi method.
EXP.Resistance (Ohm × 10−3)Average Resistance (Ohm × 10−3)Standard DeviationSN
N1N2N3N4N5
Wafer #1174.11170.99170.55172.35173.45172.2891.53141.03
Wafer #2170.22173.00170.88171.26174.75172.0241.84039.42
Wafer #3173.99170.24171.76171.68172.32171.9991.35242.09
172.1041.57440.84
Table 12. Comparison results before and after Taguchi method.
Table 12. Comparison results before and after Taguchi method.
ComparisonDepo Time
(s)
Depo O2 Flow
(sccm)
ARC-LTO Etch Time
(s)
ARC-LTO Etch Pressure
(mTorr)
Ox-SiCO Etch Time
(s)
Ox-SiCO Gas Ratio
(%)
Polish Time
(s)
Average Resistance (Ohm × 10−3)|Resistance-Target|
ACDEFGI
Before improvement301605070703090191.10014.60
TSTM2515035951003085172.1044.40
Improvement 69.89%
Table 13. Results of Monte Carlo method.
Table 13. Results of Monte Carlo method.
EXP.123456789101112131415161718
ParametersA111222333111222333
C123123123123123123
D123123231312231312
E123231123312312231
F123231312231123312
G123312231231312123
I123312312123231231
Resistance
(Ω × 10−3)
N1178.88172.46166.60184.08166.48179.81182.61187.20174.91163.96176.48178.98178.66179.63169.94185.80185.89186.32
N2179.21170.81165.81177.80168.81190.91179.65181.33184.51167.04175.28186.19175.52184.18171.84176.34183.95174.52
N3176.94176.63171.62179.70173.47182.01186.47178.89178.74168.76175.07184.16181.45189.82168.87176.06189.85177.17
N4175.53172.20168.77180.63166.24185.43183.07178.84179.64167.05174.30177.66183.50188.90173.72174.46183.68180.35
N5175.96172.89175.64180.56173.67188.41186.50181.11178.54172.11174.94189.01179.51178.26170.49174.23182.00180.22
Average resistance (Ω × 10−3)177.30173.00169.69180.55169.74185.32183.66181.48179.27167.78175.21183.20179.73184.15170.97177.38185.08179.72
Standard deviation1.6762.1754.0192.2783.6444.5312.8933.4123.4422.9740.7994.8003.0045.2401.8704.8003.0064.410
P1175.49175.64168.98181.54166.08184.77181.34178.92179.58164.53174.14184.37175.26189.12169.55178.84188.84180.67
P2175.57177.75168.86184.14165.04174.95179.59181.78177.59169.11174.11171.38175.46180.28167.70168.27182.54184.95
P3176.41172.49168.01181.16171.37183.02180.45182.13178.67163.65174.91185.76179.00183.43173.11172.70183.76179.24
P4178.07176.62172.49184.93169.12184.13177.83184.68178.62166.33174.30187.18181.42180.79171.52175.42183.60183.12
P5179.28171.66173.72182.50164.08184.81183.09180.23179.59167.09175.34181.80183.47191.01173.38175.46184.04175.49
P6178.62173.16167.40177.90166.32188.50184.34181.04184.53165.41174.92184.90180.94184.55172.23174.54179.67187.32
P7178.31173.14170.93177.77173.12179.95186.30177.83179.78168.78175.19187.13181.15177.96170.00165.00188.14177.46
P8179.34173.43173.18182.69175.49184.86185.62177.25174.56167.96173.43182.45176.91195.20172.24172.46189.22187.29
P9175.73175.99168.53179.65170.09190.29178.04178.35172.94173.19175.73183.19183.18191.82173.64176.78187.91176.57
P10177.82174.22160.50179.52167.78181.07182.16188.91179.26167.84175.19190.43182.58190.26172.03179.76185.53180.65
Table 14. Confirmation Results of Optimal Production Attributes in genetic algorithm.
Table 14. Confirmation Results of Optimal Production Attributes in genetic algorithm.
EXP.Resistance (Ohm × 10−3)Average Resistance (Ohm × 10−3)Standard DeviationSN
N1N2N3N4N5
Wafer #1177.00174.31177.16177.14178.44176.8101.51541.345
Wafer #2174.88178.07175.27177.20177.52176.5861.42441.871
Wafer #3177.03178.66177.20174.79177.94177.1241.45541.705
176.8401.46541.640
Table 15. Comparison results before and after proposed methodology.
Table 15. Comparison results before and after proposed methodology.
ComparisonDepo Time
(s)
Depo O2 Flow
(sccm)
ARC-LTO Etch Time
(s)
ARC-LTO Etch Pressure
(mTorr)
Ox-SiCO Etch Time
(s)
Ox-SiCO Gas Ratio
(%)
Polish Time
(s)
Average Resistance (Ohm × 10−3)|Resistance-Target|
ACDEFGI
Before improvement301605070703090191.1014.60
TSTM2515035951003085172.104.40
FMEA-TSTM-NNGA271514397912284176.840.34
Improvement 97.67%
Table 16. ANOVA table of average resistance of DRAM components from L12(211).
Table 16. ANOVA table of average resistance of DRAM components from L12(211).
ItemDegree of FreedomSum of SquareMean of SquareF Valuep Value
A139.33239.33243.370.096
B12.742.743.020.332
C111.96811.96813.20.171
D11.8881.8882.080.386
E1168.038168.038185.270.047
F1234.661234.661258.720.040
G16.8026.8027.50.223
H120.42120.42122.520.132
I153.76253.76259.270.082
J113.13513.13514.480.164
Error10.9070.907
Total11.0553.7
R-SqR-Sq(adj)
99.84%98.20%
Table 17. ANOVA table of SN of DRAM components from L12(211).
Table 17. ANOVA table of SN of DRAM components from L12(211).
ItemDegree of FreedomSum of SquareMean of SquareF Valuep Value
A11.9431.943103.630.062
B10.02280.02281.220.469
C114.840714.8407791.470.023
D15.6265.626300.040.037
E115.302115.3021816.080.022
F10.08070.08074.310.286
G18.32868.3286444.170.030
H10.19460.194610.380.192
I16.03176.0317321.680.035
J10.5770.57730.770.114
Error10.01880.0188
Total11.053.0
R-SqR-Sq(adj)
99.96%99.61%
Table 18. ANOVA table of average resistance of DRAM components from L18(21 × 37).
Table 18. ANOVA table of average resistance of DRAM components from L18(21 × 37).
ItemDegree of FreedomSum of SquareMean of SquareF Valuep Value
A2137.81268.905828.110.011
C20.4640.2320.090.912
D20.9790.48950.20.829
E297.15448.57719.810.019
F2174.78387.391635.650.008
G22.7521.37580.560.621
I2130.1765.08526.550.012
Error37.3552.4516
Total17551.468
R-SqR-Sq(adj)
98.67%92.44%
Table 19. ANOVA table of SN of DRAM components from L18(21 × 37).
Table 19. ANOVA table of SN of DRAM components from L18(21 × 37).
ItemDegree of FreedomSum of SquareMean of SquareF Valuep Value
A237.0618.538.050.062
C228.4414.226.180.086
D2129.38964.69428.10.011
E212.3036.1522.670.216
F211.3255.6632.460.233
G228.83314.4166.260.085
I216.1388.0693.50.164
Error36.9072.302
Total17270.395
R-SqR-Sq(adj)
97.45%85.52%
Table 20. The production attribute settings and the resistance value deviations of DRAM components after 10 times of GA.
Table 20. The production attribute settings and the resistance value deviations of DRAM components after 10 times of GA.
RunACDEFGI|Rs-Target|
#1 run27.929156.90941.90992.39698.62429.32681.7680.000131
#2 run27.014142.87635.03599.02894.05422.34288.2750.000441
#3 run27.668150.22035.26392.43691.23220.31783.0350.000900
#4 run26.264150.13043.08488.90992.29421.97183.9860.001989
#5 run25.891156.00241.06988.26698.24229.48785.5310.000139
#6 run27.042151.27438.51794.51697.59229.19085.7260.000665
#7 run27.182150.75443.12097.40791.46322.30183.7170.000120
#8 run28.948143.28741.36091.91492.03625.01580.2080.003525
#9 run26.652156.18736.95299.18287.43927.95284.4760.000272
#10 run28.877140.92544.87698.75093.28823.21082.1580.000188
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lin, C.-M.; Chen, S.-L. FMEA-TSTM-NNGA: A Novel Optimization Framework Integrating Failure Mode and Effect Analysis, the Taguchi Method, a Neural Network, and a Genetic Algorithm for Improving the Resistance in Dynamic Random Access Memory Components. Mathematics 2024, 12, 2773. https://doi.org/10.3390/math12172773

AMA Style

Lin C-M, Chen S-L. FMEA-TSTM-NNGA: A Novel Optimization Framework Integrating Failure Mode and Effect Analysis, the Taguchi Method, a Neural Network, and a Genetic Algorithm for Improving the Resistance in Dynamic Random Access Memory Components. Mathematics. 2024; 12(17):2773. https://doi.org/10.3390/math12172773

Chicago/Turabian Style

Lin, Chia-Ming, and Shang-Liang Chen. 2024. "FMEA-TSTM-NNGA: A Novel Optimization Framework Integrating Failure Mode and Effect Analysis, the Taguchi Method, a Neural Network, and a Genetic Algorithm for Improving the Resistance in Dynamic Random Access Memory Components" Mathematics 12, no. 17: 2773. https://doi.org/10.3390/math12172773

APA Style

Lin, C. -M., & Chen, S. -L. (2024). FMEA-TSTM-NNGA: A Novel Optimization Framework Integrating Failure Mode and Effect Analysis, the Taguchi Method, a Neural Network, and a Genetic Algorithm for Improving the Resistance in Dynamic Random Access Memory Components. Mathematics, 12(17), 2773. https://doi.org/10.3390/math12172773

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop