Filtern
Erscheinungsjahr
Dokumenttyp
- Zeitschriftenartikel (11)
- Beitrag zu einem Tagungsband (1)
- Vortrag (1)
Sprache
- Englisch (13)
Schlagworte
- Molecular dynamics (13) (entfernen)
Organisationseinheit der BAM
Eingeladener Vortrag
- nein (1)
In the past two decades, numerous relaxation or physical aging experiments of metallic glasses have revealed signatures of intermittent atomic-scale processes. Revealed via intensity cross-correlations from coherent scattering using X-ray photon correlation spectroscopy (XPCS), the observed abrupt changes in the time-domain of atomic motion does not fit the picture of gradual slowing down of relaxation times and their origin continues to remain unclear. Using a binary Lennard-Jones model glass subjected to microsecond-long isotherms, we show here that temporally and spatially heterogeneous atomic-cluster activity at different length-scales drive the emergence of highly non-monotonous intensity cross-correlations. The simulated XPCS experiments reveal a variety of time-dependent intensity-cross correlations that, depending on both the structural evolution and the 𝑞-space sampling, give detailed insights into the possible structural origins of intermittent aging measured with XPCS.
Rapid solidification leads to unique microstructural features, where a less studied topic is the formation of various crystalline defects, including high dislocation densities, as well as gradients and splitting of the crystalline orientation. As these defects critically affect the material’s mechanical properties and performance features, it is important to understand the defect formation mechanisms, and how they depend on the solidification conditions and alloying. To illuminate the formation mechanisms of the rapid solidification induced crystalline defects, we conduct a multiscale modelling analysis consisting of bond-order potential-based molecular dynamics (MD), phase field crystal-based amplitude expansion simulations, and sequentially coupled phase field–crystal plasticity simulations. The resulting dislocation densities are quantified and compared to past experiments. The atomistic approaches (MD, PFC) can be used to calibrate continuum level crystal plasticity models, and the framework adds mechanistic insights arising from the multiscale analysis.
This article is part of the theme issue ‘Transport phenomena in complex systems (part 2)’.
Atomistic deformation simulations in the nominally elastic regime are performed for a model binary glass with strain rates as low as 10 4 /s (corresponding to 0.01 shear strain per 1 μs). A strain rate dependent elastic softening due to a micro-plasticity is observed, which is mediated by thermally-activated localized structural transformations (LSEs). A closer inspection of the atomic-scale structure indicates the material response is distinctly different for two types of local atomic environments. A system spanning iscosahe- drally coordinated substructure responds purely elastically, whereas the remaining substructure admits both elastic and microplastic evolution. This leads to a heterogeneous internal stress distribution which, upon unloading, results in negative creep and complete residual-strain recovery. A detailed structural analysis in terms of local stress, atomic displacement, and SU(2) local bonding topology shows such mi- croscopic processes can result in large changes in local stress and are more likely to occur in geomet- rically frustrated regions characterized by higher free volume and softer elastic stiffness. The thermally- activated LSE activity also mediates structural relaxation, and in this way should be distinguished from stress-driven shear transformation activity which only rejuvenates glass structure. The frequency of LSE activity, and therefore the amount of micro-plasticity, is found to be related to the degree to which the glassy state is relaxed. These insights shed atomistic light onto the structural origins that may govern re- cent experimental observations of significant structural evolution in response to elastic loading protocols.
The sliding behavior of an amorphous silica sample between two rigid surfaces is in the focus of the present paper. Molecular Dynamics using a classical Tersoff’s potential and a recently developed ReaxFF potential was applied for simulating sliding within a thin film corresponding to a tribofilm formed from silica nanoparticles. The simulations were performed at different temperatures corresponding to moderate and severe tribological stressing conditions. Simulations with both potentials revealed the need of considering different temperatures in order to obtain a sound interpretation of experimental findings. The results show the striking differences between the two potentials not only in terms of magnitude of the resistance stress (about one order of magnitude) but also in terms of friction mechanisms. The expected smooth sliding regime under high temperature conditions was predicted by both simulations, although with Tersoff’s potential smooth sliding was obtained only at the highest temperature. On the other hand, at room temperature Tersoff-style calculations demonstrate stick-slip behavior, which corresponds qualitatively with our experimental findings. Nevertheless, comparison with a macroscopic coefficient of friction is not possible because simulated resistance stresses do not depend on the applied normal pressure.
Simulations of concrete on the mesoscale require complex geometry models that resolve single aggregates. Since they fill up to 60-80% of the virtual specimen, efficient algorithms for their placement are necessary.
Event-driven molecular dynamics (EDMD) algorithms using growing spheres are proven to create dense polydisperse sphere packings. In this paper, a modified EDMD algorithm is presented and compared to the widely used random sequential addition algorithm. It reaches denser aggregate packings and saves computation time at high volume fractions.
When meshing the geometry, a minimal distance between the aggregates ensures undistorted elements. It is obtained by increasing the aggregate size during the placement process, which makes the packing problem more challenging. For a given aggregate volume fraction, the presented algorithm maximizes this distance and meshable mesoscale geometries for concrete specimens with more than 70% aggregate content are produced.
The method of movable cellular automata (MCA) and method of molecular dynamics (MD) were applied to simulate the friction and sliding behavior of model-tribofilms formed from a nanocomposite consisting of an epoxy matrix, 10 vol % micron-sized carbon fibers and 5 vol. % silica nanoparticles. Whereas MCA considered the tribofilm as an agglomerate of silica nanoparticles released from the composite and mixed with graphite particles, MD simulated the sliding behavior of an amorphous silica layer supported by stiff crystalline substrates on both sides. The MCA model provided reasonable quantitative results which corroborate experimental findings at moderate stressing conditions. The very low coefficient of friction observed experimentally under severe stressing conditions was not explained by this model. This could be attributed to the lack of mechanical data at the high temperature expected under these conditions. Although based on a simpler assumption of the tribofilm composition, MD-modelling could be easily applied to the expected high flash temperature and was able to predict friction reduction and smooth sliding under these conditions.
The problem of polydisperse sphere packings is applied to concrete mesoscale geometries in finite sized specimens. Realistic sphere diameter distributions are derived from concrete grading curves. An event-driven molecular dynamics simulation using growing particles is introduced. Compared to the widely used random sequential addition algorithm, it reaches denser aggregate packings and saves computation time at high volume fractions.
A minimal distance between particles strongly influences the maximum aggregate content. It is essential to obtain undistorted elements when meshing the geometry for finite element simulations. The algorithm maximizes this value and produces meshable concrete mesostructures with more than 70% aggregate content.
Ergopeptides, like ergocornine and a-ergocryptine, exist in an S- and in an R-configuration. Kinetic experiments imply that certain configurations are preferred depending on the solvent. The experimental methods are explained in this article. Furthermore, computational methods are used to understand this configurational preference. Standard quantum chemical methods can predict the favored configurations by using minimum energy calculations on the potential energy landscape. However, the explicit role of the solvent is not revealed by this type of methods. In order to better understand its influence, classical mechanical molecular simulations are applied. It appears from our research that 'folding' the ergopeptide molecules into an intermediate state (between the S- and the R-configuration) is mechanically hindered for the preferred configurations.
The molecular dynamics of liquid-crystalline polymethacrylates with biphenyl in the side group forming a highly ordered smectic E phase at low temperatures is investigated by dielectric spectroscopy (DS). Although no glass transition is found by DSC, DS detects a relaxation in the smectic E phase resembling to glassy dynamics. The temperature dependence of its relaxation rate is Arrhenius like at low temperatures, which changes to a VFT law at higher temperatures. Such a behavior is found to be characteristic for glassy dynamics in confining space. It is discussed, considering the structure of the polymers as nanobiphasic consisting of liquid-crystalline layers while the backbones fill the interlayer space. A comparison is made to polymethacrylates having phenyl benzoate as mesogen.
A novel approach for the simulation of host.guest systems by systematically scanning the host moleculefs orientations within the guest cavity is presented along with a thermodynamic strategy for determining preferential binding modes and corresponding optimal interaction energies between host and guest molecules. By way of example, the elution order of hexabromocyclododecane stereoisomers from high performance liquid chromatography separation on a permethylated β-cyclcodextrin stationary phase has been computed using classical molecular dynamics simulations with the explicit solvents water and acetonitrile. Comparison of estimated with experimental separation data reveals remarkable squared coefficients of correlation with R2=0.87 and a very high correlation RLOO2 . 0:72 using the leaveone-out cross-validation method and water as solvent. In particular, the approach presented shapes up as very robust in terms of the evaluated time range under consideration, reflecting well thermodynamic equilibria. These and further observations correlating with experimental results suggest the suitability of the underlying force fields and our multi-mode approach for the estimation of relative binding affinities for host–guest systems with unknown binding modes.