Score Dynamics: Scaling Molecular Dynamics with Picoseconds Time Steps via Conditional Diffusion Model.

We propose score dynamics (SD), a general framework for learning accelerated evolution operators with large timesteps from molecular dynamics (MD) simulations. SD is centered around scores or derivatives of the transition log-probability with respect to the dynamical degrees of freedom. The latter play the same role as force fields in MD but are used in denoising diffusion probability models to generate discrete transitions of the dynamical variables in an SD time step, which can be orders of magnitude larger than a typical MD time step.

Accurate parameterization of the kinetic energy functional.

The absence of a reliable formulation of the kinetic energy density functional has hindered the development of orbital free density functional theory. Using the data-aided learning paradigm, we propose a simple prescription to accurately model the kinetic energy density of any system. Our method relies on a dictionary of functional forms for local and nonlocal contributions, which have been proposed in the literature, and the appropriate coefficients are calculated via a linear regression framework.

Accurate parameterization of the kinetic energy functional for calculations using exact-exchange.

Electronic structure calculations based on Kohn-Sham density functional theory (KSDFT) that incorporate exact-exchange or hybrid functionals are associated with a large computational expense, a consequence of the inherent cubic scaling bottleneck and large associated prefactor, which limits the length and time scales that can be accessed. Although orbital-free density functional theory (OFDFT) calculations scale linearly with system size and are associated with a significantly smaller prefactor, they are limited by the absence of accurate density-dependent kinetic energy functionals. Therefore, the development of accurate density-dependent kinetic energy functionals is important for OFDFT calculations of large realistic systems.

Atomistic insights into metal hardening.

For millennia, humans have exploited the natural property of metals to get stronger or harden when mechanically deformed. Ultimately rooted in the motion of dislocations, mechanisms of metal hardening have remained in the cross-hairs of physical metallurgists for over a century. Here, we performed atomistic simulations at the limits of supercomputing that are sufficiently large to be statistically representative of macroscopic crystal plasticity yet fully resolved to examine the origins of metal hardening at its most fundamental level of atomic motion.

Probing the limits of metal plasticity with molecular dynamics simulations.

Ordinarily, the strength and plasticity properties of a metal are defined by dislocations-line defects in the crystal lattice whose motion results in material slippage along lattice planes. Dislocation dynamics models are usually used as mesoscale proxies for true atomistic dynamics, which are computationally expensive to perform routinely. However, atomistic simulations accurately capture every possible mechanism of material response, resolving every "jiggle and wiggle" of atomic motion, whereas dislocation dynamics models do not.

Radiofrequency ablation or stripping of large-diameter incompetent great saphenous varicose veins with C2 or C3 disease.

Objective: The objectives of this study were to compare the results of radiofrequency ablation (RFA) and stripping for large-diameter varicose target veins for the period of 1 year, based on a composite end point; to analyze the pain severity on a digital rating scale for 7 days after RFA and stripping; and to detect the factors affecting the level of postoperative pain using the cluster analysis.

Methods: This was a multicenter retrospective cohort study. Two groups, stripping ≥14 mm and RFA ≥14 mm, of 129 varicose vein disease patients underwent surgical treatment in three specialized clinics.

Path factorization approach to stochastic simulations.

The computational efficiency of stochastic simulation algorithms is notoriously limited by the kinetic trapping of the simulated trajectories within low energy basins. Here we present a new method that overcomes kinetic trapping while still preserving exact statistics of escape paths from the trapping basins. The method is based on path factorization of the evolution operator and requires no prior knowledge of the underlying energy landscape.

Singular orientations and faceted motion of dislocations in body-centered cubic crystals.

Dislocation mobility is a fundamental material property that controls strength and ductility of crystals. An important measure of dislocation mobility is its Peierls stress, i.e.

First-passage kinetic Monte Carlo method.

We present an efficient method for Monte Carlo simulations of diffusion-reaction processes. Introduced by us in a previous paper [Phys. Rev.

First-passage Monte Carlo algorithm: diffusion without all the hops.

We present a novel Monte Carlo algorithm for N diffusing finite particles that react on collisions. Using the theory of first-passage processes and time dependent Green's functions, we break the difficult N-body problem into independent single- and two-body propagations circumventing numerous diffusion hops used in standard Monte Carlo simulations. The new algorithm is exact, extremely efficient, and applicable to many important physical situations in arbitrary integer dimensions.

Dislocation multi-junctions and strain hardening.

At the microscopic scale, the strength of a crystal derives from the motion, multiplication and interaction of distinctive line defects called dislocations. First proposed theoretically in 1934 (refs 1-3) to explain low magnitudes of crystal strength observed experimentally, the existence of dislocations was confirmed two decades later. Much of the research in dislocation physics has since focused on dislocation interactions and their role in strain hardening, a common phenomenon in which continued deformation increases a crystal's strength.

Adaptive importance sampling Monte Carlo simulation of rare transition events.

We develop a general theoretical framework for the recently proposed importance sampling method for enhancing the efficiency of rare-event simulations [W. Cai, M. H.

Dynamic transitions from smooth to rough to twinning in dislocation motion.

The motion of dislocations in response to stress dictates the mechanical behaviour of materials. However, it is not yet possible to directly observe dislocation motion experimentally at the atomic level. Here, we present the first observations of the long-hypothesized kink-pair mechanism in action using atomistic simulations of dislocation motion in iron.

Anomalous dislocation multiplication in FCC metals.

Direct atomistic simulations of dislocation multiplication in fcc aluminum reveal an unexpected mechanism, in which a Frank-Read source emits dislocations with Burgers vectors different from that of the source itself. The mechanism is traced to a spontaneous nucleation of partial dislocation loops within the stacking fault. Understanding and a quantitative description of this unusual process are achieved through the development of a continuum model for dislocation nucleation based on the coarse-grained dislocation dynamics approach and a minimal amount of atomistic input.

Importance sampling of rare transition events in Markov processes.

We present an importance sampling technique for enhancing the efficiency of sampling rare transition events in Markov processes. Our approach is based on the design of an importance function by which the absolute probability of sampling a successful transition event is significantly enhanced, while preserving the relative probabilities among different successful transition paths. The method features an iterative stochastic algorithm for determining the optimal importance function.

Nodal effects in dislocation mobility.

We show that, contrary to the prevailing perception, dislocations can become more mobile by zipping together to form junctions. In a series of direct atomistic simulations, the critical stress to move a junction network in a [110] plane of bcc molybdenum is found to be always smaller ( approximately 50%) than that required to move isolated dislocations. Our data support a previously proposed hypothesis about the nature of anomalous slip in bcc transition metals, yet offer a different atomistic mechanism for conservative motion of screw dislocation networks.