M.H.F. Sluiter
Please Note
20 records found
1
The thesis begins with an introduction to regression from a Bayesian perspective, establishing the connection between deterministic and probabilistic treatments, and covering models from linear to kernel and deep neural network regressions. This chapter provides the theoretical foundation for the remainder of the thesis by formalizing inference techniques, uncertainty quantification strategies, and model evaluation metrics.
Two core machine learning methods are developed to tackle the challenges of uncertainty estimation and multi-fidelity data fusion. First, a cooperative training scheme is proposed that combines a variance estimation neural network with a Bayesian mean neural network, enabling explicit disentanglement of aleatoric and epistemic uncertainties while improving mean prediction. The approach demonstrates strong scalability and generality across tasks and network architectures. Second, a practical multi-fidelity Bayesian learning framework is introduced, which fuses low- and high-fidelity data via a deterministic model, a transfer-learning module, and a Bayesian residual learner. This architecture balances expressiveness and computational efficiency, yielding robust predictions in both data-scarce and data-rich regimes.
To support data-driven mechanics, the proposed framework is extended to generalized constitutive modeling of history-dependent materials. A hierarchical learning scheme is developed that spans from single-fidelity deterministic networks to multi-fidelity Bayesian recurrent neural networks. This addresses two major limitations in data-driven modeling: the reliance on large, clean datasets and the lack of interpretability in neural network predictions.
Finally, the methodology is applied to the sustainable design of recycled composite polymers, where uncertainty arises from microstructural variability and the use of compatibilizers. Leveraging the cooperative training framework, both aleatoric and epistemic uncertainties are quantified, and a novel polymer design is optimized to achieve better expected performance and lower data variation.
Together, this thesis presents a unified and scalable framework for Bayesian learning under uncertainty and multi-fidelity conditions, with broad applicability in scientific computing, materials modeling, and sustainable engineering. The thesis concludes by summarizing key findings, discussing current limitations, and outlining future research directions. ...
The thesis begins with an introduction to regression from a Bayesian perspective, establishing the connection between deterministic and probabilistic treatments, and covering models from linear to kernel and deep neural network regressions. This chapter provides the theoretical foundation for the remainder of the thesis by formalizing inference techniques, uncertainty quantification strategies, and model evaluation metrics.
Two core machine learning methods are developed to tackle the challenges of uncertainty estimation and multi-fidelity data fusion. First, a cooperative training scheme is proposed that combines a variance estimation neural network with a Bayesian mean neural network, enabling explicit disentanglement of aleatoric and epistemic uncertainties while improving mean prediction. The approach demonstrates strong scalability and generality across tasks and network architectures. Second, a practical multi-fidelity Bayesian learning framework is introduced, which fuses low- and high-fidelity data via a deterministic model, a transfer-learning module, and a Bayesian residual learner. This architecture balances expressiveness and computational efficiency, yielding robust predictions in both data-scarce and data-rich regimes.
To support data-driven mechanics, the proposed framework is extended to generalized constitutive modeling of history-dependent materials. A hierarchical learning scheme is developed that spans from single-fidelity deterministic networks to multi-fidelity Bayesian recurrent neural networks. This addresses two major limitations in data-driven modeling: the reliance on large, clean datasets and the lack of interpretability in neural network predictions.
Finally, the methodology is applied to the sustainable design of recycled composite polymers, where uncertainty arises from microstructural variability and the use of compatibilizers. Leveraging the cooperative training framework, both aleatoric and epistemic uncertainties are quantified, and a novel polymer design is optimized to achieve better expected performance and lower data variation.
Together, this thesis presents a unified and scalable framework for Bayesian learning under uncertainty and multi-fidelity conditions, with broad applicability in scientific computing, materials modeling, and sustainable engineering. The thesis concludes by summarizing key findings, discussing current limitations, and outlining future research directions.
Precipitate development in steel
An in-depth study on formation and growth of a second solid phase
In the first chapter some theoretical background is presented to the topic.
Then in the second chapter an analytical approximation to the classical nucleation and growth model is presented.
In the third chapter a reference-free modified embedded atom method (RF-MEAM) potential is constructed for Fe, suitable for ferrite and austenite.
In chapter four the formation of Nb-C and Ti-C clusters in ferrite is studied, here formation paths with a monotone decrease in Gibbs energy are found. Moreover, the resulting clusters remain stable up to 1100K.
In the fifth chapter other group IV and V transition metal-carbon systems in ferrite are studied in small cells using density functional theory. Showing similar energy patterns as previously found for Nb and Ti. Furthermore the transformation from a metal-carbon cluster to a carbide precipitate is studied, where the matrix strain is found to be the driving force behind the transformation.
in the sixth and last chapter some recommendations are presented. ...
In the first chapter some theoretical background is presented to the topic.
Then in the second chapter an analytical approximation to the classical nucleation and growth model is presented.
In the third chapter a reference-free modified embedded atom method (RF-MEAM) potential is constructed for Fe, suitable for ferrite and austenite.
In chapter four the formation of Nb-C and Ti-C clusters in ferrite is studied, here formation paths with a monotone decrease in Gibbs energy are found. Moreover, the resulting clusters remain stable up to 1100K.
In the fifth chapter other group IV and V transition metal-carbon systems in ferrite are studied in small cells using density functional theory. Showing similar energy patterns as previously found for Nb and Ti. Furthermore the transformation from a metal-carbon cluster to a carbide precipitate is studied, where the matrix strain is found to be the driving force behind the transformation.
in the sixth and last chapter some recommendations are presented.
This thesis develops physics-guided machine learning approaches to accelerate the design, modeling, and interpretation of materials across scales and material classes. It introduces neural network architectures for metamaterial design that learn directly from data while remaining physically consistent. These models retrieve structural designs that achieve specified target properties instantly, enabling rapid exploration of the structure–property landscape and supporting scenarios with multiple design goals more efficiently than traditional optimization methods. For material modeling, this work addresses the limitation that stress fields are not directly accessible in experiments. The proposed frameworks instead learn constitutive laws from measurable quantities such as displacements and forces, while preserving essential physical principles such as thermodynamic consistency. This work further demonstrates how data-driven approaches reveal previously unknown structure–property relationships, such as the role of dangling atomic masses in the thermal transport of crystalline nanoporous materials. Finally, it introduces a chemistry-constrained generative framework that proposes synthesizable, diverse, and novel molecular structures using limited data, while providing interpretable representations of the molecular generation.
Together, these contributions establish that physics-guided machine learning can complement and extend traditional materials science by delivering reliable, interpretable, and generalizable solutions to longstanding challenges in design and modeling. ...
This thesis develops physics-guided machine learning approaches to accelerate the design, modeling, and interpretation of materials across scales and material classes. It introduces neural network architectures for metamaterial design that learn directly from data while remaining physically consistent. These models retrieve structural designs that achieve specified target properties instantly, enabling rapid exploration of the structure–property landscape and supporting scenarios with multiple design goals more efficiently than traditional optimization methods. For material modeling, this work addresses the limitation that stress fields are not directly accessible in experiments. The proposed frameworks instead learn constitutive laws from measurable quantities such as displacements and forces, while preserving essential physical principles such as thermodynamic consistency. This work further demonstrates how data-driven approaches reveal previously unknown structure–property relationships, such as the role of dangling atomic masses in the thermal transport of crystalline nanoporous materials. Finally, it introduces a chemistry-constrained generative framework that proposes synthesizable, diverse, and novel molecular structures using limited data, while providing interpretable representations of the molecular generation.
Together, these contributions establish that physics-guided machine learning can complement and extend traditional materials science by delivering reliable, interpretable, and generalizable solutions to longstanding challenges in design and modeling.
Anisotropic Stress, Plasticity, and Microstructural Evolution in Crystalline Materials
From Grain Boundaries to Nanostructures
Despite the promising performance of RRAM, a high forming voltage potentially hinders the widespread application of the device. This thesis aims to diminish and eliminate the forming voltage. To achieve this, different metals were inserted between the insulator layer and the bottom electrode of the RRAM, serving as an interface metal layer. The interface metal was expected to introduce oxygen vacancies to the insulator, thereby decreasing the forming voltage. Advanced nanofabrication processes were employed in the cleanroom, and a related recipe was developed. The influence of layer thickness and device area was also studied to gain a comprehensive understanding. Among all the samples, Ru-based devices were observed to be forming-free.
Data analysis methods were applied to model the data, with the random forest method found to be the most suitable, achieving an accuracy of 82.4%. The model was verified by measurements of 10 nm Ru-based devices. Feature importance was then calculated to interpret the model. The four most important features determining the forming voltage are the thickness, standard electrode potential, area, and work function of the interface metal. This work adopts a new approach to eliminating the forming voltage, not only providing a forming-free device but also offering a guideline for future research on forming voltage. ...
Despite the promising performance of RRAM, a high forming voltage potentially hinders the widespread application of the device. This thesis aims to diminish and eliminate the forming voltage. To achieve this, different metals were inserted between the insulator layer and the bottom electrode of the RRAM, serving as an interface metal layer. The interface metal was expected to introduce oxygen vacancies to the insulator, thereby decreasing the forming voltage. Advanced nanofabrication processes were employed in the cleanroom, and a related recipe was developed. The influence of layer thickness and device area was also studied to gain a comprehensive understanding. Among all the samples, Ru-based devices were observed to be forming-free.
Data analysis methods were applied to model the data, with the random forest method found to be the most suitable, achieving an accuracy of 82.4%. The model was verified by measurements of 10 nm Ru-based devices. Feature importance was then calculated to interpret the model. The four most important features determining the forming voltage are the thickness, standard electrode potential, area, and work function of the interface metal. This work adopts a new approach to eliminating the forming voltage, not only providing a forming-free device but also offering a guideline for future research on forming voltage.
...
Following ideas from the neuroscience literature, we create task-specific regions in the network, i.e. subnetworks, to encode information there. Thus, some parameters are responsible for solving this task, which mitigates forgetting compared to conventional training where the trainable parameters are simultaneously assigned to all tasks. A proper subnetwork should be then selected by the algorithm to make a prediction or information about the correct subnetwork must be given by the user. The subnetworks can share some connections to transfer knowledge between each other and facilitate future learning.
In the first part of the thesis, we describe the proposed methodology: task-specific subnetwork creation during training and the proper subnetwork selection during inference stages. We examine different subnetwork prediction strategies outlining their advantages and disadvantages. We validate the proposed algorithms on a series of well-known image datasets in computer vision in classification and semantic segmentation tasks. The proposed solution significantly outperforms current state-of-the-art methods by 10-20\% of accuracy.
The second part of the thesis illustrates the benefits of cooperative learning via continual learning in physical sciences and solid mechanic examples. We demonstrate that by sharing parameters, the following subnetwork can be trained either with lower prediction error, requiring fewer training data points, or both, compared to conventional training with one network per task. Importantly, the model does not forget any of the acquired knowledge since once a parameter is assigned to a subnetwork, it is not changed when training new tasks. We would like to highlight the potential importance of further development of continual learning methods in engineering to improve the generalization capabilities of the models.
The thesis concludes by discussing the main results and findings. We also outline the main limitations of the work and directions for improvement. Further development of continual learning models will lead to more advanced artificial intelligence systems that should contribute to solving a wider range of problems. ...
Following ideas from the neuroscience literature, we create task-specific regions in the network, i.e. subnetworks, to encode information there. Thus, some parameters are responsible for solving this task, which mitigates forgetting compared to conventional training where the trainable parameters are simultaneously assigned to all tasks. A proper subnetwork should be then selected by the algorithm to make a prediction or information about the correct subnetwork must be given by the user. The subnetworks can share some connections to transfer knowledge between each other and facilitate future learning.
In the first part of the thesis, we describe the proposed methodology: task-specific subnetwork creation during training and the proper subnetwork selection during inference stages. We examine different subnetwork prediction strategies outlining their advantages and disadvantages. We validate the proposed algorithms on a series of well-known image datasets in computer vision in classification and semantic segmentation tasks. The proposed solution significantly outperforms current state-of-the-art methods by 10-20\% of accuracy.
The second part of the thesis illustrates the benefits of cooperative learning via continual learning in physical sciences and solid mechanic examples. We demonstrate that by sharing parameters, the following subnetwork can be trained either with lower prediction error, requiring fewer training data points, or both, compared to conventional training with one network per task. Importantly, the model does not forget any of the acquired knowledge since once a parameter is assigned to a subnetwork, it is not changed when training new tasks. We would like to highlight the potential importance of further development of continual learning methods in engineering to improve the generalization capabilities of the models.
The thesis concludes by discussing the main results and findings. We also outline the main limitations of the work and directions for improvement. Further development of continual learning models will lead to more advanced artificial intelligence systems that should contribute to solving a wider range of problems.
DFT calculation of NMR parameters for the K doped CsPbF3 solid-state ionic conductor
A small step toward a green vigor
The road to intelligent asphalt concrete mixture design
A Data driven analysis of common asphalt concrete property prediction methods and a solution to the inverse problem
Experimental investigation on low-temperature electrical resistivity of silver thin films
The effects of fabrication, design, purity and storage
Therefore, this project investigates hydrogen diffusivity in equimolar TiVZrHfNb and the influence of hydrogen concentration on hydrogen diffusivity to gain better understanding of hydrogenation kinetics. The selected HEA has been found to absorb the highest amount of hydrogen (2.5 H/M) among other HEAs.
The investigation was done by a computational approach using ab initio molecular dynamics. BCC and face-centered cubic (FCC) supercells with different hydrogen concentrations (H/M = 0.2, 0.8, 1.4, 2, 2.4) were simulated at a temperature range of 773 – 973 K. At the same time, experimental electrochemical hydrogen charging using chronoamperometry and cyclic voltammetry was performed in order to compare computational and experimental values of hydrogen diffusivity.
The electrochemical hydrogen charging did not result in hydrogen absorption, most probably due to the passivation of the sample surface.
From the simulation results, the values of activation energy and pre-exponential factor were estimated to be in the range of 0.26 – 0.48 eV and 0.73 – 2.95 x 10-7 m2/s, respectively. Hydrogen diffusivity was found to be higher in BCC than in FCC. In BCC the hydrogen diffusivity slowly decreases linearly with increasing H/M. In the case of FCC, the hydrogen diffusivity was found to be the highest at 2.4 H/M while at 2 H/M the diffusivity was the lowest. The analysis of hydrogen occupation at 2 H/M shows that most of the hydrogen atoms are trapped inside tetrahedral sites. It is possible that the hydrogen occupation in the tetrahedral sites results in the optimum hydrogen distribution where the repulsive interaction between hydrogen atoms is the lowest. At concentrations above 2 H/M, an additional repulsive force between hydrogen atoms seems to contribute to the increase of hydrogen diffusivity. ...
Therefore, this project investigates hydrogen diffusivity in equimolar TiVZrHfNb and the influence of hydrogen concentration on hydrogen diffusivity to gain better understanding of hydrogenation kinetics. The selected HEA has been found to absorb the highest amount of hydrogen (2.5 H/M) among other HEAs.
The investigation was done by a computational approach using ab initio molecular dynamics. BCC and face-centered cubic (FCC) supercells with different hydrogen concentrations (H/M = 0.2, 0.8, 1.4, 2, 2.4) were simulated at a temperature range of 773 – 973 K. At the same time, experimental electrochemical hydrogen charging using chronoamperometry and cyclic voltammetry was performed in order to compare computational and experimental values of hydrogen diffusivity.
The electrochemical hydrogen charging did not result in hydrogen absorption, most probably due to the passivation of the sample surface.
From the simulation results, the values of activation energy and pre-exponential factor were estimated to be in the range of 0.26 – 0.48 eV and 0.73 – 2.95 x 10-7 m2/s, respectively. Hydrogen diffusivity was found to be higher in BCC than in FCC. In BCC the hydrogen diffusivity slowly decreases linearly with increasing H/M. In the case of FCC, the hydrogen diffusivity was found to be the highest at 2.4 H/M while at 2 H/M the diffusivity was the lowest. The analysis of hydrogen occupation at 2 H/M shows that most of the hydrogen atoms are trapped inside tetrahedral sites. It is possible that the hydrogen occupation in the tetrahedral sites results in the optimum hydrogen distribution where the repulsive interaction between hydrogen atoms is the lowest. At concentrations above 2 H/M, an additional repulsive force between hydrogen atoms seems to contribute to the increase of hydrogen diffusivity.
Adsorption of Phosphorus on a Ru(0001) surface
A density functional theory study
The aim of this project was to test regimes simulating steady state operation and high frequency and temperature transient pulses called ELMs (Edge Localized Modes) striking the divertor. Recent research has shown that the degradation and behavior of tungsten under these transient conditions does not consistently follow the expected parameters characterized in the literature. According to it, recrystallization, grain growth, and crack formation seem to be suppressed by the plasma loading under these regimes, thus a new understanding of the material behavior for these circumstances must be developed. To do this, ITER grade tungsten samples were subjected to a hydrogen plasma beam at DIFFER’s Magnum-PSI with temperatures at the strike point ranging from ~1000 to ~1500 °C and high frequency pulses that increased the surface temperature by ~200 to ~300 °C above the steady state temperature. The surface thermal shock response to the plasma pulses was characterized by means of infrared and pyrometer readings at the samples’ surface during exposure. Temperature and power density calculations were correlated with identified damage morphologies on the targets and a damage map for the experiments was elaborated, which showed that the most severe damage (cracks and crack networks) begin to appear in the range of the measured recrystallization temperature of the samples, which was lower than expected.
Using Vickers hardness, the recovery and recrystallization kinetics of the material were characterized by means of logarithmic decay and a modified version of JMAK recrystallization kinetics that includes an incubation time for the onset of recrystallization. Recrystallization kinetics were found to accelerate as the hydrogen exposure progresses, thus yielding lower effective activation energies for recrystallization when comparing furnace one hour exposures, plasma one hour exposures, and plasma four hour exposures. This pointed to the presence of hydrogen actively reducing the activation energy for self-diffusion. Simulations of hydrogen diffusion were performed to test this hypothesis, and even though the total concentration is low, given that the high experimental temperature does not permit trapping of the hydrogen, the diffusing atoms may still play a role in accelerating the recrystallization kinetics.
Based on the results of this research, it is proposed that interstitially diffusing hydrogen segregating to voids or grain boundaries is modifying the behavior of surrounding tungsten crystal lattice. Specifically, the mobility of the grain boundaries may be increasing because the hydrogen’s presence would be promoting the creation of ledges in the grain boundary resulting in an overall free energy reduction for grain boundary diffusion. This localized defect formation would require a lower concentration of hydrogen than that required for solute drag or other suppression mechanisms, and might be mechanism behind the behavior observed in this work’s experiments. ...
The aim of this project was to test regimes simulating steady state operation and high frequency and temperature transient pulses called ELMs (Edge Localized Modes) striking the divertor. Recent research has shown that the degradation and behavior of tungsten under these transient conditions does not consistently follow the expected parameters characterized in the literature. According to it, recrystallization, grain growth, and crack formation seem to be suppressed by the plasma loading under these regimes, thus a new understanding of the material behavior for these circumstances must be developed. To do this, ITER grade tungsten samples were subjected to a hydrogen plasma beam at DIFFER’s Magnum-PSI with temperatures at the strike point ranging from ~1000 to ~1500 °C and high frequency pulses that increased the surface temperature by ~200 to ~300 °C above the steady state temperature. The surface thermal shock response to the plasma pulses was characterized by means of infrared and pyrometer readings at the samples’ surface during exposure. Temperature and power density calculations were correlated with identified damage morphologies on the targets and a damage map for the experiments was elaborated, which showed that the most severe damage (cracks and crack networks) begin to appear in the range of the measured recrystallization temperature of the samples, which was lower than expected.
Using Vickers hardness, the recovery and recrystallization kinetics of the material were characterized by means of logarithmic decay and a modified version of JMAK recrystallization kinetics that includes an incubation time for the onset of recrystallization. Recrystallization kinetics were found to accelerate as the hydrogen exposure progresses, thus yielding lower effective activation energies for recrystallization when comparing furnace one hour exposures, plasma one hour exposures, and plasma four hour exposures. This pointed to the presence of hydrogen actively reducing the activation energy for self-diffusion. Simulations of hydrogen diffusion were performed to test this hypothesis, and even though the total concentration is low, given that the high experimental temperature does not permit trapping of the hydrogen, the diffusing atoms may still play a role in accelerating the recrystallization kinetics.
Based on the results of this research, it is proposed that interstitially diffusing hydrogen segregating to voids or grain boundaries is modifying the behavior of surrounding tungsten crystal lattice. Specifically, the mobility of the grain boundaries may be increasing because the hydrogen’s presence would be promoting the creation of ledges in the grain boundary resulting in an overall free energy reduction for grain boundary diffusion. This localized defect formation would require a lower concentration of hydrogen than that required for solute drag or other suppression mechanisms, and might be mechanism behind the behavior observed in this work’s experiments.