Analysis of granular flow in a pebble-bed nuclear reactor

Beyond Local Structures in Critical Supercooled Water through Unsupervised Learning

The presence of a second critical point in water has been a topic of intense investigation for the last few decades. The molecular origins underlying this phenomenon are typically rationalized in terms of the competition between local high-density (HD) and low-density (LD) structures. Their identification often requires designing parameters that are subject to human intervention. Herein, we use unsupervised learning to discover structures in atomistic simulations of water close to the liquid-liquid critical point (LLCP). Encoding the information on the environment using local descriptors, we do not find evidence for two distinct thermodynamic structures. In contrast, when we deploy nonlocal descriptors that probe instead heterogeneities on the nanometer length scale, this leads to the emergence of LD and HD domains rationalizing the microscopic origins of the density fluctuations close to criticality.

Self-assembled soft alloy with Frank-Kasper phases beyond metals

Soft building blocks, such as micelles, cells or soap bubbles, tend to adopt near-spherical geometry when densely packed together. As a result, their packing structures do not extend beyond those discovered in metallic glasses, quasicrystals and crystals. Here we report the emergence of two Frank-Kasper phases from the self-assembly of five-fold symmetric molecular pentagons. The μ phase, an important intermediate in superalloys, is indexed in soft matter, whereas the ϕ phase exhibits a structure distinct from known Frank-Kasper phases in metallic systems. We find a broad size and shape distribution of self-assembled mesoatoms formed by molecular pentagons while approaching equilibrium that contribute to the unique packing structures. This work provides insight into the manipulation of soft building blocks that deviate from the typical spherical geometry and opens avenues for the fabrication of 'soft alloy' structures that were previously unattainable in metal alloys.

Micromechanical origin of heat transfer to granular flow

Heat transfer across a granular flow is comprised of two resistances in series : near the wall and within the bulk particle bed, neither of which is well understood due to the lack of experimental probes to separate their respective contribution. Here, we use a frequency modulated photothermal technique to separately quantify the thermal resistances in the near-wall and the bulk bed regions of particles in flowing states. Compared to the stationary state, the flowing leads to a higher near-wall resistance and a lower thermal conductivity of bulk beds. Coupled with discrete element method simulation, we show that the near-wall resistance can be explained by particle diffusion in granular flows.

Quantitative Understanding of the Onset of Dense Granular Flows

The question of when and how dense granular materials start to flow under stress, despite many industrial and geophysical applications, remains largely unresolved. We develop and test a simple equation for the onset of quasistatic flows of granular materials which is based on the frictional aging of the granular packing. The result is a nonmonotonic stress-strain relation which-akin to classical friction-is independent of the shear rate. This relation suffices to understand the quasistatic deformations of aging granular media, and its solid-to-liquid transition. Our results also elucidate the (flow) history dependence of the mechanical properties, and the sensitivity to the initial preparation of granular media.

Study on the effect of the size irregularity gradient of metal foams on macroscopic compressive properties

Size irregularity gradient and cell wall gradient, combined as the density gradient in previous studies, affect the macroscopic mechanical properties of the gradient metal foam. More and more complex mesostructures are designed and applied in metal foams, and the density gradient becomes insufficient to describe the difference in mesostructures. To explore the effect of mesostructures carefully, this study focuses on the effect of the size irregularity gradient on the macroscopic compressive properties of metal foams. A series of metal foam models were developed using the 3D Voronoi technique. These models have the same average relative densities, the same average diameters and different size irregularity gradients. Simulation results indicated that the macroscopic mechanical properties of cell wall gradient metal foams are significantly different from those of size irregularity gradient metal foams, though these models have the same relative density gradient. To explore the effect of size irregularity gradient, a theoretical model was developed to characterize the compression process from the first cell-collapse to full condensation. Theoretical results showed a linear relationship between the nominal stress and the current relative density. These findings can provide efficient guidance for the design and applications of gradient metal foams.

Geometry-controlled phase transition in vibrated granular media

We report experiments on the dynamics of vibrated particles constrained in a two-dimensional vertical container, motivated by the following question: how to get the most out of a given external vibration to maximize internal disorder (e.g. to blend particles) and agitation (e.g. to absorb vibrations)? Granular media are analogs to classical thermodynamic systems, where the injection of energy can be achieved by shaking them: fluidization arises by tuning either the amplitude or the frequency of the oscillations. Alternatively, we explore what happens when another feature, the container geometry, is modified while keeping constant the energy injection. Our method consists in modifying the container base into a V-shape to break the symmetries of the inner particulate arrangement. The lattice contains a compact hexagonal solid-like crystalline phase coexisting with a loose amorphous fluid-like phase, at any thermal agitation. We show that both the solid-to-fluid volume fraction and the granular temperature depend not only on the external vibration but also on the number of topological defects triggered by the asymmetry of the container. The former relies on the statistics of the energy fluctuations and the latter is consistent with a two-dimensional melting transition described by the KTHNY theory.

Solute Diffusivity and Local Free Volume in Cross-Linked Polymer Network: Implication of Optimizing the Conductivity of Polymer Electrolyte

The diffusion of small molecules or ions within polymeric materials is critical for their applications, such as polymer electrolytes. Cross-linking has been one of the common strategies to modulate solute diffusivity and a polymer’s mechanical properties. However, various studies have shown different effects of cross-linking on altering the solute transports. Here, we utilized coarse-grained molecular dynamics simulation to systematically analyze the effects of cross-linking and polymer rigidity of solute diffusive behaviors. Above the glass transition temperature Tg, the solute diffusion followed the Vogel–Tammann–Fulcher (VTF) equation, D = D0 e−Ea/R(T−T0). Other than the conventional compensation relation between the activation energy Ea and the pre-exponential factor D0, we also identified a correlation between Ea and Vogel temperature T0. We further characterized an empirical relation between T0 and cross-linking density. Integrating the newly identified correlations among the VTF parameters, we formulated a relation between solute diffusion and the cross-linking density. The combined results proposed the criteria for the optimal solute diffusivity in cross-linked polymers, providing generic guidance for novel polymer electrolyte design.

Mild stratification in drying films of colloidal mixtures

Size stratification of bidisperse colloidal mixtures during vertical drying was investigated using the implicit solvent Langevin dynamics (LD) simulation and the explicit solvent lattice Boltzmann (LB) method. Simulations were performed for the Péclet number (Pe) over a wide range of 1-1000. In the case of a low size ratio of 2, mild stratification was observed in both simulation methods, in contrast to distinct stratification with thick "small-on-top" or "large-on-top" layers. The LD simulations exhibited a "small-on-top" stratification or mixed state. In contrast, the LB simulations exhibited a "large-on-top" or mixed state, according to the variation in Pe. The results demonstrated that the explicit solvent reduced the collective diffusion under moderate Pe conditions. This suppressed the steep concentration gradient of small particles in the packed region of particles near the air-solvent interface. Thus, distinguishable stratification patterns were obtained for the implicit and explicit solvent models.

Molecular Dynamics Simulations of Xe Behaviors at the Grain Boundary in UO2

In this study, we investigated the behavior of xenon (Xe) bubbles in uranium dioxide (UO2) grain boundaries using molecular dynamics simulations and compared it to that in the UO2 bulk. The results show that the formation energy of Xe clusters at the Σ5 grain boundaries (GBs) is much lower than in the bulk. The diffusion activation energy of a single interstitial Xe atom at the GBs was approximately 1 eV lower than that in the bulk. Furthermore, the nucleation and growth of Xe bubbles in the Σ5 GBs at 1000 and 2000 K were simulated. The volume and pressure of bubbles with different numbers of Xe atoms were simulated. The bubble pressure dropped with increasing temperature at low Xe concentrations, whereas the volume increased. The radial distribution function was computed to explore the configuration evolution of Xe bubbles. The bubble structures in the GB and bulk material at the same temperature were also compared. Xe atoms were more regular in the bulk, whereas multiple Xe atoms formed a planar structure at the GBs.

Chrono::GPU: An Open-Source Simulation Package for Granular Dynamics Using the Discrete Element Method

We report on an open-source, publicly available C++ software module called Chrono::GPU, which uses the Discrete Element Method (DEM) to simulate large granular systems on Graphics Processing Unit (GPU) cards. The solver supports the integration of granular material with geometries defined by triangle meshes, as well as co-simulation with the multi-physics simulation engine Chrono. Chrono::GPU adopts a smooth contact formulation and implements various common contact force models, such as the Hertzian model for normal force and the Mindlin friction force model, which takes into account the history of tangential displacement, rolling frictional torques, and cohesion. We report on the code structure and highlight its use of mixed data types for reducing the memory footprint and increasing simulation speed. We discuss several validation tests (wave propagation, rotating drum, direct shear test, crater test) that compare the simulation results against experimental data or results reported in the literature. In another benchmark test, we demonstrate linear scaling with a problem size up to the GPU memory capacity; specifically, for systems with 130 million DEM elements. The simulation infrastructure is demonstrated in conjunction with simulations of the NASA Curiosity rover, which is currently active on Mars.

Critical pore radius and transport properties of disordered hard- and overlapping-sphere models

Transport properties of porous media are intimately linked to their pore-space microstructures. We quantify geometrical and topological descriptors of the pore space of certain disordered and ordered distributions of spheres, including pore-size functions and the critical pore radius δ_{c}. We focus on models of porous media derived from maximally random jammed sphere packings, overlapping spheres, equilibrium hard spheres, quantizer sphere packings, and crystalline sphere packings. For precise estimates of the percolation thresholds, we use a strict relation of the void percolation around sphere configurations to weighted bond percolation on the corresponding Voronoi networks. We use the Newman-Ziff algorithm to determine the percolation threshold using universal properties of the cluster size distribution. The critical pore radius δ_{c} is often used as the key characteristic length scale that determines the fluid permeability k. A recent study [Torquato, Adv. Wat. Resour. 140, 103565 (2020)10.1016/j.advwatres.2020.103565] suggested for porous media with a well-connected pore space an alternative estimate of k based on the second moment of the pore size 〈δ^{2}〉, which is easier to determine than δ_{c}. Here, we compare δ_{c} to the second moment of the pore size 〈δ^{2}〉, and indeed confirm that, for all porosities and all models considered, δ_{c}^{2} is to a good approximation proportional to 〈δ^{2}〉. However, unlike 〈δ^{2}〉, the permeability estimate based on δ_{c}^{2} does not predict the correct ranking of k for our models. Thus, we confirm 〈δ^{2}〉 to be a promising candidate for convenient and reliable estimates of the fluid permeability for porous media with a well-connected pore space. Moreover, we compare the fluid permeability of our models with varying degrees of order, as measured by the τ order metric. We find that (effectively) hyperuniform models tend to have lower values of k than their nonhyperuniform counterparts. Our findings could facilitate the design of porous media with desirable transport properties via targeted pore statistics.

Twisting of a Pristine α-Fe Nanowire: From Wild Dislocation Avalanches to Mild Local Amorphization

The torsion of pristine α-Fe nanowires was studied by molecular dynamics simulations. Torsion-induced plastic deformation in pristine nanowires is divided into two regimes. Under weak torsion, plastic deformation leads to dislocation nucleation and propagation. Twisting-induced dislocations are mainly 12<111> screw dislocations in a <112>-oriented nanowire. The nucleation and propagation of these dislocations were found to form avalanches which generate the emission of energy jerks. Their probability distribution function (PDF) showed power laws with mixing between different energy exponents. The mixing stemmed from simultaneous axial and radial dislocation movements. The power-law distribution indicated strongly correlated 'wild' dislocation dynamics. At the end of this regime, the dislocation pattern was frozen, and further twisting of the nanowire did not change the dislocation pattern. Instead, it induced local amorphization at the grip points at the ends of the sample. This "melting" generated highly dampened, mild avalanches. We compared the deformation mechanisms of twinned and pristine α-Fe nanowires under torsion.

Optimized Atomic Partial Charges and Radii Defined by Radical Voronoi Tessellation of Bulk Phase Simulations

We present a novel method for the computation of well-defined optimized atomic partial charges and radii from the total electron density. Our method is based on a two-step radical Voronoi tessellation of the (possibly periodic) system and subsequent integration of the total electron density within each Voronoi cell. First, the total electron density is partitioned into the contributions of each molecule, and subsequently the electron density within each molecule is assigned to the individual atoms using a second set of atomic radii for the radical Voronoi tessellation. The radii are optimized on-the-fly to minimize the fluctuation (variance) of molecular and atomic charges. Therefore, our method is completely free of empirical parameters. As a by-product, two sets of optimized atomic radii are produced in each run, which take into account many specific properties of the system investigated. The application of an on-the-fly interpolation scheme reduces discretization noise in the Voronoi integration. The approach is particularly well suited for the calculation of partial charges in periodic bulk phase systems. We apply the method to five exemplary liquid phase simulations and show how the optimized charges can help to understand the interactions in the systems. Well-known effects such as reduced ion charges below unity in ionic liquid systems are correctly predicted without any tuning, empiricism, or rescaling. We show that the basis set dependence of our method is very small. Only the total electron density is evaluated, and thus, the approach can be combined with any electronic structure method that provides volumetric total electron densities-it is not limited to Hartree-Fock or density functional theory (DFT). We have implemented the method into our open-source software tool TRAVIS.

Isomorph Invariance of Higher-Order Structural Measures in Four Lennard-Jones Systems

In the condensed liquid phase, both single- and multicomponent Lennard-Jones (LJ) systems obey the "hidden-scale-invariance" symmetry to a good approximation. Defining an isomorph as a line of constant excess entropy in the thermodynamic phase diagram, the consequent approximate isomorph invariance of structure and dynamics in appropriate units is well documented. However, although all measures of the structure are predicted to be isomorph invariant, with few exceptions only the radial distribution function (RDF) has been investigated. This paper studies the variation along isomorphs of the nearest-neighbor geometry quantified by the occurrence of Voronoi structures, Frank-Kasper bonds, icosahedral local order, and bond-orientational order. Data are presented for the standard LJ system and for three binary LJ mixtures (Kob-Andersen, Wahnström, NiY2). We find that, while the nearest-neighbor geometry generally varies significantly throughout the phase diagram, good invariance is observed along the isomorphs. We conclude that higher-order structural correlations are no less isomorph invariant than is the RDF.

Unravelling the role of phoretic and hydrodynamic interactions in active colloidal suspensions

Active fluids comprise a variety of systems composed of elements immersed in a fluid environment which can convert some form of energy into directed motion; as such they are intrinsically out-of-equilibrium in the absence of any external force. A fundamental problem in the physics of active matter concerns the understanding of how the characteristics of autonomous propulsion and agent-agent interactions determine the collective dynamics of the system. We study numerically the suspensions of self-propelled diffusiophoretic colloids, in (quasi)-2d configurations, accounting for both dynamically resolved solute-mediated phoretic interactions and solvent-mediated hydrodynamic interactions. Our results show that the system displays different scenarios at changing the colloid-solute affinity and it develops a cluster phase in the chemoattractive case. We study the statistics of cluster sizes and cluster morphologies for different magnitudes of colloidal activity. Finally, we provide evidences that hydrodynamics plays a relevant role in the aggregation kinetics and cluster morphology, significantly hindering cluster growth.

Dynamic metal-polymer interaction for the design of chemoselective and long-lived hydrogenation catalysts

Metal catalysts are generally supported on hard inorganic materials because of their high thermochemical stabilities. Here, we support Pd catalysts on a thermochemically stable but "soft" engineering plastic, polyphenylene sulfide (PPS), for acetylene partial hydrogenation. Near the glass transition temperature (~353 K), the mobile PPS chains cover the entire surface of Pd particles via strong metal-polymer interactions. The Pd-PPS interface enables H2 activation only in the presence of acetylene that has a strong binding affinity to Pd and thus can disturb the Pd-PPS interface. Once acetylene is hydrogenated to weakly binding ethylene, re-adsorption of PPS on the Pd surface repels ethylene before it is further hydrogenated to ethane. The Pd-PPS interaction enables selective partial hydrogenation of acetylene to ethylene even in an ethylene-rich stream and suppresses catalyst deactivation due to coke formation. The results manifest the unique possibility of harnessing dynamic metal-polymer interaction for designing chemoselective and long-lived catalysts.

Particle-mesh two-dimensional pattern reverse Monte Carlo analysis on filled-gels during uniaxial expansion

A particle-mesh-based two-dimensional pattern reverse Monte Carlo (RMC) analysis method (PM-2DpRMC) is proposed for analyzing two-dimensional small-angle-scattering (2D-SAS) patterns. The results of analyzing such patterns in expanded gel networks filled with spherical nano-particles during uniaxial elongation are provided. Previously, characteristic 2D-SAS patterns, such as the two-point bar plate and figure-8 patterns, were observed to possess high stretching ratios in a coarse-grained molecular dynamics (CGMD) simulation, and the patterns were found to depend on the properties of the networks, such as whether they were topologically percolated or not. To establish real-space visualizations of the changes in the particle-configurations corresponding to the 2D-SAS patterns, the on-the-fly PM-2DpRMC method was employed to model the morphology changes of the filler particles during stretching in a series of 2D-SAS patterns. The use of quasi-dynamics to express deformation was validated by evaluating the accuracy of the elongation speed and the number of RMC trials predicted from the behavior of the converged χ2. The results confirm that the obtained sequences of 3D configurations were similar to the original configurations obtained in the CGMD simulations. Further, a demonstrative test using actual experimental data was conducted to verify the completeness of the confirmation.

Protocol Dependence and State Variables in the Force-Moment Ensemble

Stress-based ensembles incorporating temperaturelike variables have been proposed as a route to an equation of state for granular materials. To test the efficacy of this approach, we perform experiments on a two-dimensional photoelastic granular system under three loading conditions: uniaxial compression, biaxial compression, and simple shear. From the interparticle forces, we find that the distributions of the normal component of the coarse-grained force-moment tensor are exponential tailed, while the deviatoric component is Gaussian distributed. This implies that the correct stress-based statistical mechanics conserves both the force-moment tensor and the Maxwell-Cremona force-tiling area. As such, two variables of state arise: the tensorial angoricity (α[over ^]) and a new temperaturelike quantity associated with the force-tile area which we name keramicity (κ). Each quantity is observed to be inversely proportional to the global confining pressure; however, only κ exhibits the protocol independence expected of a state variable, while α[over ^] behaves as a variable of process.

Shockwave generates < 100 > dislocation loops in bcc iron

The formation mechanism of < 100 > interstitial dislocation loops in ferritic steels stemming from irradiation remains elusive, as their formations are either too short for experiments, or too long for molecular dynamics simulations. Here, we report on the formation of both interstitial and vacancy dislocation loops in high energy displacement cascades using large-scale molecular dynamics simulations with up to 220 million atoms. Riding the supersonic shockwave generated in the cascade, self-interstitial atoms are punched out to form < 100 > dislocation loops in only a few picoseconds during one single cascade event, which is several orders of magnitude faster than any existing mechanisms. The energy analysis suggests that the formation of the interstitial loops depends on kinetic energy redistribution, where higher incidence energy or larger atom mass could improve the probability of the direct nucleation of interstitial dislocation loops.

Irradiating iron introduces defects such as interstitial dislocation loops, whose exact formation mechanism remains unclear. Here, the authors use large scale molecular dynamics simulations to reveal a punch out mechanism that can directly create < 100 > interstitial dislocation loops.

Simulating structure and dynamics in small droplets of 1-ethyl-3-methylimidazolium acetate

To investigate the structure and dynamics of small ionic liquid droplets in gas phase, we performed a DFT-based ab initio molecular dynamics study of several 1-ethyl-3-methylimidazolium acetate clusters in vacuum as well as a bulk phase simulation. We introduce an unbiased criterion for average droplet diameter and density. By extrapolation of the droplet densities, we predict the experimental bulk phase density with a deviation of only a few percent. The hydrogen bond geometry between cations and anions is very similar in droplets and bulk, but the hydrogen bond dynamics is significantly slower in the droplets, becoming slower with increasing system size, with hydrogen bond lifetimes up to 2000 ps. From a normal mode analysis of the trajectories, we identify the modes of the ring proton C-H stretching, which are strongly affected by hydrogen bonding. From analyzing these, we find that the hydrogen bond becomes weaker with increasing system size. The cations possess an increased concentration inside the clusters, whereas the anions show an excess concentration on the outside. Almost all anions point towards the droplet center with their carboxylic groups. Ring stacking is found to be a very important structural motif in the droplets (as in the bulk), but side chain interactions are only of minor importance. By using Voronoi tessellation, we define the exposed droplet surface and find that it consists mainly of hydrogen atoms from the cation's and anion's methyl and ethyl groups. Polar atoms are rarely found on the surface, such that the droplets appear completely hydrophobic on the outside.

Chromatin Loop Extrusion and Chromatin Unknotting

It has been a puzzle how decondensed interphase chromosomes remain essentially unknotted. The natural expectation is that in the presence of type II DNA topoisomerases that permit passages of double-stranded DNA regions through each other, all chromosomes should reach the state of topological equilibrium. The topological equilibrium in highly crowded interphase chromosomes forming chromosome territories would result in formation of highly knotted chromatin fibres. However, Chromosome Conformation Capture (3C) methods revealed that the decay of contact probabilities with the genomic distance in interphase chromosomes is practically the same as in the crumpled globule state that is formed when long polymers condense without formation of any knots. To remove knots from highly crowded chromatin, one would need an active process that should not only provide the energy to move the system from the state of topological equilibrium but also guide topoisomerase-mediated passages in such a way that knots would be efficiently unknotted instead of making the knots even more complex. We perform coarse-grained molecular dynamics simulations of the process of chromatin loop extrusion involving knotted and catenated chromatin fibres to check whether chromatin loop extrusion may be involved in active unknotting of chromatin fibres. Our simulations show that the process of chromatin loop extrusion is ideally suited to actively unknot, decatenate and demix chromatin fibres in interphase chromosomes.

Coarse-grained modeling of multiphase interactions at microscale

The objective of this study is to develop and test a coarse-grained molecular dynamics framework to model microscale multiphase systems with different inter-particle interactions and recover emerging thermodynamic and mechanical properties at the microscale. A water-vapor model and a fused silica model are developed to demonstrate the capability of our framework. The former can reproduce the water density and surface tension over a wide range of temperatures; the latter can reproduce experimental density, tensile strength, and Young's modulus of fused silica. Therefore, the deformable solid model is implemented in the proposed framework. Validations of spatial scaling methods for solid, liquid, and multiphase systems suggest that the proposed framework can be calibrated at an arbitrary microscale and used at a different length scale without recalibration. Different values of wettability for a solid-liquid-vapor system that is characterized by the contact angle can be achieved by changing the solid-liquid inter-particle potential. Thanks to these features, the proposed coarse-grained molecular dynamics framework can potentially find applications in modeling systems in which multiple phases coexist and have substantial interactions.

Configuration correlation governs slow dynamics of supercooled metallic liquids

The origin of dramatic slowing down of dynamics in metallic glass-forming liquids toward their glass transition temperatures is a fundamental but unresolved issue. Through extensive molecular dynamics simulations, here we show that, contrary to the previous beliefs, it is not local geometrical orderings extracted from instantaneous configurations but the intrinsic correlation between configurations that captures the structural origin governing slow dynamics. More significantly, it is demonstrated by scaling analyses that it is the correlation length extracted from configuration correlation rather than dynamic correlation lengths that is the key to determine the drastic slowdown of supercooled metallic liquids. The key role of the configuration correlation established here sheds important light on the structural origin of the mysterious glass transition and provides an essential piece of the puzzle for the development of a universal theoretical understanding of glass transition in glasses.

Study of Pressure Drop in Fixed Bed Reactor Using a Computational Fluid Dynamics (CFD) Code

Pressure drops of water and critical steam flowing in the fixed bed of mono-sized spheres are studied using SolidWorks 2017 Flow Simulation CFD code. The effects of the type of bed formation, flow velocity, density, and pebble size are evaluated. A new equation is concluded from the data, which is able to estimate pressure drop of a packed bed for high particle Reynolds number, from 15,000 to 1,000,000.

The Quantified Characterization Method of the Micro-Macro Contacts of Three-Dimensional Granular Materials on the Basis of Graph Theory

We have attempted a multiscale and quantified characterization method of the contact in three-dimensional granular material made of spherical particles, particularly in cemented granular material. Particle contact is defined as a type of surface contact with voids in its surroundings, rather than a point contact. Macro contact is a particle contact set satisfying the restrictive condition of a two-dimensional manifold with a boundary. On the basis of graph theory, two dual geometrical systems are abstracted from the granular pack. The face and the face set, which satisfies the two-dimensional manifold with a boundary in the solid cell system, are extracted to characterize the particle contact and the macro contact, respectively. This characterization method is utilized to improve the post-processing in DEM (Discrete Element Method) from a micro perspective to describe the macro effect of the cemented granular material made of spherical particles. Since the crack has the same shape as its corresponding contact, this method is adopted to characterize the crack and realize its visualization. The integral failure route of the sample can be determined by a graph theory algorithm. The contact force is assigned to the weight value of the face characterizing the particle contact. Since the force vectors can be added, the macro contact force can be solved by adding the weight of its corresponding faces.

Dislocation loop formation by swift heavy ion irradiation of metals

A coupled two-temperature, molecular dynamics methodology is used to simulate the structural evolution of bcc metals (Fe and W) and fcc metals (Cu and Ni) following irradiation by swift heavy ions. Electronic temperature dependent electronic specific heat capacities and electron-phonon coupling strengths are used to capture the full effects of the variation in the electronic density of states. Tungsten is found to be significantly more resistant to damage than iron, due both to the higher melting temperature and the higher thermal conductivity. Very interesting defect structures, quite different from defects formed in cascades, are found to be created by swift heavy ion irradiation in the bcc metals. Isolated vacancies form a halo around elongated interstitial dislocation loops that are oriented along the ion path. Such configurations are formed by rapid recrystallization of the molten cylindrical region that is created by the energetic ion. Vacancies are created at the recrystallization front, resulting in excess atoms at the core which form interstitial dislocation loops on completion of crystallization. These unique defect structures could, potentially, be used to create metal films with superior mechanical properties and interesting nanostructures.