Source author record

Mihai-Cosmin Marinica

Mihai-Cosmin Marinica appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.stat-mech cond-mat.mtrl-sci physics.chem-ph math-ph math.MP

Catalog footprint

What is connected

6works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2018arXiv

Relaxation volumes of microscopic and mesoscopic irradiation-induced defects in tungsten

The low energy structures of irradiation-induced defects have been studied in detail, as these determine the available modes by which a defect can diffuse or relax. As a result, there are many studies concerning the relative energies of possible defect structures, and empirical potentials are commonly fitted to or evaluated with respect to these energies. But recently [Dudarev et al Nuclear Fusion 2018], we have shown how to determine the stresses, strains and swelling of reactor components under irradiation from the elastic properties of ensembles of irradiation-induced defects. These elastic properties have received comparatively little attention. Here we evaluate relaxation volumes of irradiation-induced defects in tungsten computed with empirical potentials, and compare to density functional theory results where available. Different empirical potentials give different results, but some potential-independent trends in relaxation volumes can be identified. We show that the relaxation volume of small defects can be predicted to within 10% from their point-defect count. For larger defects we provide empirical fits for the relaxation volume of as a function of size. We demonstrate that the relaxation volume associated with a single primary-damage cascade can be estimated from the primary knock-on atom (PKA) energy. We conclude that while annihilation of vacancy- and interstitial- character defects will invariably reduce the total relaxation volume of the cascade debris, empirical potentials disagree whether coalescence of defects will reduce or increase the total relaxation volume.

preprint2015arXiv

Using Bayes formula to estimate rates of rare events in transition path sampling simulations

Transition path sampling is a method for estimating the rates of rare events in molecular systems based on the gradual transformation of a path distribution containing a small fraction of reactive trajectories into a biased distribution in which these rare trajectories have become frequent. Then, a multistate reweighting scheme is implemented to postprocess data collected from the staged simulations. Herein, we show how Bayes formula allows to directly construct a biased sample containing an enhanced fraction of reactive trajectories and to concomitantly estimate the transition rate from this sample. The approach can remediate the convergence issues encountered in free energy perturbation or umbrella sampling simulations when the transformed distribution insufficiently overlaps with the reference distribution.

preprint2014arXiv

Free energy calculations from adaptive molecular dynamics simulations with adiabatic reweighting

We propose an adiabatic reweighting algorithm for computing the free energy along an external parameter from adaptive molecular dynamics simulations. The adaptive bias is estimated using Bayes identity and information from all the sampled configurations. We apply the algorithm to a structural transition in a cluster and to the migration of a crystalline defect along a reaction coordinate. Compared to standard adaptive molecular dynamics, we observe an acceleration of convergence. With the aid of the algorithm, it is also possible to iteratively construct the free energy along the reaction coordinate without having to differentiate the gradient of the reaction coordinate or any biasing potential.

preprint2013arXiv

Point defect modeling in materials: coupling ab initio and elasticity approaches

Modeling point defects at an atomic scale requires careful treatment of the long-range atomic relaxations. This elastic field can strongly affect point defect properties calculated in atomistic simulations because of the finite size of the system under study. This is an important restriction for ab initio methods which are limited to a few hundred atoms. We propose an original approach coupling ab initio calculations and linear elasticity theory to obtain the properties of an isolated point defect for reduced supercell sizes. The reliability and benefit of our approach are demonstrated for three problematic cases: the self-interstitial in zirconium, clusters of self-interstitials in iron, and the neutral vacancy in silicon.

preprint2012arXiv

Estimating time-correlation functions by sampling and unbiasing dynamically activated events

Transition path sampling is a rare-event method that estimates state-to-state timecorrelation functions in many-body systems from samples of short trajectories. In this framework, it is proposed to bias the importance function using the lowest Jacobian eigenvalue moduli along the dynamical trajectory. A lowest eigenvalue modulus is related to the lowest eigenvalue of the Hessian matrix and is evaluated here using the Lanczos algorithm as in activation-relaxation techniques. This results in favoring the sampling of activated trajectories and enhancing the occurrence of the rare reactive trajectories of interest, those corresponding to transitions between locally stable states. Estimating the time-correlation functions involves unbiasing the sample of simulated trajectories which is done using the multi-state Bennett acceptance ratio (MBAR) method. To assess the performance of our procedure, we compute the time-correlation function associated with the migration of a vacancy in α-iron. The derivative of the estimated time-correlation function yields a migration rate in agreement with the one given by transition state theory. Besides, we show that the information relative to rejected trajectories can be recycled within MBAR, resulting in a substantial speed-up. Unlike original transition path-sampling, our approach does not require computing the reversible work to confine the trajectory endpoints to a reactive state.

preprint2010arXiv

Free energy reconstruction from steered dynamics without post-processing

Various methods achieving importance sampling in ensembles of nonequilibrium trajectories enable to estimate free energy differences and, by maximum-likelihood post-processing, to reconstruct free energy landscapes. Here, based on Bayes theorem, we propose a more direct method in which a posterior likelihood function is used both to construct the steered dynamics and to infer the contribution to equilibrium of all the sampled states. The method is implemented with two steering schedules. First, using non-autonomous steering, we calculate the migration barrier of the vacancy in Fe-alpha. Second, using an autonomous scheduling related to metadynamics and equivalent to temperature-accelerated molecular dynamics, we accurately reconstruct the two-dimensional free energy landscape of the 38-atom Lennard-Jones cluster as a function of an orientational bond-order parameter and energy, down to the solid-solid structural transition temperature of the cluster and without maximum-likelihood post-processing.