Multiscale modeling of stochastic dynamics processes with MBN Explorer.

Stochastic dynamics describes processes in complex systems having the probabilistic nature. They can involve very different dynamical systems and occur on very different temporal and spatial scale. This paper discusses the concept of stochastic dynamics and its implementation in the popular program MBN Explorer. Stochastic dynamics in MBN Explorer relies on the Monte Carlo approach and permits simulations of physical, chemical, and biological processes. The paper presents the basic theoretical concepts underlying stochastic dynamics implementation and provides several examples highlighting its applicability to different systems, such as diffusing proteins seeking an anchor point on a cell membrane, deposition of nanoparticles on a surface leading to structures with fractal morphologies, and oscillations of compounds in an autocatalytic reaction. The chosen examples illustrate the diversity of applications that can be modeled by means of stochastic dynamics with MBN Explorer.

Atomistic simulation of the FEBID-driven growth of iron-based nanostructures.

The growth of iron-containing nanostructures in the process of focused electron beam-induced deposition (FEBID) of Fe(CO)5 is studied by means of atomistic irradiation-driven molecular dynamics (IDMD) simulations. The geometrical characteristics (lateral size, height and volume), morphology and metal content of the grown nanostructures are analyzed at different irradiation and precursor replenishment conditions corresponding to the electron-limited and precursor-limited regimes (ELR & PLR) of FEBID. A significant variation of the deposit's morphology and elemental composition is observed with increasing the electron current from 1 to 4 nA. At low beam current (1 nA) corresponding to the ELR and a low degree of Fe(CO)5 fragmentation, the nanogranular structures are formed which consist of isolated iron clusters embedded into an organic matrix. In this regime, metal clusters do not coalesce with increasing electron fluence, resulting in relatively low metal content of the nanostructures. A higher beam current of 4 nA corresponding to the PLR facilitates the precursor fragmentation and the coalescence of metal clusters into a dendrite-like structure with the size corresponding to the primary electron beam. The IDMD simulations enable atomistic-level predictions on the nanoscopic characterization of the initial phase of nanostructure growth in the FEBID process. These predictions can be verified in high-resolution transmission electron microscopy experiments.

Irradiation-driven molecular dynamics simulation of the FEBID process for Pt(PF3)4.

This paper presents a detailed computational protocol for the atomistic simulation of formation and growth of metal-containing nanostructures during focused electron beam-induced deposition (FEBID). The protocol is based upon irradiation-driven molecular dynamics (IDMD), a novel and general methodology for computer simulations of irradiation-driven transformations of complex molecular systems by means of the advanced software packages MBN Explorer and MBN Studio. Atomistic simulations performed following the formulated protocol provide valuable insights into the fundamental mechanisms of electron-induced precursor fragmentation and the related mechanism of nanostructure formation and growth using FEBID, which are essential for the further advancement of FEBID-based nanofabrication. The developed computational methodology is general and applicable to different precursor molecules, substrate types, and irradiation regimes. The methodology can also be adjusted to simulate the nanostructure formation by other nanofabrication techniques using electron beams, such as direct electron beam lithography. In the present study, the methodology is applied to the IDMD simulation of the FEBID of Pt(PF3)4, a widely studied precursor molecule, on a SiO2 surface. The simulations reveal the processes driving the initial phase of nanostructure formation during FEBID, including the nucleation of Pt atoms and the formation of small metal clusters on the surface, followed by their aggregation and the formation of dendritic platinum nanostructures. The analysis of the simulation results provides spatially resolved relative metal content, height, and growth rate of the deposits, which represents valuable reference data for the experimental characterization of the nanostructures grown by FEBID.

Multiscale simulation of the focused electron beam induced deposition process.

Focused electron beam induced deposition (FEBID) is a powerful technique for 3D-printing of complex nanodevices. However, for resolutions below 10 nm, it struggles to control size, morphology and composition of the structures, due to a lack of molecular-level understanding of the underlying irradiation-driven chemistry (IDC). Computational modeling is a tool to comprehend and further optimize FEBID-related technologies. Here we utilize a novel multiscale methodology which couples Monte Carlo simulations for radiation transport with irradiation-driven molecular dynamics for simulating IDC with atomistic resolution. Through an in depth analysis of [Formula: see text] deposition on [Formula: see text] and its subsequent irradiation with electrons, we provide a comprehensive description of the FEBID process and its intrinsic operation. Our analysis reveals that simulations deliver unprecedented results in modeling the FEBID process, demonstrating an excellent agreement with available experimental data of the simulated nanomaterial composition, microstructure and growth rate as a function of the primary beam parameters. The generality of the methodology provides a powerful tool to study versatile problems where IDC and multiscale phenomena play an essential role.

Modeling MesoBioNano systems with MBN Studio made easy.

This paper introduces MesoBioNano (MBN) Studio - a graphical user interface for a popular multiscale simulation package MBN Explorer. MBN Studio has been developed to facilitate setting up and starting MBN Explorer calculations, monitoring their progress and examining the calculation results. It is tailored for any calculations that are supported by MBN Explorer, such as for example the single-point energy calculations, structure optimization, molecular dynamics, and kinetic Monte Carlo simulations. Apart from that MBN Studio has built-in tools allowing the calculation and analysis of specific characteristics that are determined by the output of the simulations, such as the diffusion coefficients of molecular species, melting temperatures and associated heat capacities, radial distribution function; a dedicated modeling plug-in allows constructing molecular systems in a quick and efficient manner. Employing this plug-in, one can easily construct molecular systems of different geometries (e.g., spherical or ellipsoidal nanoparticles, cubic crystalline samples) with various atomic composition. The paper presents the first public release of MBN Studio and provides an overview of its significant capabilities, as well as the reference point for further extensions.

Molecular dynamics simulation of self-diffusion processes in titanium in bulk material, on grain junctions and on surface.

The process of self-diffusion of titanium atoms in a bulk material, on grain junctions and on surface is explored numerically in a broad temperature range by means of classical molecular dynamics simulation. The analysis is carried out for a nanoscale cylindrical sample consisting of three adjacent sectors and various junctions between nanocrystals. The calculated diffusion coefficient varies by several orders of magnitude for different regions of the sample. The calculated values of the bulk diffusion coefficient correspond reasonably well to the experimental data obtained for solid and molten states of titanium. Investigation of diffusion in the nanocrystalline titanium is of a significant importance because of its numerous technological applications. This paper aims to reduce the lack of data on diffusion in titanium and describe the processes occurring in bulk, at different interfaces and on surface of the crystalline titanium.

Validation of classical force fields for the description of thermo-mechanical properties of transition metal materials.

It is demonstrated that classical force fields validated through the density functional theory (DFT) calculations of small titanium and nickel clusters can be applied for the description of thermo-mechanical properties of corresponding materials. This has been achieved by means of full-atom molecular dynamics simulations of nanoindentation of amorphous and nanostructured Ti and Ni-Ti materials. The theoretical analysis performed and comparison with experimental data demonstrate that the utilized classical force fields for Ti-Ti, Ni-Ni and Ni-Ti interactions describe reasonably well hardness and the Young's modulus of these materials. This observation is of the general nature and can be utilized for similar numerical exploration of thermo-mechanical properties of a broad range of materials.