• Title/Summary/Keyword: Molecular-Dynamics Simulation

Search Result 540, Processing Time 0.027 seconds

Damage studies on irradiated tungsten by helium ions in a plasma focus device

  • Seyyedhabashy, Mir mohammadreza;Tafreshi, Mohammad Amirhamzeh;bidabadi, Babak Shirani;Shafiei, Sepideh;Nasiri, Ali
    • Nuclear Engineering and Technology
    • /
    • v.52 no.4
    • /
    • pp.827-834
    • /
    • 2020
  • Damage of tungsten due to helium ions of a PF device was studied. The tungsten was analyzed by SEM and AFM after irradiation. SEM revealed fine bubbles of helium atoms with diameters of a few nanometers, which join and form larger bubbles and blisters on the surface of tungsten. This observation confirmed the results of molecular dynamics simulation. SEM analysis after etching of the irradiated surface indicated cavities with depth range of 35-85 nm. The average fluence of helium ion of the PF device was calculated about 5.2 × 1015 cm-2 per shot, using Lee code. Energy spectrum of helium ions was estimated using a Thomson parabola spectrometer as a function of dN/dE ∝ E-2.8 in the energy range of 10-200 keV. The characteristics of helium ion beam was imported to SRIM code. SRIM revealed that the maximum DPA and maximum helium concentration occur in the depth range of 20-50 nm. SRIM also showed that at depth of 30 nm, all of the tungsten atoms are displaced after 20 shots, while at depth of higher than 85 nm the destruction is insignificant. There is a close match between SRIM results and the measured depths of cavities in SEM images of tungsten after etching.

Characterization of Thickness and Thermoelastic Properties of Interphase in Polymer Nanocomposites using Multiscale Analysis (멀티스케일 해석을 통한 고분자 나노복합재의 계면 상 두께와 열탄성 물성 도출)

  • Choi, Joonmyung;Cho, Maenghyo
    • Journal of the Computational Structural Engineering Institute of Korea
    • /
    • v.29 no.6
    • /
    • pp.577-582
    • /
    • 2016
  • In this study, a multiscale method for solving a thermoelasticity problem for interphase in the polymeric nanocomposites is developed. Molecular dynamics simulation and finite element analysis were numerically combined to describe the geometrical boundaries and the local mechanical response of the interfacial region where the polymer networks were highly interacted with the nanoparticle surface. Also, the micrmechanical thermoelasticity equations were applied to the obtained equivalent continuum unit to compute the growth of interphase thickness according to the size of nanoparticles, as well as the thermal phase transition behavior at a wide range of temperatures. Accordingly, the equivalent continuum model obtained from the multiscale analysis provides a meaningful description of the thermoelastic behavior of interphase as well as its nanoparticle size effect on thermoelasticity at both below and above the glass transition temperature.

MD simulation of structural change of polyethylene induced by high energy ion bombardment

  • Kim, Chan-Soo;Ahmed, Sk. Faruque;Moon, Myoung-Woon;Lee, Kwang-Ryeol
    • Proceedings of the Korean Vacuum Society Conference
    • /
    • 2010.02a
    • /
    • pp.358-358
    • /
    • 2010
  • Ion beam bombardment at low energy forms nanosize patterns such as ripples, dots or wrinkles on the surface of polymers in ambient temperature and pressure. It has been known that the ion beam can alter the polymer surface that induces skins stiffer or the density higher by higher compressive stress or strain energies associated with chain scissions and crosslinks of the polymer. Atomic scale structure evolution in polymers is essential to understand a stress generation mechanism during the ion beam bombardment, which governs the nanoscale surface structure evolution. In this work, Molecular Dynamics (MD) simulations are employed to characterize the phenomenon occurred in bombardment between the ion beam and polymers that forms nanosize patterns. We investigate the structure evolution of Low Density Polyethylene (LDPE) at 300 K as the polymer is bombarded with Argon ions having various kinetic energies ranging from 100 eV to 1 KeV with 50 eV intervals having the fluence of $1.45\;{\times}\;1014 #/cm2$. These simulations use the Reactive Force Field (ReaxFF), which can mimic chemical covalent bonds and includes van der Waals potentials for describing the intermolecular interactions. The results show the details of the structural evolution of LDPE by the low energy Ar ion bombardment. Analyses through kinetic and potential energy, number of crosslinks and chain scissions, level of local densification and motions of atoms support that the residual strain energies on the surface is strongly associated with the number of crosslinks or scissored chains. Also, we could find an optimal Ar ion beam energy to make crosslinks well.

  • PDF

MULTI-SCALE MODELING AND ANALYSIS OF CONVECTIVE BOILING: TOWARDS THE PREDICTION OF CHF IN ROD BUNDLES

  • Niceno, B.;Sato, Y.;Badillo, A.;Andreani, M.
    • Nuclear Engineering and Technology
    • /
    • v.42 no.6
    • /
    • pp.620-635
    • /
    • 2010
  • In this paper we describe current activities on the project Multi-Scale Modeling and Analysis of convective boiling (MSMA), conducted jointly by the Paul Scherrer Institute (PSI) and the Swiss Nuclear Utilities (Swissnuclear). The long-term aim of the MSMA project is to formulate improved closure laws for Computational Fluid Dynamics (CFD) simulations for prediction of convective boiling and eventually of the Critical Heat Flux (CHF). As boiling is controlled by the competition of numerous phenomena at various length and time scales, a multi-scale approach is employed to tackle the problem at different scales. In the MSMA project, the scales on which we focus range from the CFD scale (macro-scale), bubble size scale (meso-scale), liquid micro-layer and triple interline scale (micro-scale), and molecular scale (nano-scale). The current focus of the project is on micro- and meso-scales modeling. The numerical framework comprises a highly efficient, parallel DNS solver, the PSI-BOIL code. The code has incorporated an Immersed Boundary Method (IBM) to tackle complex geometries. For simulation of meso-scales (bubbles), we use the Constrained Interpolation Profile method: Conservative Semi-Lagrangian $2^{nd}$ order (CIP-CSL2). The phase change is described either by applying conventional jump conditions at the interface, or by using the Phase Field (PF) approach. In this work, we present selected results for flows in complex geometry using the IBM, selected bubbly flow simulations using the CIP-CSL2 method and results for phase change using the PF approach. In the subsequent stage of the project, the importance of effects of nano-scale processes on the global boiling heat transfer will be evaluated. To validate the models, more experimental information will be needed in the future, so it is expected that the MSMA project will become the seed for a long-term, combined theoretical and experimental program.

Catalytic Ability Improvement of Phenylalanine Hydroxylase from Chromobacterium violaceum by N-Terminal Truncation and Proline Introduction

  • Liu, Zhongmei;Cheng, Zhongyi;Ye, Shuangshuang;Zhou, Li;Zhou, Zhemin
    • Journal of Microbiology and Biotechnology
    • /
    • v.29 no.9
    • /
    • pp.1375-1382
    • /
    • 2019
  • Phenylalanine hydroxylase from Chromobacterium violaceum (CvPAH) is a monomeric enzyme that converts phenylalanine to tyrosine. It shares high amino acid identity and similar structure with a subunit of human phenylalanine hydroxylase that is a tetramer, resulting in the latent application in medications. In this study, semirational design was applied to CvPAH to improve the catalytic ability based on molecular dynamics simulation analyses. Four N-terminal truncated variants and one single point variant were constructed and characterized. The D267P variant showed a 2.1-fold increased thermal stability compared to the wild type, but lower specific activity was noted compared with the wild type. The specific activity of all truncated variants was a greater than 25% increase compared to the wild type, and these variants showed similar or slightly decreased thermostability with the exception of the $N-{\Delta}9$ variant. Notably, the $N-{\Delta}9$ variant exhibited a 1.2-fold increased specific activity, a 1.3-fold increased thermostability and considerably increased catalytic activity under the neutral environment compared with the wild type. These properties of the $N-{\Delta}9$ variant could advance medical and pharmaceutical applications of CvPAH. Our findings indicate that the N-terminus might modulate substrate binding, and are directives for further modification and functional research of PAH and other enzymes.

Large cylindrical deflection analysis of FG carbon nanotube-reinforced plates in thermal environment using a simple integral HSDT

  • Djilali, Nassira;Bousahla, Abdelmoumen Anis;Kaci, Abdelhakim;Selim, Mahmoud M.;Bourada, Fouad;Tounsi, Abdeldjebbar;Tounsi, Abdelouahed;Benrahou, Kouider Halim;Mahmoud, S.R.
    • Steel and Composite Structures
    • /
    • v.42 no.6
    • /
    • pp.779-789
    • /
    • 2022
  • This work presents a non-linear cylindrical bending analysis of functionally graded plate reinforced by single-walled carbon nanotubes (SWCNTs) in thermal environment using a simple integral higher-order shear deformation theory (HSDT). This theory does not require shear correction factors and the transverse shear stresses vary parabolically through the thickness. The material properties of SWCNTs are assumed to be temperature-dependent and are obtained from molecular dynamics simulations. The material properties of functionally graded carbon nanotube-reinforced composites (FG-CNTCRs) are considered to be graded in the thickness direction, and are estimated through a micromechanical model. The non-linear strain-displacement relations in the Von Karman sense are used to study the effect of geometric non-linearity and the solution is obtained by minimization of the total potential energy. The numerical illustrations concern the nonlinear bending response of FG-CNTRC plates under different sets of thermal environmental conditions, from which results for uniformly distributed CNTRC plates are obtained as benchmarks.

Nucleotide Sequence, Structural Investigation and Homology Modeling Studies of a Ca2+-independent α-amylase with Acidic pH-profile

  • Sajedi, Reza Hassan;Taghdir, Majid;Naderi-Manesh, Hossein;Khajeh, Khosro;Ranjbar, Bijan
    • BMB Reports
    • /
    • v.40 no.3
    • /
    • pp.315-324
    • /
    • 2007
  • The novel $\alpha$-amylase purified from locally isolated strain, Bacillus sp. KR-8104, (KRA) (Enzyme Microb Technol; 2005; 36: 666-671) is active in a wide range of pH. The enzyme maximum activity is at pH 4.0 and it retains 90% of activity at pH 3.5. The irreversible thermoinactivation patterns of KRA and the enzyme activity are not changed in the presence and absence of $Ca^{2+}$ and EDTA. Therefore, KRA acts as a $Ca^{2+}$-independent enzyme. Based on circular dichroism (CD) data from thermal unfolding of the enzyme recorded at 222 nm, addition of $Ca^{2+}$ and EDTA similar to its irreversible thermoinactivation, does not influence the thermal denaturation of the enzyme and its Tm. The amino acid sequence of KRA was obtained from the nucleotide sequencing of PCR products of encoding gene. The deduced amino acid sequence of the enzyme revealed a very high sequence homology to Bacillus amyloliquefaciens (BAA) (85% identity, 90% similarity) and Bacillus licheniformis $\alpha$-amylases (BLA) (81% identity, 88% similarity). To elucidate and understand these characteristics of the $\alpha$-amylase, a model of 3D structure of KRA was constructed using the crystal structure of the mutant of BLA as the platform and refined with a molecular dynamics (MD) simulation program. Interestingly enough, there is only one amino acid substitution for KRA in comparison with BLA and BAA in the region involved in the calcium-binding sites. On the other hand, there are many amino acid differences between BLA and KRA at the interface of A and B domains and around the metal triad and active site area. These alterations could have a role in stabilizing the native structure of the loop in the active site cleft and maintenance and stabilization of the putative metal triad-binding site. The amino acid differences at the active site cleft and around the catalytic residues might affect their pKa values and consequently shift its pH profile. In addition, the intrinsic fluorescence intensity of the enzyme at 350 nm does not show considerable change at pH 3.5-7.0.

Functional Defects of Hb Kempsey (${\beta}99Asp{\rightarrow}Asn$) Can be Compensated by Insertion of a New Intersubunit Hydrogen Bond at the ${\alpha}_1{\beta}_2$ Subunit Interface

  • Yeh, Byung-Il;Choi, Jong-Whan;Sohn, Joon-Hyung;Lee, Hyean-Woo;Han, Dong-Pyou;Jung, Seun-Ho;Kim, Hyun-Won
    • BMB Reports
    • /
    • v.31 no.6
    • /
    • pp.590-594
    • /
    • 1998
  • X-ray crystallographic studies of the deoxy form of human adult hemoglobin (Hb A) have shown that ${\beta}99Asp$ is hydrogen bonded to both ${\alpha}42Tyr$ and ${\alpha}97Asn$ in the ${\alpha}_1{\beta}_2$ subunit interface, suggesting that the essential role of ${\beta}99Asp$ is to stabilize the deoxy-Hb by creating the intersubunit hydrogen bond. In particular, for Hb Kempsey (${\beta}99Asp{\rightarrow}Asn$), molecular dynamics simulation indicated that a new hydrogen bond involving ${\beta}99Asn$ can be induced by replacing ${\alpha}42Tyr$ with a strong hydrogen-bond acceptor such as Asp. Designed mutant recombinant (r) Hb (${\beta}99Asp{\rightarrow}Asn$, ${\alpha}42Tyr{\rightarrow}Asp$) have been produced in the Escherichia coli expression system and have shown that functional defects of Hb Kempsey could be compensated by the ${\alpha}42Tyr{\rightarrow}Asp$ substitution. However, as the ${\alpha}42 Tyr{\rightarrow}Asp$ mutation has never been reported before, it is still possible that the functional properties of r Hb (${\beta}99Asp{\rightarrow}Asn$, ${\alpha}42Tyr{\rightarrow}Asp$) may be due to the mutation itself. Thus, it is required to produce r Hb (${\alpha}42Tyr{\rightarrow}Asp$) and r Hb Kempsey (${\beta}99Asp{\rightarrow}AsnX$( as controls, and to compare their properties with those of r Hb (${\beta}99Asp{\rightarrow}Asn$, ${\alpha}42Tyr{\rightarrow}Asp$). r Hb (${\alpha}42Tyr{\rightarrow}Asp$) could not be purified because it is an unstable hemoglobin which forms Heinz bodies. r Hb Kempsey (${\beta}99Asp{\rightarrow}Asn$) exhibits very high oxygen affinity and greatly reduced cooperativity. Thus, r Hb (${\beta}99Asp{\rightarrow}Asn$) and r Hb (${\alpha}42Tyr{\rightarrow}Asp)$ compensate each other.

  • PDF

Design Sensitivity Analysis of Coupled MD-Continuum Systems Using Bridging Scale Approach (브리징 스케일 기법을 이용한 분자동역학-연속체 연성 시스템의 설계민감도 해석)

  • Cha, Song-Hyun;Ha, Seung-Hyun;Cho, Seonho
    • Journal of the Computational Structural Engineering Institute of Korea
    • /
    • v.27 no.3
    • /
    • pp.137-145
    • /
    • 2014
  • We present a design sensitivity analysis(DSA) method for multiscale problems based on bridging scale decomposition. In this paper, we utilize a bridging scale method for the coupled system analysis. Since the analysis of full MD systems requires huge amount of computational costs, a coupled system of MD-level and continuum-level simulation is usually preferred. The information exchange between the MD and continuum levels is taken place at the MD-continuum boundary. In the bridging scale method, a generalized Langevin equation(GLE) is introduced for the reduced MD system and the GLE force using a time history kernel is applied at the boundary atoms in the MD system. Therefore, we can separately analyze the MD and continuum level simulations, which can accelerate the computing process. Once the simulation of coupled problems is successful, the need for the DSA is naturally arising for the optimization of macro-scale design, where the macro scale performance of the system is maximized considering the micro scale effects. The finite difference sensitivity is impractical for the gradient based optimization of large scale problems due to the restriction of computing costs but the analytical sensitivity for the coupled system is always accurate. In this study, we derive the analytical design sensitivity to verify the accuracy and applicability to the design optimization of the coupled system.

Multiscale Analysis on Expectation of Mechanical Behavior of Polymer Nanocomposites using Nanoparticulate Agglomeration Density Index (나노 입자의 군집밀도를 이용한 고분자 나노복합재의 기계적 거동 예측에 대한 멀티스케일 연구)

  • Baek, Kyungmin;Shin, Hyunseong;Han, Jin-Gyu;Cho, Maenghyo
    • Composites Research
    • /
    • v.30 no.5
    • /
    • pp.323-330
    • /
    • 2017
  • In this study, multiscale analysis in which the information obtained from molecular dynamics simulation is applied to the continuum mechanics level is conducted to investigate the effects of clustering of silicon carbide nanoparticles reinforced into polypropylene matrix on mechanical behavior of nanocomposites. The elastic behavior of polymer nanocomposites is observed for various states of nanoparticulate agglomeration according to the model reflecting the degradation of interphase properties. In addition, factors which mainly affect the mechanical behavior of the nanocomposites are identified, and new index 'clustering density' is defined. The correlation between the clustering density and the elastic modulus of nanocomposites is understood. As the clustering density increases, the interfacial effect decreased and finally the improvement of mechanical properties is suppressed. By considering the random distribution of the nanoparticles, the range of elastic modulus of nanocomposites for same value of clustering density can be investigated. The correlation can be expressed in the form of exponential function, and the mechanical behavior of the polymer nanocomposites can be effectively predicted by using the nanoparticulate clustering density.