Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular dynamic (MD) simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Ion Size Effects on the Thermodynamic, Kinetic, and Mechanical Properties during Ion Exchange in Solid-State Electrolytes

Ion exchange offers a pathway to impose residual compressive stresses to mitigate the electro-chemo-mechanical cracking of solid-state electrolytes such as lithium lanthanum zirconium oxide. This study uses a coupled multiscale framework (integrating density functional theory (DFT), molecular dynamics (MD), and continuum modeling) to examine how exchange ion size influences stress, diffusion, fracture toughness, and electronic properties. Larger isovalent ions (Na + , Ag + , K + ) were exchanged with Li + , with DFT confirming their preference for octahedral 96h sites and a linear relationship between ion size and chemical free expansion coefficient. MD simulations reveal stress and concentration effects on exchange ion diffusivity at elevated temperatures, with Na + and Ag + maintaining favorable mobility while K + showing concentration-dependent clustering. Continuum modeling predicts the range of fracture strength improvements and the required ion exchange concentration profile. It was shown that a 5% surface exchange concentration can induce ∼0.6 GPa of surface compressive stress using Na + and ∼1.0 GPa of surface compressive stress using Ag + . On the other hand, larger ion exchange species may penalize Li + diffusivity by increasing the activation volume and activation energy. Interestingly, Na + has a negligible penalty on Li-ion diffusivity. The room temperature Li + ion diffusivity is reduced by ∼40% with Ag + ion exchange. Electronic band structure analysis shows no size-dependent change in the bandgap, though Ag + introduces localized defect states near the valence band maximum. This study highlights ion size as a key factor in optimizing LLZO properties, offering a framework to improve the solid-state battery performance.

Jagad, Harsh D. [Brown Univ., Providence, RI (Unit↗
DOE OSTI2022

Exploration of the Application of Machine Learning to the improvement of Interatomic Potentials

Current methods for atomistic simulations of material systems suffer from limitations which restrict the ability of the simulation to correctly characterize certain material behavior and physical phenomena. Small scale ab initio (AIMD) modeling is highly accurate but is computationally expensive. Classical molecular dynamics (CMD) simulations use interatomic potentials (IPs) to describe larger systems at a reduced computational cost, but with a reduced accuracy. Creating more robust IPs enable more precise CMD simulations. In this thesis, the application of a specific machine learning process,artifical neural networks (ANNs), to the improvement of IPs and MD simulation is discussed.

O'Connor, Sean↗
NASA NTRS1991

Disordering and H-embrittlement of Pb(110) surface using embedded atom method and molecular dynamics

The dynamics and structure of the crystal melt interface are of great importance in crystal growth, melting, and other properties of the solid-liquid interface. The Embedded Atom Method (EAM), a semiempirical method developed by Daw and Baskes for determining the energetics of atoms in a bulk environment, is a useful method in the computer simulation. The EAM functions of PB in conjunction with the MD code are employed to study melting and disordering of the Pb(100) surface. The EAM functions are fitted to the bulk experimental data at zero temperature, and it is not clear if it will behave correctly at higher temperatures. In fact, it is shown that use of a similar type of potential at higher temperatures might result in a negative coefficient of thermal expansion. The primary purpose here is to perform MD simulation combined with the EAM functions of PB to study melting of the Pb(100) surface. This provides an accurate test of the EAM functions at higher temperatures. In particular, we calculate the two dimensional structure factor, the two dimensional distribution function, and density profile as a function of temperature up to the melting point of lead.

Karimi, Majid↗
DOE OSTI2021

Bayesian automated weighting of aggregated DFT, MD, and experimental data for candidate thermodynamic models of aluminum with uncertainty quantification

Atomic-scale modeling methods such as density functional theory (DFT) and molecular dynamics (MD) can predict the thermodynamic properties of materials at a lower cost than experimental measurements. However, their regular usage in thermodynamic model construction is hampered by the lack of quantitative agreement with experimental measurements and the lack of uncertainty estimates on the data. To make regular usage of this atomistic simulation data, it is important to assess whether the atomistic simulation datasets, by themselves or in combination with experimental measurements, result in the same physics-informed models best supported by experimental measurements alone. Here, models of aluminum thermodynamic properties are discussed using three data sources: atomistic calculations (DFT and MD), experiments, and a combination of atomistic calculations and experiments. The study shows that, after ensuring self-consistency in predicting key invariant points, both experimental measurements and atomistic calculations can significantly contribute to an optimal model.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Zinc binding alters the conformational dynamics and drives the transport cycle of the cation diffusion facilitator YiiP

YiiP is a secondary transporter that couples Zn 2+ transport to the proton motive force. Structural studies of YiiP from prokaryotes and Znt8 from humans have revealed three different Zn 2+ sites and a conserved homodimeric architecture. These structures define the inward-facing and outward-facing states that characterize the archetypal alternating access mechanism of transport. To study the effects of Zn 2+ binding on the conformational transition, we use cryo-EM together with molecular dynamics simulation to compare structures of YiiP from Shewanella oneidensis in the presence and absence of Zn 2+ . To enable single-particle cryo-EM, we used a phage-display library to develop a Fab antibody fragment with high affinity for YiiP, thus producing a YiiP/Fab complex. To perform MD simulations, we developed a nonbonded dummy model for Zn 2+ and validated its performance with known Zn 2+ -binding proteins. Using these tools, we find that, in the presence of Zn 2+ , YiiP adopts an inward-facing conformation consistent with that previously seen in tubular crystals. After removal of Zn 2+ with high-affinity chelators, YiiP exhibits enhanced flexibility and adopts a novel conformation that appears to be intermediate between inward-facing and outward-facing states. This conformation involves closure of a hydrophobic gate that has been postulated to control access to the primary transport site. Comparison of several independent cryo-EM maps suggests that the transition from the inward-facing state is controlled by occupancy of a secondary Zn 2+ site at the cytoplasmic membrane interface. This work enhances our understanding of individual Zn 2+ binding sites and their role in the conformational dynamics that govern the transport cycle.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A prediction of the thermodynamic, thermophysical, and mechanical properties of CrTaO 4 from first principles

It is reported that the self-forming CrTaO 4 oxide scale can protect refractory high-entropy alloys from oxidation, superior to Cr 2 O 3 . In this paper, the phase stability, mechanical, and thermal properties of three polymorphous phases of CrTaO 4 are systematically investigated from first-principles density functional theory calculations. The mechanical properties predicted using the strain–energy methods indicated that all three phases are mechanically stable. The temperature dependence of elastic constants and polycrystalline moduli of three phases demonstrated the thermal softening as temperature increase. The Helmholtz-free energies as a function of volume and temperature are derived from phonon dispersions within the quasi-harmonic approximation at six strained volumes. The calculated apparent bulk coefficients of thermal expansion of these three phases are evaluated, the highest value approximately 13.4× 10 –6 K –1 within a temperature range of 500–2000 K for the rutile I4 1 md phase. The lattice thermal conductivity calculated by the Debye–Callaway model suggested that the rutile type I4 1 md phase has the lowest value of approximately 2.1 W/m/K at 1800 K. The other two phases, C2/m and P2/c, exhibit higher values due to relatively lower Grüneisen parameters and larger phonon velocities. The melting point of CrTaO 4 is predicted to be between 1975 and 2449 K using ab initio molecular dynamics simulations. Finally, this work provides a comprehensive theoretical understanding of the thermodynamic, mechanical, and thermal properties for the new material CrTaO 4 and serves as an example of a viable computational design strategy for improved oxidation resistance of refractory alloys at high temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A First Principles Approach to Spectral Phonon Transport in Heterostructures

Understanding thermal transport across interfaces which give rise to a thermal resistance (also known as Kapitza resistance) is a critical issue affecting the development of nanotechnologies. Much modern and emergent nanotechnology consist of adjacent materials, and phonon mediated heat transfer governs thermal behavior across internal interfaces in these devices. The physics of thermal transport in solids are governed both by phenomena occurring at the atomic scale and interactions with the material's microstructure. The forecasting of fundamental quantities such as temperature, heat flux and thermal conductivity typically employs the semi-classical Boltzmann transport equation to predict the macroscopic behavior of materials in terms of the microscopic dynamics of its heat carriers. Kapitza resistance was first discovered in liquid helium experiments and has led to a fundamental research thrust in micro and nano-scale heat transport, the behavior of thermal carriers across internal interfaces. Thermal interfacial resistance (TIR) is a widely studied phenomenon, first engaged by Swartz and Pohl through their development of the acoustic and diffuse mismatch methods, then continued through myriad efforts with varying methods and approaches in an attempt to resolve carrier behavior at thermal interfaces. Many of the fundamental approaches to TIR have been at the nanoscale, and research is conducted with molecular dynamics (MD) and density functional theory (DFT) methods. The limitations of these methods is system size; atomistic methods tend to be limited to system sizes of 100,000 atoms or less. Larger length-scale methods have also been pursued, based on the principles of acoustic or diffuse mismatch, but not all include simulation of TIR using a full phonon band spectrum, or temperature dependent methods. Our approach to enabling phonon transport in layered materials draws upon our previous work of demonstrating spectrally coupled phonon transport in homogeneous and heterogeneous materials. We use a semi-analytical approach in which the Bose-Einstein (B-E) statistics set the strength of the phonon radiance in a frequency group, but the B-E statistics are informed with information from the transport system. The B-E statistics in a single frequency group feels the influence of all the groups through the spatial temperature. We also include a new field term which is an indicator of the amount of non-equilibrium behavior of the phonon spectrum---this is added to the phonon source term in all groups to ensure closure and conservation of energy, as the phonon groups in the transport system and the analytical systems are coupled. This work builds upon our previous approach by adding a phonon coupling term at an internal interface, using the principles of the DMM through transmission and reflection coefficients. In this work, the coefficients are determined through computing a common temperature at the interface, influenced by the phonon band structure of both materials, in effect, providing mixing between the two material systems and using the common temperature to set the strength of the phonon radiance at the boundaries on either side of the interface. Our approach uses material properties computed along various crystallographic orientations, and while some isotropy is built into the interface condition, the material properties weight the phonon distributions in the proper crystalline direction. Greater resolution of phonon behavior in proximity to an interface, and more accurate predictions of TIR are obtained. While it is true the assumption of diffuse mismatch can yield inconsistent results compared to experiment especially at low temperatures, this work focuses on room temperature and beyond effects, for future applications in nuclear fuel, or thermoelectric devices; a modified mismatch approach may be feasible if applied properly. Additionally, our methods focus on bridging mesoscale to engineering scale

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling bicarbonate formation in an alkaline solution with multi-level quantum mechanics/molecular dynamics simulations

Understanding carbonate speciation and how it may be modulated is essential for the advancement of carbon dioxide (CO 2 ) capture and storage technologies, which often rely on the transformation of CO 2 into carbonate, e.g. via the formation of carbonate minerals. To date, few atomic-level, quantum-mechanics-based simulations have been carried out to characterize how carbonic acid (H 2 CO 3 ) and bicarbonate ($HCO^{-}_{3}$) form in aqueous solution, and how pH affects this process. Recently, Martirez and Carter utilized rare-event sampling density functional theory molecular dynamics simulations in combination with multi-level embedded correlated wavefunction theory, thus accounting for both solvent dynamics and electron correlation accurately, to elucidate the mechanism of H 2 CO 3 formation in neutral solution (J. Am. Chem. Soc., 145, 12561, 2023). Here, we perform a complementary simulation using the same method to map out the energetics of $HCO^{-}_{3}$ formation from dissolved CO 2 in basic solution. We find that, as in H 2 CO 3 formation, including water dynamics is important to obtain an accurate prediction of the energetics for the aforementioned reaction. Furthermore, only with MD did we identify the correct pathway for the reaction, in which water – not hydroxide – acts as the initial nucleophile and only at the transition state does it lose a proton.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Molecular Insights into Geochemical Reactions of Iron-Bearing Minerals: Implications for Hydrogen Geo-Storage

Here, this study investigates the reaction of hydrogen (H 2 ) with pyrite (FeS 2 ), focusing on how temperature and the presence of water influence the reaction pathways and kinetics. Utilizing computational molecular simulations and kinetic analyses, we explore the impact of these factors on the formation of hydrogen sulfide (H 2 S) and related species. First, grand canonical Monte Carlo/molecular dynamics (GCMC/MD) simulations reveal that physical H 2 adsorption occurs in distinct layers on the pyrite surface. In addition, increased temperatures reduce the absolute adsorption capacities. Reactive MD simulations demonstrate that H 2 interacts differently with pyrite under varying conditions. At 298 K, H 2 reacts with pyrite to form HS – , leading to the formation of HS – through covalent bonding with sulfur of pyrite. However, no H 2 S is produced at this temperature, suggesting that a kinetic barrier (i.e., activation energy) may prevent this reaction. At higher temperatures, H 2 S production significantly increases. The presence of water introduces additional complexity to the reaction mechanism. Unlike dry conditions, water enhances H 2 S generation, even at low temperatures. Water also facilitates the formation of additional products, such as SOH, indicating a more intricate chemical environment on the pyrite surface. Our findings identify the association of HS – ions to form H 2 S as the rate-limiting step, with temperature influencing this process. This finding suggests that while the presence of water can create a more dynamic reaction environment, the overall mechanisms leading to H 2 S formation remain consistent. These outcomes suggest the need for developing targeted strategies to manage and control H 2 S emissions within the context of underground hydrogen storage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Effect of Solvent on the Interaction of Lignin with a Zeolite Nanosheet in the Condensed Phase

Lignin is the essential building block of lignocellulosic biomass, an excellent renewable source of different aromatic monomers for the polymer and biofuel industry. The depolymerization of lignin into value-added chemicals and fuels through the catalytic process poses a significant challenge due to the complex structure of lignin. Understanding lignin’s conformational diversity and dynamics in the liquid phase is crucial to designing an effective depolymerization process. Here, we conducted all-atom molecular dynamics simulations to understand the conformation and dynamics of softwood lignin on the all-silica zeolite nanosheet based on the MFI topology in a binary mixture of water–methanol at three different molar compositions (0%, 50%, and 100% methanol). We observed that the methanol–surface interaction is stronger than the water–surface interaction, and methanol readily diffused into the MFI core. Lignin surface contacts decrease with increasing methanol composition due to higher solubility and dynamics. Lignin dynamics on the surface in neat water is an order of magnitude smaller than methanol. Here, we also found that lignin adopts a slightly extended conformation when it stays on the surface than in the bulk solution phase for the pure water case, whereas for pure methanol and the binary solution structures are statistically similar.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Free Energy and Flexibility Analysis of Autoinhibited Human BRAF

The RAF serine/threonine protein kinases function as direct effectors of RAS in the intracellular transmission of extracellular growth signals, and they are key targets for drug discovery, given the high incidence of oncogenic mutations in RAF and other components of this signaling pathway. In its inactive state, RAF is held in an autoinhibited conformation in the cytosol through a combination of intramolecular interactions and binding to a regulatory 14−3−3 protein dimer. Activation of RAF is initiated by its interaction with membrane-localized GTP-bound RAS, which induces conformational changes that release RAF from its autoinhibited state. However, the molecular mechanisms governing RAF activation remain incomplete, largely due to the challenges in experimentally capturing the intermediate conformational states in this process. To address this gap, we developed a comprehensive all-atom model of BRAF based on existing cryo-EM structures. Using this model, we performed extensive molecular dynamics simulations to evaluate the stability and free energy landscape of autoinhibited BRAF in solution. Our analysis reveals conformational flexibility within the autoinhibited complex, suggesting that this dynamic behavior may play a role in facilitating BRAF activation upon engagement with the membrane-bound RAS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Limiting velocities and transonic dislocations in Mg

To accurately predict the mechanical response of materials, especially at high strain rates, it is important to account for dislocation velocities in these regimes. Under these extreme conditions, it has been hypothesized that dislocations can move faster than the speed of sound. However, the presence of such dislocations remains elusive due to challenges associated with measuring these experimentally. In this work, molecular dynamics simulations were used to investigate the dislocation velocities for the basal edge, basal screw, prismatic edge, and prismatic screw dislocations in Mg in the sub-, trans-, and supersonic regimes. Also, our results show that only prismatic edge dislocations achieve supersonic velocities. Furthermore, this work demonstrates that the discrepancy between the theoretical limiting velocity and the MD results for Mg is due to its sensitivity to large hydrostatic stress around the dislocation core, which was not the case for fcc metals such as Cu.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Towards a unified nonlocal, peridynamics framework for the coarse-graining of molecular dynamics data with fractures

Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. In this work, we propose a learning framework to extract a peridynamics model as a mesoscale continuum surrogate from MD simulated material fracture data sets. Firstly, we develop a novel coarse-graining method, to automatically handle the material fracture and its corresponding discontinuities in the MD displacement data sets. Inspired by the weighted essentially non-oscillatory (WENO) scheme, the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil, then reconstruct the coarse-grained material displacement field as the piecewise smooth solutions containing discontinuities. Then, based on the coarse-grained MD data, a two-phase optimization-based learning approach is proposed to infer the optimal peridynamics model with damage criterion. In the first phase, we identify the optimal nonlocal kernel function from the data sets without material damage to capture the material stiffness properties. Then, in the second phase, the material damage criterion is learnt as a smoothed step function from the data with fractures. As a result, a peridynamics surrogate is obtained. As a continuum model, our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training, and hence allows for substantial reductions in computational cost compared with MD. We illustrate the efficacy of the proposed approach with several numerical tests for the dynamic crack propagation problem in a single-layer graphene. Our tests show that the proposed data-driven model is robust and generalizable, in the sense that it is capable of modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Clustering Effects on the Structure of Ionomer Solutions: A Combined SANS and Simulations Study

Ionic assemblies, or clusters, determine the structure and dynamics of ionizable polymers and enable their many applications. Fundamental to attaining well-defined materials is controlling the balance between van der Waals interactions that govern the backbone behavior and the forces that drive the formation of ionic clusters. Here, using small angle neutron scattering and fully atomistic molecular dynamics simulations, the structure of a model ionomer, sulfonated polystyrene in toluene solutions, was investigated as the cluster cohesion was tweaked by the addition of ethanol. The static structure factor was measured by both techniques and correlated with the size of the ionic clusters as the polymer concentration was varied. The conjunction of SANS results and molecular insight from MD simulations enabled the determination of the structure in these inhomogeneous networks on multiple length scales. Additionally, we find that across the entire concentration range studied, a network driven by the formation of ionic clusters was formed, where the size of the clusters drives the inhomogeneity of these systems. Tweaking the ionic clusters through the addition of ethanol impacts the packing of the sulfonated groups, their shape, and their size distribution, which in turn, affects the structure of these networks.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Predicted Fracture Tendency of Naturally Occurring Aluminum Surface Coatings under Tensile Loading

Naturally occurring coatings on aluminum metal, such as its oxide or hydroxide, serve to protect the material from corrosion. Understanding the conditions under which these coatings mechanically fail is therefore expected to be an important aspect of predictive models for aluminum component lifetimes. To this end, we develop and apply a molecular dynamics (MD) modeling framework for conducting tension tests that is capable of isolating factors governing the mechanical strength as a function of coating chemistry, defect morphology, and variables associated with the loading path. We consider two representative materials, including γ-Al 2 O 3 and γ-Al(OH) 3 (i.e., oxide and hydroxide), both of which form readily as aluminum surface coatings. Our results indicate that defects have a significant bearing on the strength of aluminum oxide, with grain boundaries serving to reduce the strain at failure from ε zz = 0.300 to 0.219, relative to perfect single crystal. Our simulations also predict that porosity lowers the elastic stiffness and yield strength of the oxide. Relative to perfect crystal, we find porosity factors of 5%, 10% and 20% decrease the yield stress by 26%, 36% and 53%, respectively. MD predicts that perfect hydroxide and oxide single crystal have respective strains at failure of 0.08 and 0.31 under tensile uniaxial strain loading, and that the corresponding yield stresses are respectively 1.6 and 11.1 GPa. These data indicate that the hydroxide is substantially more susceptible to mechanical failure than the oxide. Our results, coupled with literature findings that indicate hot and humid conditions favor formation of hydroxide and defective oxide coatings, indicate the potential for a complicated dependence of aluminum corrosion susceptibility and stress corrosion cracking on aging history.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Insights on the dissolution of water in an albite melt at high pressures and temperatures from a direct structural analysis

The water dissolution mechanism in silicate melts under high pressures is not well understood. Here we present the first direct structure investigation of a water-saturated albite melt to monitor the interactions between water and the network structure of silicate melt at the molecular level. In situ high-energy X-ray diffraction was carried out on the NaAlSi 3 O 8 -H 2 O system at 800 °C and 300 MPa, at the Advanced Photon Source synchrotron facility. The analysis of the X-ray diffraction data was augmented with classical Molecular Dynamics simulations of a hydrous albite melt, incorporating accurate water-based interactions. The results show that metal–oxygen bond breaking at the bridging sites occurs overwhelmingly at the Si site upon reaction with H 2 O, with subsequent Si–OH bond formation and negligible Al–OH formation. Furthermore, we see no evidence for the dissociation of the Al 3+ ion from the network structure upon breaking of the Si–O bond in the hydrous albite melt. The results also indicate that the Na + ion is an active participant in the modifications of the silicate network structure of the albite melt upon water dissolution at high P–T conditions. We do not find evidence for the Na + ion dissociating from the network structure upon depolymerization and subsequent formation of NaOH complexes. Instead, our results show that the Na + ion persists as a structure modifier with a shift away from Na–BO bonding to an increase in the extent of Na-NBO bonding, in parallel with pronounced depolymerization of the network. Our MD simulations show that the Si–O and Al–O bond lengths are expanded by about 6% in the hydrous albite melt compared to those of the dry melt at high P–T conditions. The changes in the network silicate structure of a hydrous albite melt at high pressure and temperature, as revealed in this study, must be considered in the advancement of water dissolution models of hydrous granitic (or alkali aluminosilicate) melts.

36 MATERIALS SCIENCE↗
DOE OSTI2022

First-principles equation of state of CHON resin for inertial confinement fusion applications

A wide-range (0 to 1044.0 g/ cm 3 and 0 to 10 9 K) equation-of-state (EOS) table for a CH 1.72 O 0.37 N 0.086 quaternary compound has been constructed based on density-functional theory (DFT) molecular-dynamics (MD) calculations using a combination of Kohn-Sham DFT MD, orbital-free DFT MD, and numerical extrapolation. The first-principles EOS data are compared with predictions of simple models, including the fully ionized ideal gas and the Fermi-degenerate electron gas models, to chart their temperature-density conditions of applicability. The shock Hugoniot, thermodynamic properties, and bulk sound velocities are predicted based on the EOS table and compared to those of C-H compounds. The Hugoniot results show the maximum compression ratio of the C-H-O-N resin is larger than that of CH polystyrene due to the existence of oxygen and nitrogen; while the other properties are similar between CHON and CH. Radiation hydrodynamic simulations have been performed using the table for inertial confinement fusion targets with a CHON ablator and compared with a similar design with CH. Here, the simulations show CHON outperforms CH as the ablator for laser-direct-drive target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗